Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,814,166 T→C R100R (CGT→CGC katE → catalase HPII, heme d‑containing

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,814,1660TC85.7% 34.0 / ‑0.7 14R100R (CGT→CGCkatEcatalase HPII, heme d‑containing
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (6/6);  total (7/7)
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

CTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTACGCAAAGGCAGTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCC  >  NC_000913/1814039‑1814285
                                                                                                                               |                                                                                                                       
cTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTACGCAAAGGCAGTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGCGGTCCAACgct                                                                                                              >  1:520656/1‑139 (MQ=255)
        cgTAACGAAAAACTTAATTCTCTGGAAGACGTACGCAAAGGCAGTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGCGGTCCAACGCTGCTGGAAg                                                                                                      >  2:339036/1‑139 (MQ=255)
                                 aaGACGTACGCAAAGGCAGTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCAc                                                                             <  2:490173/139‑1 (MQ=255)
                                                          ttATGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGCGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCg                                                    >  1:168517/1‑139 (MQ=255)
                                                                  tGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGCGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTAt                                            >  1:140466/1‑139 (MQ=255)
                                                                    aCCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTg                                          >  2:541589/1‑139 (MQ=255)
                                                                    aCCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGCGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCAtt                                                       >  1:40498/1‑126 (MQ=255)
                                                                    aCCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGCGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCAtt                                                       <  2:40498/126‑1 (MQ=255)
                                                                         tAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGCGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCAt                                     <  2:140466/139‑1 (MQ=255)
                                                                                      cgcATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGCGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAg                        <  2:168517/139‑1 (MQ=255)
                                                                                                          cTCACTGCGTGCCGGTAGCCGCGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAg    <  2:520656/139‑1 (MQ=255)
                                                                                                            cACTGCGTGCCGGTAGCCGCGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGcc  <  1:339036/139‑1 (MQ=255)
                                                                                                             aCTGCGTGCCGGTAGCCGCGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTcc      <  1:481506/134‑1 (MQ=255)
                                                                                                             aCTGCGTGCCGGTAGCCGCGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTcc      >  2:481506/1‑134 (MQ=255)
                                                                                                                               |                                                                                                                       
CTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTACGCAAAGGCAGTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCC  >  NC_000913/1814039‑1814285

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