Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 248,882 C→T intergenic (+748/‑2016) rayT → / → dinB REP‑associated tyrosine transposase/DNA polymerase IV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913248,8820CT100.0% 16.3 / NA 7intergenic (+748/‑2016)rayT/dinBREP‑associated tyrosine transposase/DNA polymerase IV
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (4/3);  total (4/3)

TCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGATCACGCTGGCGCTGTCGATCACCTGATACCCCATATTCAGCGCCTTCGCTTTCTGCGCCGCCTGGATCCAGGTGACC  >  NC_000913/248641‑249069
                                                                                                                                                                                                                                                 |                                                                                                                                                                                           
tCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTGCGTCACCGTGCTACTGGCGACCAGCACGGTGGGGATGGTGACGATATCGCGCTGGTAAACGCCTTCGGTCAGCAGCGTCCGGTACACTTTCAGCATCGGTCTGTAATTGAGCGCCGCGCTCAAATCTTTCGCCAGan                                                                                                                                                                                     <  1:19620/250‑2 (MQ=255)
                                                                                 ngcACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTTCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAAt                                                                                                   >  1:103569/2‑251 (MQ=255)
                                                                                             nccAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTTCGCCAGACGCGGAGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATg                                                                                       >  1:41873/2‑251 (MQ=255)
                                                                                               cagcagAATATGATCTTTCGTCACCGNGCTACTNNCGNCCAGCACGGTGGCGATGGTNACNNNANNNCNCAGGGNGANNCCTTCGGTCNGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAANCTTTCGCCANNCGCGGTGCCGTCGACGATAAANGGTNNTNNANCTGNGtn                                                                                                                                               >  2:38396/1‑192 (MQ=255)
                                                                                                                                   gACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAGACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTTCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGATCACGCTGGCGCTGTCGATCACCTGATAcc                                                 <  1:38396/251‑1 (MQ=255)
                                                                                                                                                                cgcgCANNGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTTCGCCAGACGCGGtc                                                                                                                                                                               <  2:8530/96‑2 (MQ=255)
                                                                                                                                                                                  nggTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTTCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGATCACGCTGGCGCTGTCGATCACCTGATACCCCATATTCAGCGCCTTCGCTTTCTGCGCCGCCTGGATCCAGGTGAcc  >  1:103240/2‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                           
TCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGATCACGCTGGCGCTGTCGATCACCTGATACCCCATATTCAGCGCCTTCGCTTTCTGCGCCGCCTGGATCCAGGTGACC  >  NC_000913/248641‑249069

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: