Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,104,492 G→A V6V (GTC→GTT)  wbbK ← putative glycosyltransferase WbbK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,104,4920GA87.5% 11.7 / ‑3.4 8V6V (GTC→GTT) wbbKputative glycosyltransferase WbbK
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (5/2);  total (5/3)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

GAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGACTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGAACCTCTAACAACAGAATTGGCGCCGACTACGACTCCATTACCAATAATTGTTCCAGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCC  >  NC_000913/2104268‑2104701
                                                                                                                                                                                                                                |                                                                                                                                                                                                                 
gAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGACACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACAACTATGCTTTTTCCCATAATTTGGtc                                                                                                                                                                                         >  1:123274/1‑251 (MQ=255)
                                           gCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACAACTATGCTTTTTCCCATAATTTGGTCTCATGATTg                                                                                                                                                                                <  1:6733/217‑1 (MQ=255)
                                           gCCACGAACCTTTGACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACAACTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCa                                                                                                                                               >  1:46093/1‑250 (MQ=255)
                                                                                       ggATAACTTTCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGACTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGAACCTCTAACAACAGAATTGgcg                                                                                                  <  1:53167/251‑1 (MQ=255)
                                                                                                 nttttAACACTTTAGCAGAATGGACTAATGCGACAAAACTGACATTTTCTTTATTATTACTTGCTGCCAAAAATTTTTTCCAACTGGTACATGGACCGCCAGTGGGAAAATTGACAGCAGAAACAACAACTATGCTTTTTCCCATAATTTGGtc                                                                                                                                                                                         >  1:15889/2‑154 (MQ=255)
                                                                                                                       naCTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACAACTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGAACCTCTAACAACAGAATTGGCGCCGACTACGACTCCATTACCAATAATTGTTcc                                                                  >  1:41015/2‑251 (MQ=255)
                                                                                                                                                                       aaaTTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACAACTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGAACCTCTAACAACAGAATTGGCGCCGACTACGACTCCATTACCAATAATTGTTCCAGGCAAAACCGTCACATTCTCACCCAACCAAACTCTCTGGCCAATTac                  <  1:46323/251‑1 (MQ=255)
                                                                                                                                                                                       gTAAATGGACCACCAGTGGTAAAATTGACCGCAGAAACAACAACTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGAACCTCTAACAACAGAATTGGCGCCGACTACGACTCCATTACCAATAATTGTTCCAGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTcc  >  1:10503/1‑251 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                 
GAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGACTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGAACCTCTAACAACAGAATTGGCGCCGACTACGACTCCATTACCAATAATTGTTCCAGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCC  >  NC_000913/2104268‑2104701

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑