Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 405,632 A→G 19.4% T2A (ACA→GCA)  aroL → shikimate kinase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb405,6320AG19.4% 74.4 / 13.0 36T2A (ACA→GCA) aroLshikimate kinase II
Reads supporting (aligned to +/- strand):  ref base A (16/13);  new base G (5/2);  total (21/15)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCTCATGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCTTGCCGATT  >  W3110S.gb/405577‑405701
                                                       |                                                                     
ttCTCATGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGgcac                                                                    <  1:2524007/59‑1 (MQ=255)
  ctcATGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGGCACaa                                                                  >  1:221117/1‑59 (MQ=255)
  ctcATGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGGCACAACCTCttttt                                                         >  1:3487738/1‑68 (MQ=255)
  ctcATGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCt                                                       <  1:2788388/70‑1 (MQ=255)
     aTGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGGCACAACCTCTTTTTCTg                                                      <  1:2485056/68‑1 (MQ=255)
     aTGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGGCACAACCTCTTTTTCTg                                                      <  1:1511783/68‑1 (MQ=255)
     aTGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTg                                                      <  1:3391263/68‑1 (MQ=255)
     aTGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTg                                                      <  1:3169897/68‑1 (MQ=255)
     aTGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTg                                                      <  1:1664992/68‑1 (MQ=255)
     aTGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTg                                                      <  1:821638/68‑1 (MQ=255)
               cTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGgca                                                                     >  1:3118074/1‑43 (MQ=255)
                  tCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGc                                       <  1:1633181/70‑1 (MQ=255)
                       gcATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTg                                                      <  1:1000915/50‑1 (MQ=255)
                         aTTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTaa                               <  1:1987805/71‑1 (MQ=255)
                         aTTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTaa                               <  1:513987/71‑1 (MQ=255)
                            gCGACCTATTGGGGAAAACCCACGATGACACAACCTCttttt                                                         >  1:3087748/1‑42 (MQ=255)
                            gCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGa                                                     >  1:310985/1‑46 (MQ=255)
                            gCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGAt                                                    >  1:1490487/1‑47 (MQ=255)
                            gCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGAt                                                    >  1:2372418/1‑47 (MQ=255)
                            gCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGAt                                                    >  1:855786/1‑47 (MQ=255)
                            gCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGAt                                                    >  1:2536503/1‑47 (MQ=255)
                            gCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTaaaa                             >  1:1090635/1‑70 (MQ=255)
                            gCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTaaaa                             >  1:244907/1‑70 (MQ=255)
                                   aTTGGGGAAAACCCACGATGGCACAACCTCttttt                                                         >  1:2393117/1‑35 (MQ=255)
                                   aTTGGGGAAAACCCACGATGGCACAACCTCTTTTTCTGAt                                                    >  1:2867570/1‑40 (MQ=255)
                                      ggggAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCgg                   <  1:2811872/70‑1 (MQ=255)
                                      ggggAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCgg                   <  1:2425266/70‑1 (MQ=255)
                                      ggggAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCgg                   >  1:1298047/1‑70 (MQ=255)
                                         gAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATg               >  1:3629263/1‑71 (MQ=255)
                                          aaaaCCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCtg                                     >  1:2548552/1‑48 (MQ=255)
                                          aaaaCCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTaa                               <  1:711388/54‑1 (MQ=255)
                                            aaCCCACGATGGCACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGcc            >  1:939197/1‑71 (MQ=255)
                                            aaCCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTaaaa                             >  1:1994743/1‑54 (MQ=255)
                                             aCCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGccc           <  1:974577/71‑1 (MQ=255)
                                              cccACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATgg              >  1:2832447/1‑67 (MQ=255)
                                              cccACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGc             >  1:1913676/1‑68 (MQ=255)
                                               ccACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCtt         >  1:3068813/1‑71 (MQ=255)
                                                      gACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCTTGCCGAtt  >  1:415551/1‑71 (MQ=255)
                                                       acacAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTaa                               <  1:1482170/41‑1 (MQ=255)
                                                       acacAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTaa                               <  1:632895/41‑1 (MQ=255)
                                                       acacAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTaa                               <  1:404688/41‑1 (MQ=255)
                                                       acacAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTaa                               <  1:1224598/41‑1 (MQ=255)
                                                       acacAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGaa                 >  1:2927815/1‑55 (MQ=255)
                                                       |                                                                     
TTCTCATGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCTTGCCGATT  >  W3110S.gb/405577‑405701

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: