Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,748,215 G→A 31.7% S256S (AGC→AGT tnaB ← tryptophan transporter of low affinity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,748,2150GA31.7% 52.8 / 33.1 41S256S (AGC→AGTtnaBtryptophan transporter of low affinity
Reads supporting (aligned to +/- strand):  ref base G (10/18);  new base A (8/5);  total (18/23)
Fisher's exact test for biased strand distribution p-value = 1.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGTTTGGTGCCGAGGAACGATTTCACCAGCGAATCAACGTTGCCGCCTGAGGAGATAATCGCCTTAAAGCTTTCTCGCGGAATATTCCCCATGGTGCAATAGAGCCAGAAGAGATAAATCACCAGCGCC  >  W3110S.gb/3748146‑3748275
                                                                     |                                                            
cTGTTTGGTGCCGAGGAACGATTTCACCAGCGAATCAACGTTGCCGCCTGAGGAGATAATCGCCTTAAAGc                                                             >  1:2352442/1‑71 (MQ=255)
   tttGGTGCCGAGGAACGATTTCACCAGCGAATCAACGTTGCCGCCTGAGGAGATAATCGCCTTAAAGCtt                                                           <  1:101062/70‑1 (MQ=255)
   tttGGTGCCGAGGAACGATTTCACCAGCGAATCAACGTTGCCGCCTGAGGAGATAATCGCCTTAAAACt                                                            <  1:25252/69‑1 (MQ=255)
     tGGTGCCGAGGAACGATTTCACCAGCGAATCAACGTTGCCGCCTGAGGAGATAATCGCCTTAAAACTTtct                                                        >  1:3105960/1‑71 (MQ=255)
     tGGTGCCGAGGAACGATTTCACCAGCGAATCAACGTTGCCGCCTGAGGAGATAATCGCCTTAAAACTTtct                                                        >  1:3656872/1‑71 (MQ=255)
          ccGAGGAACGATTTCACCAGCGAATCAACGTTGCCGCCTGAGGAGATAATCGCCTTAAAGCTTTCTCGCgg                                                   <  1:3697316/71‑1 (MQ=255)
          ccGAGGAACGATTTCACCAGCGAATCAACGTTGCCGCCTGAGGAGATAATCGCCTTAAAACTTTCTCGCgg                                                   <  1:1018569/71‑1 (MQ=255)
          ccGAGGAACGATTTCACCAGCGAATCAACGTTGCCGCCTGAGGAGATAATCGCCTTAAAACTTTCTCGCgg                                                   <  1:1193202/71‑1 (MQ=255)
            gAGGAACGATTTCACCAGCGAATCAACGTTGCCGCCTGAGGAGATAATCGCCTTAAAGCTTTCTCGCGGaa                                                 <  1:953192/71‑1 (MQ=255)
             aGGAACGATTTCACCAGCGAATCAACGTTGCCGCCTGAGGAGATAATCGCCTTAAAGCTTTCTCGCGGaa                                                 <  1:3846054/70‑1 (MQ=255)
              ggAACGATTTCACCAGCGAATCAACGTTGCCGCCTGAGGAGATAATCGCCTTAAAGCTTTCTCGCGGAata                                               <  1:3069670/71‑1 (MQ=255)
               gAACGATTTCACCAGCGAATCAACGTTGCCGCCTGAGGAGATAATCGCCTTAAAGc                                                             <  1:33845/56‑1 (MQ=255)
               gAACGATTTCACCAGCGAATCAACGTTGCCGCCTGAGGAGATAATCGCCTTAAAGCTTTCTCGCGGAata                                               >  1:2934759/1‑70 (MQ=255)
                aaCGATTTCACCAGCGAATCAACGTTGCCGCCTGAGGAGATAATCGCCTTAAAACTTTCTCGCGGAATAtt                                             <  1:3588269/71‑1 (MQ=255)
                    aTTTCACCAGCGAATCAACGTTGCCGCCTGAGGAGATAATCGCCTTAAAGCTTTCTCGCGGAATATTccc                                          >  1:2556103/1‑70 (MQ=255)
                    aTTTCACCAGCGAATCAACGTTGCCGCCT‑‑GGAGATAATCGCCTTAAAGCtt                                                           <  1:91114/51‑1 (MQ=37)
                     tttCACCAGCGAATCAACGTTGCCGCCTGAGGAGATAATCGCCTTAAAGCTTTCTCGCGGaa                                                 <  1:2709265/62‑1 (MQ=255)
                           cAGCGAATCAACGTTGCCGCCTGAGGAGATAATCGCCTTAAAGCTTTCTCGCGGAATATTCCCCATGGTGc                                  <  1:249777/71‑1 (MQ=255)
                            aGCGAATCAACGTTGCCGCCTGAGGAGATAATCGCCTTAAAGCTTTCTCGCGGaa                                                 <  1:966057/55‑1 (MQ=255)
                            aGCGAATCAACGTTGCCGCCTGAGGAGATAATCGCCTTAAAGCTTTCTCGCGGAATATTCCCCATGGTGc                                  >  1:3962973/1‑70 (MQ=255)
                              cGAATCAACGTTGCCGCCTGAGGAGATAATCGCCTTAAAACTTTCTCGCGGAata                                               >  1:1228353/1‑55 (MQ=255)
                              cGAATCAACGTTGCCGCCTGAGGAGATAATCGCCTTAAAACTTTCTCGCGGAata                                               >  1:628938/1‑55 (MQ=255)
                              cGAATCAACGTTGCCGCCTGAGGAGATAATCGCCTTAAAACTTTCTCGCGGAata                                               >  1:2374875/1‑55 (MQ=255)
                               gAATCAACGTTGCCGCCTGAGGAGATAATCGCCTTAAAGCTTTCTCGCGGAATATTCCCCATGGTGCaa                                >  1:1902328/1‑69 (MQ=255)
                               gAATCAACGTTGCCGCCTGAGGAGATAATCGCCTTAAAGCTTTCTCGCGGAATATTCCCCATGGTGCaa                                >  1:1359720/1‑69 (MQ=255)
                                  tCAACGTTGCCGCCTGAGGAGATAATCGCCTTAAAGCTTTCTCGCGGaa                                                 <  1:3302424/49‑1 (MQ=255)
                                  tCAACGTTGCCGCCTGAGGAGATAATCGCCTTAAAACTTTCTCGCGTAATATTCCCCATGGTGCAATagag                           >  1:2872248/1‑71 (MQ=255)
                                      cGTTGCCGCCTGAGGAGATAATCGCCTTAAAGCTTTCTCGCGGAATATTCCCCATGGTGCAATAGAGCCag                       <  1:2161752/71‑1 (MQ=255)
                                      cGTTGCCGCCTGAGGAGATAATCGCCTTAAAGCTTTCTCGCGGAATATTCCCCATGGTGCAATAGAGCCag                       <  1:1715273/71‑1 (MQ=255)
                                                    gagaTAATCGCCTTAAAACTTTCTCGCGGAATATTCCCCATGGTGCAATAGAGCCAGAAGAGATAAATCAc         <  1:1869114/71‑1 (MQ=255)
                                                     agaTAATCGCCTTAATGCTTTCTCGCGGAATATTCCCCATGGTGCAATAGAGCCAGAAGAGATAAATCAc         <  1:1823345/70‑1 (MQ=255)
                                                       aTAATCGCCTTAAAGCTTTCTCGCGGAATATTCCCCATGGTGCAAt                               >  1:307251/1‑46 (MQ=255)
                                                       aTAATCGCCTTAAAGCTTTCTCGCGGAATATTCCCCATGGTGCAATAGAGCCagaaga                   <  1:3041461/58‑1 (MQ=255)
                                                         aaTCGCCTTAAAGCTTTCTCGCGGAATATTCCCCATGGTGCAATAGAGCCAGAAgaga                 >  1:479435/1‑58 (MQ=255)
                                                          aTCGCCTTAAAGCTTTCTCGCGGAATATTCCCCATGGTGCAATAGAGCCAGAAGAGATa               >  1:1181762/1‑59 (MQ=255)
                                                           tggCCTTAAAGCTTTCTCGCGGAATATTCCCCATGGTGCAATAGAGCCAGAAGAGATAAATCACCAGCGcc  <  1:3323177/69‑1 (MQ=255)
                                                             gCCTTAAAACTTTCTCGCGGAATATTCCCCATGGTGCAATa                              >  1:642452/1‑41 (MQ=255)
                                                             gCCTTAAAACTTTCTCGCGGAATATTCCCCATGGTGCAATa                              >  1:2550455/1‑41 (MQ=255)
                                                              ccTTAAAGCTTTCTCGCGGAATATTCCCCATGGTGCa                                 <  1:1301168/37‑1 (MQ=255)
                                                              ccTTAAAGCTTTCTCGCGGAATATTCCCCATGGTGCa                                 <  1:2614374/37‑1 (MQ=255)
                                                                  aaaaCTTTCTCGCGGAATATTCCCCATGGTGCAAt                               >  1:3237445/1‑35 (MQ=255)
                                                                  aaaaCTTTCTCGCGGAATATTCCCCATGGTGCAAt                               >  1:1935817/1‑35 (MQ=255)
                                                                   aaGATTTCTCGCGGAATATTCCCCATGGTGCAATa                              >  1:2252837/1‑35 (MQ=255)
                                                                     |                                                            
CTGTTTGGTGCCGAGGAACGATTTCACCAGCGAATCAACGTTGCCGCCTGAGGAGATAATCGCCTTAAAGCTTTCTCGCGGAATATTCCCCATGGTGCAATAGAGCCAGAAGAGATAAATCACCAGCGCC  >  W3110S.gb/3748146‑3748275

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