Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,196,387 C→T E123K (GAA→AAA)  hldE ← fused heptose 7‑phosphate kinase/heptose 1‑phosphate adenyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,196,3870CT91.7% 24.1 / ‑3.2 12E123K (GAA→AAA) hldEfused heptose 7‑phosphate kinase/heptose 1‑phosphate adenyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (7/4);  total (7/5)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

ATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGATCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTCGGCCAGAGATTTACTCAGCGCGCGCGCTGCATCGTCAATGCCCGTCAACCCGACCAGGCGTGCATTAGCACCGAGAGAAGCGATATTCATCGCCACGTTAGCCGCGCCGCCCGGACGTTCTT  >  NC_000913/3196159‑3196614
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                   
aTTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGATCAACACCTTCGAAACCTTTTTCAAAATCCAGACGGATCAgn                                                                                                                                                                                                               <  1:203781/251‑2 (MQ=255)
                                                     gTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGATCAACACCTTCGAAACCTTTTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTAc                                                                                                                                                          >  1:290152/1‑251 (MQ=255)
                                                         nttttGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGATCAACACCTTCGAAACCTTTTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTAccga                                                                                                                                                       >  1:198372/2‑247 (MQ=255)
                                                                                                                                   ccTTTTGCGTAGTCAGAAAGCACCAGCGCGCCACTCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGATCAACACCTTCGAAACCTTTTTCAAAATCCAGACGGATCAGCTgtt                                                                                                                                                                                                           >  1:250155/1‑124 (MQ=255)
                                                                                                                                                  gAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGATCAACACCTTCGAAACCTTTTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTCGGCCAGAGATTTACTCAGCGCGCGCGCTGCATCGTCAATGCCCGTCAACCCGACCAGGCGTGc                                                             >  1:17107/1‑251 (MQ=255)
                                                                                                                                                      gCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGATCAACACCTTCGAAACCTTTTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTCGGCCAGAGATTTACTCAGCGCGCGCGCTGCATCGTCAATGc                                                                                   <  1:266369/225‑1 (MQ=255)
                                                                                                                                                       cACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGATCGACACCTTCGAAACCTTTTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTCGGCCAGAGATTTACTCAGCGCGCGCGCTGCATCGTc                                                                                        >  1:172118/1‑219 (MQ=255)
                                                                                                                                                           aGCGCGCCAAGCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGATCAACACCTTCGAAACCTTTTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCTGAATTGGGTAAGGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTTACGTTGACGTCGGCCAGAGATTGACTCAGCGCGTGCGCTGCATCGTCAATGCCATTCATCCCGACCAGGCGTGCATTAGCAcc                                                    <  1:214490/251‑1 (MQ=255)
                                                                                                                                                                       gAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGATCAACACCTTCGAAACCTTTTTCACCATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTAccc                                                                                                                                                                                         >  1:53237/1‑106 (MQ=255)
                                                                                                                                                                                                cGTGCAGCGGCTGCGGATCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCTTAATTTGGTAATGGTCGGATGCGGCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTCGGCCAGAGAGTTACTCAGCGCGCGCGCTGCATCGTCAATGCCCGTCAACCCGACCAGGCGTGCATTAGCGCCGAGAGAAGCGATATTCATCGCCACGTTAgccgcgc                 <  1:181586/249‑1 (MQ=255)
                                                                                                                                                                                                  ngCANCGGCTGCGGATCAACACCTTCGAAACCTTTTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATttgacgttgacg                                                                                                                              <  2:198372/137‑1 (MQ=255)
                                                                                                                                                                                                             cGGATCAACACCTTAGAAACCTTTTTCAAAATCCAGAAGGCTCAGCTGTTGGTTGCGGAAAAGTACCCGTAATTTGGTAATGGTCAGATGCGTCGGTACAGAAACAAAGTAGCACTTGACGTTGACGTCGGCCAGAGATTTACTCAGCGCGCGCGCTGCATCGTCAATGCCCGTCAACCCGACCAGGCGTGCATTAGCACCGAGAGAAGCGATATTAATCGACACGTTAGCCGCGCCGCCCGGACGttctt  >  1:306641/1‑251 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                   
ATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGATCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTCGGCCAGAGATTTACTCAGCGCGCGCGCTGCATCGTCAATGCCCGTCAACCCGACCAGGCGTGCATTAGCACCGAGAGAAGCGATATTCATCGCCACGTTAGCCGCGCCGCCCGGACGTTCTT  >  NC_000913/3196159‑3196614

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