Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,774,603 G→A G272D (GGC→GAC)  ydiB → quinate/shikimate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,774,6030GA85.7% 10.9 / ‑4.8 7G272D (GGC→GAC) ydiBquinate/shikimate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (4/2);  minor base T (1/0);  total (5/2)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

TTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGGGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAATGAAAACCGTAACTGTAAAAGATCTCGTCATTGGTACGGGCGCACCTAAAATCATCGTCTCGCTGATGGCGAAAGATATCGCCAGCGTGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATATTCTGGAATGG  >  NC_000913/1774386‑1774826
                                                                                                                                                                                                                         |                                                                                                                                                                                                                               
ttAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGGGGCTGAACAGTTCACATTATGGACTGACAAAGATTTCCCTCTGGAATATGTTAAACAgg                                                                                                                                                                                                 <  1:62010/250‑1 (MQ=255)
                                                                              ggACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGGGGCTGAACAGTTCACATTATGGACTGACAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAATGAAAACCGTAACTGTAAAAGATCTCGt                                                                                                                  <  1:58329/251‑1 (MQ=255)
                                                                                           cTGAATGCGTGTATAACCCGCATATGAAGAAGTTATTGCAGCAGGCGCAACCAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGGGGCTGAACAGTTCACATTATGGACTGACAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAATGAAAACCGTAACTGTAAAAGATCTCGTCATTGGTACGGgc                                                                                                     >  1:72004/1‑251 (MQ=255)
                                                                                                                       gAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGGGGCTGAACAGTTCACATTATGGACTGACAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGAAGAAAGGGAAAAAAATGAAAACCGAAACTGTAAAAGATCTCGTCATTGGTACGGGCGCACCTAAAATCATCGTCTCGCTGATgg                                                                         >  1:4234/1‑251 (MQ=255)
                                                                                                                                                                          cGGCATGTTGTTGTGGCAAGGGGCTGAACAGTTCACATTATGGACTGTCAAAGATTTCCCACTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCATGACAGGCTGACCGCGTGCAGAAAGGGAAAAAAATGAAAACCGAAACTGAAAAAGATCTCGCCATTGGTACGGGCGCAACTAAAATCATCGTCTCGCTGATGGCGAAAGAAATCGCCAGCGTGAAATCCGAAGCTCTCGACTATCGTGAAGCgg                      >  1:116015/1‑251 (MQ=255)
                                                                                                                                                                                          cAAGGGGCTGAACAGTTCACATTATGGACTGACAAAGATTTCCCTCTGGACTATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAATGAAAACCGTAACTGTAAAAGATCTCGTCCTTGGTACGGGCGCACCTAAAATCATCGTCTCGCTGATGGCGAAAGATATCGCCAGCGTGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATATTCTgg       >  1:27682/1‑250 (MQ=255)
                                                                                                                                                                                              gggCTGAACAGTTCACATTATGGACTGACAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAATGAAAACCGTAACTGTAAAAGATCTCGTCATTGGTACGGGCGCACCTAAAATCATCGTCTCGCTGATGGCGAAAGATATCGCCAGCGTGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATCTTCTGGAATgg  >  1:130458/1‑251 (MQ=255)
                                                                                                                                                                                                                         |                                                                                                                                                                                                                               
TTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGGGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAATGAAAACCGTAACTGTAAAAGATCTCGTCATTGGTACGGGCGCACCTAAAATCATCGTCTCGCTGATGGCGAAAGATATCGCCAGCGTGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATATTCTGGAATGG  >  NC_000913/1774386‑1774826

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