Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,285,696 C→T F73F (TTC→TTT agaD → galactosamine‑specific PTS enzyme IID component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,285,6960CT85.7% 26.5 / ‑1.6 14F73F (TTC→TTTagaDgalactosamine‑specific PTS enzyme IID component
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (5/7);  total (6/8)
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.09e-01

CGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCT‑‑ATTTTTCGCCGTTTACCTGCTTA  >  NC_000913/3285527‑3285907
                                                                                                                                                                         |                                                                                                                                                                                                                       
cgtcgCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGTCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCTTATTCCTGATGTTTTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTgg                                                                                                                                         <  1:506438/250‑1 (MQ=255)
           cAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCGGGTCGGATTTCTGATGGGGTTATTAATGTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCgg                                                                                                                             <  1:423458/251‑1 (MQ=255)
               cGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACATGGGCGATGTTGACGATCCTGAAAAAGCTTTATAAGGAAGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGCATTTATTAATACCCACCCGACTCTGGTCGGATTTCTGATGGGGGTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACCCCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGAtt                                                                                                                         >  1:337285/1‑251 (MQ=255)
                                                      ccTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAAGTTAGTAATACCCACCCGAATCTGGTCGGATTTCTGATGGGGTTATTAATTTc                                                                                                                                                                                                     <  1:228904/136‑1 (MQ=255)
                                                                 ttGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTTCTGATTGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGc                                                                                                                                                  <  1:434680/176‑1 (MQ=255)
                                                                          cTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTTCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGcc                                                                                                                                       >  1:216273/1‑178 (MQ=255)
                                                                                       ataAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTAGTAATACCCACCCGAAGCTGGTCGGATTTCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCTCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTGCTGGTTTACTTTGTTGCCGAGTATGGCTGGAATTTGCTCATCATTCGCCAGCCAggg                                                <  1:167447/250‑1 (MQ=255)
                                                                                               gacAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTTCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGATTTGCTCATCATTTTCCAGCCCggggaaaacgct                                       >  1:220978/1‑242 (MQ=255)
                                                                                                                     cgATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTAc                                                                                                  >  1:443522/1‑172 (MQ=255)
                                                                                                                          aaaGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTTCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGccagccag                                                  <  1:294021/213‑1 (MQ=255)
                                                                                                                                  ccTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTTCTGATGGGGATATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGAATTTTTCATTCCTTTTCCGCCAGGGGACCCTGCCGGGGGCGATTCTTTATTTTTCCCGTTTCCTGCttt  >  1:198124/1‑250 (MQ=255)
                                                                                                                                     cGAATTTCTTAATACCCCCCCGAATCTGGTCGGATTTCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGAACTGTTTGGccc                                                                                                                                      >  1:153844/1‑120 (MQ=255)
                                                                                                                                     cGAATTTATTAATACCCACCCGAATCTGGTCGGATTTCTGATGGGGTTATTAAGTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTTTTTGGCCCAATCGCCGGTATTGGCGATGCGATTTTCTGGt                                                                                                      <  1:426793/152‑1 (MQ=255)
                                                                                                                                                                 gTCGGATTTCTGATGGGTTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGgcgatgc                                                                                                                 <  1:188540/113‑1 (MQ=255)
                                                                                                                                                                         |                                                                                                                                                                                                                       
CGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCT‑‑ATTTTTCGCCGTTTACCTGCTTA  >  NC_000913/3285527‑3285907

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