Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 390,641 C→T 15.9% F411F (TTC→TTT proY → predicted cryptic proline transporter

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606390,6410C→T15.9% 2.0 ‑2.5 ‑0.0361233proYpredicted cryptic proline transporter
Reads supporting (aligned to +/- strand):  new base (6/0):  ref base (18/12):  total (24/12)
Fisher's exact test strand distribution p-value = 7.93e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.89e-01

GGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCC  >  REL606/390606‑390676
                                   |                                   
ggTGGGGTAGCAACGACCATCGGCGGGCTGATTTTc                                     <  1:2721394/36‑1 (MQ=255)
 gTGGGGTAGCAACGACCATCGGCTGGCTGATTTTcc                                    <  1:5759743/36‑1 (MQ=255)
 gTGGGGTAGCAACGACCATCGGCGGGCTGATTTTcc                                    >  1:7101167/1‑36 (MQ=255)
 gTGGGGTAGCAACGACCATCGGCGGGCTGATTTTcc                                    <  1:5595889/36‑1 (MQ=255)
  tGGGGTAGCAACGACCATCGGCGGGGTGATTTTTCt                                   >  1:4973893/1‑36 (MQ=37)
  tGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCt                                   <  1:1548633/36‑1 (MQ=255)
  tGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCt                                   <  1:6612375/36‑1 (MQ=255)
   ggggTAGCAACGACCATCGGCGGGCTTATTTTTCTg                                  >  1:1715262/1‑36 (MQ=37)
     ggTAGCAACGACCATCGGCGGGCTGATTTTCCTTct                                >  1:299476/1‑36 (MQ=255)
      gTAGCAACGACCATCGGCGGGCTGGTTTTCCTGctc                               >  1:6696948/1‑36 (MQ=255)
      gTAGCAACGACCATCGGCGGGCTGATTTTCCTGctc                               <  1:4634170/36‑1 (MQ=255)
        aGCAACGACCATCGGCGGGCTGATTTTCCTGCTCtt                             <  1:6125075/36‑1 (MQ=255)
        aGCAACGACCATCGGCGGGCTGATTTTCCTGCCCtt                             >  1:1329857/1‑36 (MQ=255)
         gCAACGACCATCGGCGGGGCGATTTTTCCGCTCttt                            >  1:5684846/1‑36 (MQ=16)
            aCGTCCATCGGCGGGCTGATTTTCCTGCTCTttatt                         <  1:837283/36‑1 (MQ=255)
             cGACCATCGGCGGTCCGTTTTTTCTCCTCCttttta                        >  1:4340126/1‑36 (MQ=2)
             cGACCATCGGCGGGCTGATTTTCCTTCTCTttatta                        >  1:4492526/1‑36 (MQ=255)
             cGACCATCGGCGGGCTGATTTTCCTGCTCTttatta                        <  1:6645816/36‑1 (MQ=255)
             cGACCATCGGCGGGCTGATTTTCCTGCTCTttagat                        >  1:6121618/1‑33 (MQ=255)
               aCCATCGGCGGGCCCAATTTCCTGCTCTTTATTCTc                      >  1:7221718/1‑36 (MQ=16)
                 cATCGGCGGGCTGATTTTTCTTCTCTTTTTTATctg                    >  1:673458/1‑34 (MQ=21)
                 cATCGGCGGGCTGATTTTTCTTCTCTTTATTATctt                    >  1:1820394/1‑34 (MQ=25)
                  aTCGGCGGGCTGATTTTCCTGCTCTTTATTATCggg                   <  1:1824809/36‑1 (MQ=255)
                    cggcggGCTGATTTTCCTGCTCTTTATTATCGGGtt                 <  1:493227/36‑1 (MQ=255)
                     ggcggGCTGATTTTCCTGCTCTTTATTATCGGGTTg                >  1:4874173/1‑36 (MQ=255)
                     ggcggGCTGATTTTCCTGCTCTTTATTATCGGGTTg                >  1:3497480/1‑36 (MQ=255)
                      gcggGCTGATTTTCCTGCTCTTTATTATCGGGTTGa               >  1:7086485/1‑36 (MQ=255)
                      gcggGCTGATTTTCATGCTCTTTATTATCGGGTTGa               >  1:1268237/1‑36 (MQ=255)
                       cggGCTGATTTTCCTGCTCTTTATTATCGGGTTGAt              >  1:6511589/1‑36 (MQ=255)
                       cggGCTGATTTTCCTGCTCTTTATTATCGGGTTGAt              >  1:5474920/1‑36 (MQ=255)
                        gggCTGATTTTCCTGCTCTTTATTATCGGGTTGAtt             >  1:5522503/1‑36 (MQ=255)
                        gggCTGATTTTCCTGCTCTTTATTATCGGGTTGAtt             <  1:6457370/36‑1 (MQ=255)
                        gggCTGATTTTCCTGCTCTTTATTATCGGCTTGAtt             >  1:3144931/1‑36 (MQ=255)
                            tGATTTTCCTGCTCTTTATTATCGGGTTGATTGGtt         >  1:6171143/1‑36 (MQ=255)
                            tGATTTTCCTGCTCTTTATTATCGGGTTGATTGGtt         <  1:2994377/36‑1 (MQ=255)
                                 ttCCTGCTCTTTATTATCGGGTTGATTGGTTATCAc    >  1:5220577/1‑36 (MQ=255)
                                 ttCCTGCTCTTTATTATCGGGTTGATTGGTTATCAc    >  1:5051127/1‑36 (MQ=255)
                                 ttCCTGCTCTTTATTATCGGGTTGATTGGTTATCAc    <  1:1277946/36‑1 (MQ=255)
                                  tCCTGCTCTTTATTATCGGGTTGATTGGTTATCAcc   >  1:6336770/1‑36 (MQ=255)
                                  tCCTGCTCTTTATTATCGGGTTGATTGGTTATCAcc   <  1:1391342/36‑1 (MQ=255)
                                   ccTGCTCTTTATTATCGGGTTGATTTGTTATCAccc  >  1:7359506/1‑36 (MQ=255)
                                   |                                   
GGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCC  >  REL606/390606‑390676

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 16 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG