Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 374,464 C→T H199H (CAC→CAT mhpE → 4‑hydroxy‑2‑oxovalerate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913374,4640CT95.0% 47.8 / ‑5.9 20H199H (CAC→CATmhpE4‑hydroxy‑2‑oxovalerate aldolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (8/11);  minor base G (0/1);  total (8/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

CCCGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAGCTGATGGAAGGCTACGGTGCGACCTGTATTTATGTGGTGGATTCTGGCGGTGCGATGAACATGAGCGATATCCGTGACCGTTTCCGCGCCCTGAAAGCAGAGCTGAAACCAGAAACGCAAACTGGCATGCACGCTCACCATAACCTGAGTCTTGGCGTGGCGAACTCTATCGCGGCGGTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAG  >  NC_000913/374244‑374690
                                                                                                                                                                                                                            |                                                                                                                                                                                                                                  
cccGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAGCTGATGGAAGGCTACGGTGCGACCTGTATTTATGTGGTGGATTCTGGCGGTGCGATGAACATGAGCGATATCCGTGACCGTTTCCGCGCCCTGAAAGCAGAGCTGAAACCAGAAACGCAAACTGGCATGCACGCTCATCATAACCTGAGTCTTGGCGTGGCGAActct                                                                                                                                                                                                      >  1:177549/1‑251 (MQ=255)
                    acCGTTGGTGTTCTGATGATGAGCCATATGACCACGCCGGAGAAGCTCGCAAAGCAGGCAAAGCTGATGGAAGGCTACGGTGCGACCGGTATTTATGTGGTGGATTCTGGCGGTGCGATGAACATGAGCGATAGCCGTGACCGTTTCCGCGCCCTGAAAGCAGAGCTGAAACCAGAAACGCAAACTGGCATGCACGCTCATCATAACCTGATTCTTGGCGTGGCGAACTCTATCGCGGCGGTGGAAGAggn                                                                                                                                                                                  <  1:46917/251‑2 (MQ=255)
                                                tGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAGCTGATGCAAGGCTAAGGTACCAAATGCCATCATGAAGGGCATGACGTCGGAGCGATGAAAATGAGCGATATTCGAGACACTTTCCGAGACCTGAAAGCAGAGCTGAAAACAGAACAGCAAACTGCCATGCAACCTAATCTTAACCTGAGTCTTGGAGTGGCGAACTCTATCGCGGAGGAGGAAGTGGGCTGTGACCGAATCGACGCACGCCTcgcg                                                                                                                                                      >  1:376513/1‑251 (MQ=255)
                                                     aCGCCGGAGAATCTCGCCAAGCAGGCAAAGCGGATGGAAGGCTACGGTGCGACCTGTATTTATGTGGTGGATTCTGGCGGTGCGATGAACATGAGCGATATCCGTGACCGTTTCCGCGCCCTGAAAGCAGAGCTGAAACCAGAAACGCAAACTGGCATGCACGCTCATCATAACCTGAGTCTTGGCGTGGCGAACTCTATCGCGGCGGTGGAAGAGGGCtn                                                                                                                                                                               <  1:149093/221‑2 (MQ=255)
                                                                                                aCGGTGCGACCTGTATTTATGTGGTGGATTCTGGCGGTGCGATGAACATGAGCGATATCCGTGACCGTTTCCGCGCCCTGAAAGCAGAGCTGAAACCAGAAACGCAAACTGGCATGCACGCTCATCATAACCTGAGTCTTGGCGTGGCGAACTCTATCGCGGCGGTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGTAGGTAACGCACCGCTGGAAGTGTTTATTgcc                                                                                                      <  1:418871/251‑1 (MQ=255)
                                                                                                aCGGTGCGACCTGTATTTATGTGGTGGATTCTGGCGGTGCGATGAACATGAGCGATATCCGTGACCGTTTCCGCGCCCTGAAAGCAGAGCTGAAACCAGAAACGCAAACTGGCATGCACGCTCATCATAACCTGAGTCTTGGCGTGGCGAACTCTATCGCGGCGGTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTgcc                                                                                                      <  1:418864/251‑1 (MQ=255)
                                                                                                     gCGACCTGTATTTATGTGGTGGATTCGGGCGGTGCGATAAACATGAGCGATATCCGTGACCGTTTCCGCGCCCTAAAAGCAGAGCTTAAACCAGAAACGCAAACTGGCATGCACGCTCAGCGTAACCTGATTCTTGGCGTGGCGATCTCTATCGCGGCGGTGGAAGAGGGCTGCGACCGAATCGACg                                                                                                                                                                 <  1:503908/187‑1 (MQ=255)
                                                                                                                    gtggtgGATTCTGGCGGTGCGATGAACATGAGCGATATCCGTGACCGTTTCCGCGCCCTGAAAGCAGAGCTGAAACCAGAAACGCAAACTGGCATGCACGCTCATCATAACCTGAGTCTTGGCGTGGCGAACTCTATCGCGGCGGTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTg                                                                                  <  1:399713/251‑1 (MQ=255)
                                                                                                                      ggtggATTCTGGAGGTGCGATGAACATGAGCGATATCCGTGATCGTTTCCGCGCCCTGAAAGCAGAGCTGAAACCAGAAACGCAAACTGGCATGCACGCTCATCATAACCTGAGTCTTGGCGTGGCGAACTCTATCGCGGCGGGGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGgc                                                                                >  1:112298/1‑251 (MQ=255)
                                                                                                                       gtggATTCTGGCGGTGCGATGAACATGAGCGATATCCGTGACCGTTTCCGCGCCCTGAAAGCAGAGCTGAAACCAGAAACGCAAACTGGCATGCACGCTCATCATAACCTGAGTCTtg                                                                                                                                                                                                                    <  1:10590/118‑1 (MQ=255)
                                                                                                                                        gATGAACATGAGCGATATCCGTGACCGTTTCCGCGCCCTGAAAGCAGAGCTGAAACCAGAAACGCAAACTGGCATGCACGCTCATCATAACCTGAGTCTTGGCGTGGCGAACTCTATCGCGGCGGTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATctct                                                              >  1:501092/1‑251 (MQ=255)
                                                                                                                                                            gTGACCGTTTCCGCGCCCTGAAAGCAGAGCTGAAACCAGAAACGCAAACTGGCATGCACGCTCATCATAACCTGAGTCTTGGCGTGGCGAACTCTATCGCGGCGGTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTgg                                                                                 <  1:359419/212‑1 (MQ=255)
                                                                                                                                                                ccGTTTCCGCGCCCTGAAAGCAGAGCTGAAACCAGAAACGCAAACTGGCATGCACGCTCATCATAACCTGAGTCTTGGCGTGGCGAACTCTATCGCGGCGGTGGNAGANGGCGNCGNCCGAATCGACGCCAGNCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTNNTNNNTNCCGCcgcg                                                                                                <  2:345212/193‑1 (MQ=255)
                                                                                                                                                                ccGTTTCCGCGCCCTGAAAGCAGAGCTGAAACCAGAAACCCAAACTGGCATGCACGCTCATCATAACCTGAGTCTTGGCGTGGCGAACTCTATCGCGGCGGTGGAAGAGGGGTGCGACCGAATCGACGCCAACCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACAGCTGGAAGTGTTTATTGCCGCcgcg                                                                                                >  1:345212/1‑193 (MQ=255)
                                                                                                                                                                      ccGCGCCCTGAAAGCAGAGCTGAAACCAGAAACGCAAACTGGCATGCCCGCTCATCATAACCTGAGTCTTGGCGTGGCGAACTCTATCGCGGCGGTGGAAGAGGGCGGCGACCGAAGCGACGCCAGCCTCGCGGGAATGGGCGCGGGTGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCGGCGGATAAACTGGGCGGTCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTgg                                <  1:328085/251‑1 (MQ=255)
                                                                                                                                                                           cccTGAAAGCAGAGCTGAAACCAGAAACGCAAACTGGCATGCACGCTCATCATAACCTGAGTCTTGGCGTGGCGAACTCTATCGCGGCGGCGGCAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGCCGACCTGGTGCGt                           >  1:452759/1‑251 (MQ=255)
                                                                                                                                                                             cTGAAAGCAGAGCTGAAACCAGAAACGCAAACTGGCATGCACGCTCATCCTAACCTGAGTCTTGGCGTGGCGAACTCTCTcgc                                                                                                                                                                                                 >  1:244383/1‑83 (MQ=255)
                                                                                                                                                                                        gcnaaanCNNGNNGCGNAAACGGGCATGCACGCACATCNNAACCTGAGTNTTGGCGTGGCGAACTCTATCGCGGCGGTGGAAGAGGGCTGCGACCGAATCGACNCCAGCCTCGNNGGNATGGTNGNNNGNNNAGNTAACGCACCGCTGGAAGTGNTTNNTGCCGCNGCGGATAAACTGGGCTTTCAGCATgg                                                                         <  2:112298/188‑1 (MQ=255)
                                                                                                                                                                                                    aaCGCAAACTGGCATGCACGCTCATCATAACCTGAGTCTTGGCGTGGCGAACTCTATCGCGGCGGTGGAAGAGGGCTGCGACCGAATCGGCGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGTTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAg  <  1:218637/251‑1 (MQ=255)
                                                                                                                                                                                                               gCATGCACGCTCATCATAACCTGAGTCTTGGCGTGGCGCACTCTATCGCGCCGGTGGAAGAGGGCTGCGACCGAATAGACGCCAGCCTCGCGGCAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATctct                                                              >  1:466519/1‑180 (MQ=255)
                                                                                                                                                                                                                            |                                                                                                                                                                                                                                  
CCCGCGAGCTCGGAATGGACACCGTTGGTTTTCTGATGATGAGCCATATGACCACGCCGGAGAATCTCGCCAAGCAGGCAAAGCTGATGGAAGGCTACGGTGCGACCTGTATTTATGTGGTGGATTCTGGCGGTGCGATGAACATGAGCGATATCCGTGACCGTTTCCGCGCCCTGAAAGCAGAGCTGAAACCAGAAACGCAAACTGGCATGCACGCTCACCATAACCTGAGTCTTGGCGTGGCGAACTCTATCGCGGCGGTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAG  >  NC_000913/374244‑374690

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