Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 344,381 G→A intergenic (+166/‑247) yahL → / → yahM hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb344,3810GA98.2% 172.1 / ‑4.1 55intergenic (+166/‑247)yahL/yahMhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (21/33);  total (22/33)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CCATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACAT  >  W3110S.gb/344314‑344448
                                                                   |                                                                   
ccATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGttt                                                                  >  1:1340764/1‑71 (MQ=255)
ccATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGAttt                                                                  >  1:1015195/1‑71 (MQ=255)
ccATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGAttt                                                                  >  1:372785/1‑71 (MQ=255)
ccATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGAttt                                                                  >  1:3556545/1‑71 (MQ=255)
ccATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGAttt                                                                  >  1:2863900/1‑71 (MQ=255)
ccATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGAttt                                                                  >  1:2726299/1‑71 (MQ=255)
ccATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGAttt                                                                  >  1:4128689/1‑71 (MQ=255)
ccATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGAttt                                                                  >  1:245307/1‑71 (MQ=255)
ccATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGAttt                                                                  >  1:2069882/1‑71 (MQ=255)
ccATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGAttt                                                                  >  1:1680195/1‑71 (MQ=255)
ccATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGAttt                                                                  >  1:1508896/1‑71 (MQ=255)
     tctTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAg                                                             >  1:3396046/1‑71 (MQ=255)
     tctTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAg                                                             >  1:1460799/1‑71 (MQ=255)
     tctTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAg                                                             >  1:411515/1‑71 (MQ=255)
     tctTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAg                                                             >  1:1578777/1‑71 (MQ=255)
               ccGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGgtg                                                   <  1:4056651/71‑1 (MQ=255)
                             tgATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGAtc                                     >  1:2040253/1‑71 (MQ=255)
                             tgATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGAtc                                     >  1:1391847/1‑71 (MQ=255)
                                tCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCa                                  >  1:3765034/1‑71 (MQ=255)
                                tCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCa                                  >  1:1826288/1‑71 (MQ=255)
                                  ggCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATg                                <  1:1599113/71‑1 (MQ=255)
                                  ggCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATg                                <  1:2548806/71‑1 (MQ=255)
                                  ggCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATg                                <  1:2723588/71‑1 (MQ=255)
                                  ggCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATg                                <  1:1776503/71‑1 (MQ=255)
                                  ggCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATg                                <  1:3904021/71‑1 (MQ=255)
                                   gCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATg                                <  1:2041108/70‑1 (MQ=255)
                                      acacTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCgg                                                     <  1:3187533/46‑1 (MQ=255)
                                      acacTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCgg                                                     <  1:555870/46‑1 (MQ=255)
                                               ttAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGt                         <  1:3514556/65‑1 (MQ=255)
                                               ttAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGt                         <  1:2426281/65‑1 (MQ=255)
                                               ttAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGt                         <  1:2326061/65‑1 (MQ=255)
                                               ttAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGt                         <  1:2091506/65‑1 (MQ=255)
                                               ttAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGt                         <  1:820485/65‑1 (MQ=255)
                                                tAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATg                                <  1:1362110/57‑1 (MQ=255)
                                                 aaaaaaGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTtgtg                   <  1:931486/69‑1 (MQ=255)
                                                 aaaaaaGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTtgtg                   <  1:3735562/69‑1 (MQ=255)
                                                 aaaaaaGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTtgtg                   <  1:37816/69‑1 (MQ=255)
                                                 aaaaaaGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTtgtg                   <  1:700664/69‑1 (MQ=255)
                                                 aaaaaaGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTtgtg                   <  1:3433280/69‑1 (MQ=255)
                                                 aaaaaaGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTtgtg                   <  1:1073992/69‑1 (MQ=255)
                                                 aaaaaaGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTtgtg                   <  1:189942/69‑1 (MQ=255)
                                                 aaaaaaGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTtgtg                   <  1:2652837/69‑1 (MQ=255)
                                                 aaaaaaGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTtgtg                   <  1:2837184/69‑1 (MQ=255)
                                                 aaaaaaGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTtgtg                   <  1:3121184/69‑1 (MQ=255)
                                                 aaaaaaGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTtgtg                   <  1:3666771/69‑1 (MQ=255)
                                                 aaaaaaGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGATTtgtg                   <  1:6247/69‑1 (MQ=255)
                                                         aCCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTgttgg         <  1:3819198/71‑1 (MQ=255)
                                                         aCCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTgttgg         <  1:2703034/71‑1 (MQ=255)
                                                         aCCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTgttgg         <  1:2971083/71‑1 (MQ=255)
                                                         aCCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTgttgg         <  1:2525706/71‑1 (MQ=255)
                                                         aCCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTgttgg         <  1:1636121/71‑1 (MQ=255)
                                                          ccAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTtgtg                   <  1:717434/60‑1 (MQ=255)
                                                                gtGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACAt  >  1:4116040/1‑71 (MQ=255)
                                                                gtGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACAt  >  1:4120717/1‑71 (MQ=255)
                                                                gtGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACAt  >  1:1758780/1‑71 (MQ=255)
                                                                   |                                                                   
CCATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACAT  >  W3110S.gb/344314‑344448

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: