Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,111,836 G→A S325N (AGC→AAC)  opgH → osmoregulated periplasmic glucans biosynthesis protein H

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,111,8360GA83.3% 10.7 / NA 6S325N (AGC→AAC) opgHosmoregulated periplasmic glucans biosynthesis protein H
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (1/4);  minor base C (1/0);  total (2/4)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACCCGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGCCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTAC  >  NC_000913/1111648‑1112064
                                                                                                                                                                                            |                                                                                                                                                                                                                                    
gCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACCCGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAACGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTg                                                                                                                                                                         <  1:15231/250‑1 (MQ=255)
                                               ttGATGTCTACATTCTTAGTGACAGTTATAACCCGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAACGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTgcn                                                                                                                          <  1:24715/250‑2 (MQ=255)
                                                               tAGTGACAGTTATAACCCGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAACGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGGGCTGGATGCGGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAAc                                                                                                          <  1:81139/250‑1 (MQ=255)
                                                                                                 nnnnanaannccttngtggaGCTTATCGCTGAAGTCGGGNGCGAAGGTCAGATTTTCTAGCGCCGCCTCCGNCNCCGCGTGAAGCGTAAAAACGGTAATATCGATGACTTCTGCCGTCGCTGNGGCAGCCAGTNCAGNNACATNGNGNTGCAGGATGCTGACTCGGTAATGACCGGTNATTGTTTGTGCGGGCTGGTGCGCNNGANNNAANNCAAcc                                                                                                         <  2:156768/201‑1 (MQ=255)
                                                                                                                                  gTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAACCGGTAATCTCGATGACTTATGCCGTCACTGGGGCAGCCAGTACCGCTACCTGGTGGTGATGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGCCGGGATCATTCAGTCGTCGCCGAAAg                                                                        >  1:200020/1‑217 (MQ=255)
                                                                                                                                                                      cgTCGCCGCGTGAAGCGTAAAAACGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGCCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTAc  >  1:99820/1‑251 (MQ=255)
                                                                                                                                                                                            |                                                                                                                                                                                                                                    
GCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACCCGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGCCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTAC  >  NC_000913/1111648‑1112064

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