Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 344,381 G→A 33.3% 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,3810GA33.3% 36.7 / 25.5 34intergenic (+166/‑247)yahL/yahMhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (12/10);  new base A (6/5);  total (18/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCA  >  W3110S.gb/344314‑344451
                                                                   |                                                                      
ccATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGAttt                                                                     >  1:706240/1‑71 (MQ=255)
     tctTTTTTCGCCGATTTTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAg                                                                >  1:2142749/1‑71 (MQ=255)
     tctTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAg                                                                >  1:818267/1‑71 (MQ=255)
     tctTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAg                                                                >  1:3298562/1‑71 (MQ=255)
     tctTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAg                                                                >  1:2970511/1‑71 (MQ=255)
     tctTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAg                                                                >  1:434916/1‑71 (MQ=255)
     tctTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAg                                                                >  1:2498473/1‑71 (MQ=255)
     tctTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAg                                                                >  1:2413455/1‑71 (MQ=255)
     tctTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAg                                                                >  1:2049195/1‑71 (MQ=255)
     tctTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAg                                                                >  1:1618211/1‑71 (MQ=255)
     tctTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAg                                                                >  1:1612585/1‑71 (MQ=255)
     tctTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAg                                                                >  1:1441410/1‑71 (MQ=255)
     tctTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAg                                                                >  1:660690/1‑71 (MQ=255)
                    atTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCAt                                                            >  1:2219863/1‑60 (MQ=255)
                    atTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGtata                                                  >  1:411427/1‑70 (MQ=255)
                    atTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGtata                                                  >  1:3221213/1‑70 (MQ=255)
                    atTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGtata                                                  >  1:394547/1‑70 (MQ=255)
                    atTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGtata                                                  >  1:1543625/1‑70 (MQ=255)
                       ttgttgTGATCGGCGCCACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATaaaa                                               <  1:3075490/70‑1 (MQ=255)
                       ttgttgTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATaaaa                                               <  1:2545898/70‑1 (MQ=255)
                       ttgttgTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATaaaa                                               <  1:2734515/70‑1 (MQ=255)
                       ttgttgTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATaaaa                                               <  1:2462708/70‑1 (MQ=255)
                                    cGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCa                                 <  1:2887864/71‑1 (MQ=255)
                                         cTTCGCTTAAAAAAGCACCAGTAGTGCTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGt                            <  1:3132975/71‑1 (MQ=255)
                                         cTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGt                            <  1:48531/71‑1 (MQ=255)
                                         cTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGt                            <  1:507315/71‑1 (MQ=255)
                                         cTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGt                            <  1:2862093/71‑1 (MQ=255)
                                                          gcAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGt                            <  1:2699345/53‑1 (MQ=255)
                                                             gtagtGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGa        <  1:510557/71‑1 (MQ=255)
                                                             gtagtGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGa        <  1:707296/71‑1 (MQ=255)
                                                             gtagtGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGa        <  1:851985/71‑1 (MQ=255)
                                                             gtagtGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGa        <  1:1871051/71‑1 (MQ=255)
                                                             gtagtGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGa        <  1:1843606/71‑1 (MQ=255)
                                                             gtagtGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGa        <  1:1548117/71‑1 (MQ=255)
                                                                   gTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCa  <  1:2543238/71‑1 (MQ=255)
                                                                   gTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCa  <  1:1392356/71‑1 (MQ=255)
                                                                   gTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCa  <  1:1174099/71‑1 (MQ=255)
                                                                   |                                                                      
CCATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCA  >  W3110S.gb/344314‑344451

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