Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,575,050 C→T intergenic (+1119/+1862) yjiT → / ← mcrC putative uncharacterized protein YjiT/Type IV methyl‑directed restriction enzyme EcoKMcrBC subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,575,0500CT87.5% 15.9 / ‑6.1 8intergenic (+1119/+1862)yjiT/mcrCputative uncharacterized protein YjiT/Type IV methyl‑directed restriction enzyme EcoKMcrBC subunit
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (3/4);  total (4/4)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGCCAATGCCGATCTGGAAAAATTAGCGGCCAGGGCGAAGCAAGCCCGTGAGCAGGCAGAAGAATATCGAAGCTGGGGAATGCCAGACCAGGCGAAATTGTCACAAGCACAGGCTGAACAGCTTGAACAGGCTTATCAGGCTGCAACCGGTGGGAAAGCCGCGACTATCCTGGTATCCCGAACCTGGACGGAGATGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCAGGAAAAAGGATTAACGCTGGATGACTG  >  NC_000913/4574817‑4575273
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                               
atgccaatgccGATCTGGAAAAATTAGCGGCCAGGGCGAAGCAAGCCCGTGAGCAGGCAGAAGAATATCGAAGCTGGGGAATGCCAGACCAGGCGAAATTGTCACAAGCACAGGCTGAACAGCTTGAACAGGCTTATCAGGCGGCAACCGGTGTGAAAGCCGCGACTATCCTGGTATCCCGAGCCTGGACGGATATGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGTCGGTTCTGCAAtatacc                                                                                                                                                                                                                <  1:203946/251‑3 (MQ=255)
                              ccAGGGCGAAGCAAGCCCGTGAGCAGGCAGAAGAATATCGAAGCTGGGGAATGCCAGACCAGGCGACATTGTCACAAGCACAGGCTGAACAGCTTGAACAGGCTTATCAGGCTGCAACCGGTGGGAAAGCCGCGACTATCCTGGTATCCCGAACCTGGACGGAGATGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGTCGGTTCTGCAATATAACCATTa                                                                                                                                                                                                           >  1:153821/1‑226 (MQ=255)
                                    cGAAGCAAGCCCGTGAGCAGGCAGAAGAATATCGAAGCTGGGGAATGCCAGACCAGGCGAAATTGTCACAAGCACAGGCTGAACAGCTTGAACAGGCTTATCAGGCTGCAACCGGTGGGAAAGCCGCGACTATCCTGGTATCCCGAACCTTGACGGAGATGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGTCGGTTCt                                                                                                                                                                                                                          <  1:189203/205‑1 (MQ=255)
                                                                        gCTGGGGAATGCCAGACCAGGCGAAATTGTCACAAGCACAGGCTGAACAGCTTGAACAGGCTTATCAGGCTGCAACCGGTGGGAAAGCCGCGACTATCCTGGTATCCCGAACCTGGACGGAGATGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGTCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCgg                                                                                                                                        <  1:176830/251‑1 (MQ=255)
                                                                                      gACCAGGCGAAATTGTCACAAGCACAGGCTGAACAGCTTGAACAGGCTTATCAGGCTGCAACCGGTGGGAAAGCAGCGACAATCCTGGTATCCCGAACCTGGACGGAGATGGATAACGAGATTAAAGAGAGAGCCGATATACGCGGGCCGGTTCTGCAATCTAACCATTATCTTAAGCATGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCg                                                                                                                          >  1:151900/1‑251 (MQ=255)
                                                                                                                              aaCAGGCTTATCAGGCTGCAACCGGTGGGAAAGCCGCGACTATCCTGGTATCCCGAACCTGGACGGAGATGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGTCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGtt                                                                                  >  1:177087/1‑251 (MQ=255)
                                                                                                                                                                                                              nttAAAGAGAGAGCCGATATCCGCGGGTCGGTTCTGCAATATAACCATTCTCTTAAGACTGAAGAGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGCAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTACCAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCAGGAAAAAGGAGTAACGCTGGATGACTg  >  1:20412/2‑251 (MQ=255)
                                                                                                                                                                                                                tAAAGAGAGATCCGATATCCGCGGGTCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCAGGAAAATTGATTATCGCTGGATGACTg  <  1:167562/249‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                               
ATGCCAATGCCGATCTGGAAAAATTAGCGGCCAGGGCGAAGCAAGCCCGTGAGCAGGCAGAAGAATATCGAAGCTGGGGAATGCCAGACCAGGCGAAATTGTCACAAGCACAGGCTGAACAGCTTGAACAGGCTTATCAGGCTGCAACCGGTGGGAAAGCCGCGACTATCCTGGTATCCCGAACCTGGACGGAGATGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCAGGAAAAAGGATTAACGCTGGATGACTG  >  NC_000913/4574817‑4575273

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