Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,818,128	A→T	A123A (GCA→GCT)	dinD →	DNA damage‑inducible protein D

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,818,128	0	A	T	100.0%	48.3 / NA	16	A123A (GCA→GCT)	dinD	DNA damage‑inducible protein D
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (9/7);  total (9/7)

AAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGTGAAGACGAAAAGCGTCTGTTTCTGAGAAATGAGTTGAAAGAACATAACAAACAATTGGTTGAGGCAGCACAGCAAGCTGGGGTAGCGACAGCTACTGATTTCGCCATATTTCAGAATCATGGTTACCAGGGGCTGTATGGTGGATTAGATCAGAAAGCTATTCATCAGCGGAAGGGGCTGAAAAAGAATCAGAAGATCCTGGATCATATGGGTTCAACAGAACT  >  NC_000913/3817982‑3818367                                                                                                                                                   |                                                                                                                                                                                                                                                caacgagaacttgaagatgttcatctctCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTTATTGCGGCAGGGCAAACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTGGCTGACGATGAAGCTTTCAGGCAACTTCGTGAAGACGAAAAGCGTCTTTTTCTGAGAAATGAATTGAAAGAACATAACAAACAATTGGTTGAGGCAGCGCAGCAAGCTgg                                                                                                                                                  <  1:307198/218‑1 (MQ=255) caacgagaacttgaagatgttcatctctCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTTATTGCGGCAGGGCAAACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTGGCTGACGATGAAGCTTTCAGGCAACTTCGTGAAGACGAAAAGCGTCTTTTTCTGAGAAATGAATTGAAAGAACATAACAAACAATTGGTTGAGGCAGCGCAGCAAGCTgg                                                                                                                                                  >  2:307198/25‑242 (MQ=255)                   gttcatctctCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTTATTGCGGCAGGGCAAACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTGGCTGACGATGAAGCTTTCAGGCAACTTCGTGAAGACGAAAAGCGTCTTTTTCTGAGAAATGAATTGAAAGAACATAACAAACAATTGGTTGATGCAGCGCAGCAAGCTGGGGGAGCGACAGCTATCGATTTCGCCat                                                                                                                       <  1:279019/245‑1 (MQ=255)                   ggtcatctctCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTTATTGCGGCAGGGCAAACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTGGCTGACGATGAAGCTTTCAGGCAACTTCGTGAAGACGAAAAGCGTCTTTTTCTGAGAAATGAATTGAAAGAACATAACAAACAATTGGTTGAGGCAGCGCAGCAAGCTGGGGGAGCGACAGCTATCGATTTCGCCa                                                                                                                        <  2:76138/244‑1 (MQ=255)                         ctctCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTTATTGCGGCAGGGCAAACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTGGCTGACGATGAAGCTTTCAGGCAACTTCGTGAAGACGAAAAGCGTCTTTTTCTGAGAAATGAATTGAAAGAACATAACAAACAATTGGTTGAGGCAGCGCAGCAAGCTGGGGGAGCGACAGCTATCGATTTCGCCATATTTc                                                                                                                  >  1:244652/1‑250 (MQ=255)                               cGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTTATTGCGGCAGGGCAAACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTGGCTGACGATGAAGCTTTCAGGCAACTTCGTGAAGACGAAAAGCGTCTTTTTCTGAGAAATGAATTGAAAGAACATAAAAAACAATTGGTTGAGGCAGCGCAGAAAGCTGGGGGAGCGACAGCTATCGATTTCGCCATATTTCAGAATc                                                                                                            >  2:432086/1‑250 (MQ=255)                                        tGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTTATTGCGGCAGGGCAAACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTGGCTGACGATGAAGCTTTCAGGCAACTTCGTGAAGACGAAAAGCGTCTTTTTCTGAGAAATGAATTGAAAGAACATAACAAACAATTGGTTGAGGCAGCGCAGCAAGCTGGGGGAGCGACAGCTATCGATTTCGCCATATTTCAGAATCATGGTTAccg                                                                                                  >  2:107762/1‑250 (MQ=255)                                          ttACCTGGTAGTACAAAACGGCGACCCTGAGAAACCGGTTATTGCGGCAGGGCAAACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTGGCTGACGATGAAGCTTTCAGGCAACTTCGTGAAGACGAAAAGCGTCTTTTTCTGAGAAATGAATTGAAAGAACATAACAAACAATTGGTTGAGGCAGCGCAGCAAGCTGGGGGAGCGACAGCTATCGATTTCGCCATATTTCAGAATCATGGTTAccggg                                                                                                >  1:97317/1‑248 (MQ=255)                                                                                 taTTGCGGCATGGCAAACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTGGCTTACGATGAAGCTTTCAGGCAATTTCGTGAAGACGAAAAGCGTCTTTTTCTGAGAACTGAATTGAAAGAACATAACAAACAATTGGTTGAGGCATCGCAGCAAGCTGGTGGAGCGACAGCTATCGATTTCGCCATATTTCAGAATCATGGTTACCGGGGGCTGTATGGTGGATTAGATCAGAAAGCCATCCATCAg                                                          <  1:231213/249‑1 (MQ=255)                                                                                 taTTGCGGCAGGGCAAACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTGACTGACGATGAAGCTTTCAGGCAACTTCGTGAAGACGAAAAGCGTCTTTTTCTGAGAAATGAATTGAAAGAACATAACAAACAATTGGTTGAGGCAGCGCAGCAAGATGAAGGAGCGCCAGCTATCGATTACGACATATTTCAGAATCATGGTTACAGGGGGCTGTATGGAGGATTAGATCAGAAAGCGATCCATCAGc                                                         >  1:287923/2‑251 (MQ=255)                                                                                            ggcaaACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTGGCTGACGATGAAGCTTTCAGGCAACTTCGTGAAGACGAAAAGCGTCTTTTTCTGAGAAATGAATTGAAAGAACATAACAAACAATTGGTTGAGGCAGCGCAGCAAGCTGGGGGAGCGACAGCTATCGATTTCGCCATATTTCAGAATCATGGTTACCGGGGGCTGTATGGTGGATTAGATCAGAAAGCCATCCATCAGCGGAAGGGGCTg                                              >  1:108379/1‑251 (MQ=255)                                                                                                       tttGCTAGCCAGACCCGACGGCAGGAGCTGGCTGACGATGAAGCTTTCAGGCAACTTCGTGAAGACGAAAAGCGTCTTTTTCTGAGAAATGAATTGAAAGAACATAACAAACAATTGGTTGAGGCAGCGCAGCAAGCTGGGGGAGCGACAGCTATCGATTGCGCCATATTTCAGAATCATGGTTACCGGGGGCTGTATGGTGGATTAGATCAGAAAGCCATCCATCAGCGGAAGGGGCTGAAAAAGagtc                                    <  2:109449/250‑4 (MQ=255)                                                                                                           cTATCCAGACCCGACGGCAGGAGCTGGCTGACGATGAAGCTTTCAGGCAACTTCGTGAAGACGAAAAGCGTCTTTTTCTGAGAAATGAATTGAAAGAACATAACAAACAATTGGTTGAGGCAGCGCAGCAAGCTGGGGGAGCGACAGCTATCGATTTCGCCATATTTCAGAATCATGGTTACCGGGGGCTGTATGGTGGATTAGATCAGAAAGCCATCCATCAGCGGAAGGGGCTGAAAAAGAGTCaga                                 >  2:137511/1‑249 (MQ=255)                                                                                                                    cccGACGGCAGGAGCTGGCGGACGATGAAGCTTTCAGGCAACTTCGTGAAGACGAAAAGCGTCTTTTTCTGAGAAATGAATTGAAAGAACATAACAAACAATTGGTTGAGGCAGCGCAGCAAGCTGGGGGAGCGACAGCTATCGATTTCGCCATATTTCAGAATCATGGTTACCGGGGGCTGTATGGTGGATTAGATCAGAAAGCCATCCATCAGCGGAATGGGCTGAAAAAGAGTCAGAAGAGCCTGGAt                      <  2:161031/251‑1 (MQ=255)                                                                                                                                 gctggCTGACGATGAAGCTTTCAGGCAACTTCGTGAAGACGAAAAGCGTCTTTTTCTGAGAAATGAATTGAAAGAACATAACAAACAATTGGTTGAGGCAGCGCAGCAAGCTGGGGGAGCGACAGCTATCGATTTCGCCATATTTCAGAATCATGGTTACCGGGGGCTGTATGGTGGATTAGATCAGAAAGCCATCCATCAGCGGAAGGGGCTGAAAAAGAGTCAGAAGATCCTGGATCATATGGGTTCaa         <  2:12647/247‑1 (MQ=255)                                                                                                                                        gACGATGAAGCTTTCAGGCAACTTCGTGAAGACGAAAAGCGTCTTTTTCTGAGAAATGACTTGAAAGAACATAACAAACAATTGGTTGAGGCAGCGCAGCAAGCTGGGGGAGCGCCAGCTATCGATTTCGCCATATTTCAGAATCATGGTTAACGGGGGCTGTATGGTGGATTAGATCAGAAAGCCATCCATCAGCGGAAGGGGCTGAAAAAGAGTCAGAAGATCCTGGATCATATGGGTTCAACAGAACt  >  1:72087/1‑251 (MQ=255)                                                                                                                                                   |                                                                                                                                                                                                                                                AAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGTGAAGACGAAAAGCGTCTGTTTCTGAGAAATGAGTTGAAAGAACATAACAAACAATTGGTTGAGGCAGCACAGCAAGCTGGGGTAGCGACAGCTACTGATTTCGCCATATTTCAGAATCATGGTTACCAGGGGCTGTATGGTGGATTAGATCAGAAAGCTATTCATCAGCGGAAGGGGCTGAAAAAGAATCAGAAGATCCTGGATCATATGGGTTCAACAGAACT  >  NC_000913/3817982‑3818367

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