Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,794,102 G→A 75.0% L398L (TTG→TTA gabP → gamma‑aminobutyrate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,794,1020GA75.0% 117.0 / 44.7 68L398L (TTG→TTAgabPgamma‑aminobutyrate transporter
Reads supporting (aligned to +/- strand):  ref base G (6/11);  new base A (25/26);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.03e-01

TCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTT  >  W3110S.gb/2794038‑2794161
                                                                |                                                           
tCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGa                                                          <  1:1041419/68‑1 (MQ=255)
 cATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGaa                                                         <  1:308779/68‑1 (MQ=255)
  aTCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCg                                                           <  1:193842/65‑1 (MQ=255)
     gCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTgg                                                    >  1:3631986/1‑69 (MQ=255)
     gCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTgg                                                    >  1:3336315/1‑69 (MQ=255)
     gCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTgg                                                    >  1:1484086/1‑69 (MQ=255)
       cttttCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCt                                                  <  1:2875755/67‑1 (MQ=255)
       cGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCt                                                  <  1:369766/69‑1 (MQ=255)
       cGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCt                                                  <  1:2016547/69‑1 (MQ=255)
        gTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCtt                                                 <  1:3512147/69‑1 (MQ=255)
        gTTTCACAGTTGCGGATGCGCAAAATTCTGCGACCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCtt                                                 <  1:1915458/69‑1 (MQ=255)
          ttCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTa                                               <  1:1088019/69‑1 (MQ=255)
          ttCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTa                                               <  1:2766079/69‑1 (MQ=255)
          ttCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTa                                               <  1:2552407/69‑1 (MQ=255)
           tCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTAt                                              >  1:936607/1‑69 (MQ=255)
           tCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTAt                                              <  1:554232/69‑1 (MQ=255)
           tCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTAt                                              >  1:1109784/1‑69 (MQ=255)
             acaGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTATc                                             <  1:3239328/68‑1 (MQ=255)
                   gCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGc                                                   >  1:3371104/1‑56 (MQ=255)
                   gCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGc                                                   >  1:2427180/1‑56 (MQ=255)
                   gCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCtt                                                 <  1:3474199/58‑1 (MQ=255)
                   gCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCtt                                                 <  1:76967/58‑1 (MQ=255)
                   gCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCtt                                                 <  1:2688126/58‑1 (MQ=255)
                   gCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCt                                      >  1:2198905/1‑69 (MQ=255)
                   gCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCt                                      >  1:1529404/1‑69 (MQ=255)
                     ggATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTc                                     <  1:1919147/68‑1 (MQ=255)
                     ggATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTc                                     <  1:1579422/68‑1 (MQ=255)
                     ggATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTc                                     <  1:1642780/68‑1 (MQ=255)
                     ggATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTc                                     <  1:827689/68‑1 (MQ=255)
                      gATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCtt                                                 <  1:1194893/55‑1 (MQ=255)
                        tGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCtt                                                 <  1:3602813/53‑1 (MQ=255)
                        tGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCtt                                                 <  1:212039/53‑1 (MQ=255)
                        tGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCtt                                                 <  1:2643026/53‑1 (MQ=255)
                         gcgcAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTg                                >  1:2892376/1‑69 (MQ=255)
                         gcgcAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTg                                >  1:1884192/1‑69 (MQ=255)
                         gcgcAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGc                                                   <  1:380423/50‑1 (MQ=255)
                         gcgcAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGc                                                   <  1:2704685/50‑1 (MQ=255)
                         gcgcAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGc                                                   <  1:3435605/50‑1 (MQ=255)
                         gcgcAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGc                                                   <  1:1489496/50‑1 (MQ=255)
                         gcgcAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGc                                                   <  1:2875764/50‑1 (MQ=255)
                         gcgcAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTg                                >  1:1755225/1‑69 (MQ=255)
                          cgcAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTgg                                        >  1:2885837/1‑60 (MQ=255)
                          cgcAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCt                                      >  1:3491636/1‑62 (MQ=255)
                           gcAAAATTCTGCGAGCAGAAGGAAGCGAATTTCGCTTACGAATGTGGGtt                                                 <  1:3175823/50‑1 (MQ=255)
                                aTTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTc                         >  1:2625723/1‑69 (MQ=255)
                                aTTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTc                         >  1:1500585/1‑69 (MQ=255)
                                        aGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCtt                 <  1:668672/69‑1 (MQ=255)
                                        aGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCtt                 <  1:1170046/69‑1 (MQ=255)
                                          cAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCAc                                   <  1:2948160/49‑1 (MQ=255)
                                           aGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTttat              <  1:195315/69‑1 (MQ=255)
                                           aGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTttat              <  1:873559/69‑1 (MQ=255)
                                           aGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTttat              <  1:1983171/69‑1 (MQ=255)
                                           aGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATAGGCTttat              <  1:1254576/69‑1 (MQ=255)
                                               ggaagCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTAcc          >  1:1983045/1‑69 (MQ=255)
                                               ggaagCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTAcc          >  1:2134716/1‑69 (MQ=255)
                                               ggaagCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTAcc          >  1:3189648/1‑69 (MQ=255)
                                               ggaagCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTAcc          >  1:664326/1‑69 (MQ=255)
                                               ggaagCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTAcc          >  1:1156509/1‑69 (MQ=255)
                                               ggaagCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTAcc          >  1:1557874/1‑69 (MQ=255)
                                                  agCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGc                              >  1:1595936/1‑46 (MQ=255)
                                                     gAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTc                                     >  1:512691/1‑36 (MQ=255)
                                                     gAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTc                                     >  1:1191485/1‑36 (MQ=255)
                                                     gAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTc                                     >  1:1895680/1‑36 (MQ=255)
                                                     gAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTc                                     >  1:852128/1‑36 (MQ=255)
                                                        aTTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGtt  >  1:2351820/1‑68 (MQ=255)
                                                        aTTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGtt  >  1:1537246/1‑68 (MQ=255)
                                                        aTTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGtt  >  1:2222416/1‑68 (MQ=255)
                                                         ttCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCt                  >  1:3176390/1‑51 (MQ=255)
                                                                |                                                           
TCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTT  >  W3110S.gb/2794038‑2794161

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: