Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,818,271	A→G	Q171R (CAG→CGG)	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,271	0	A	G	100.0%	59.4 / NA	20	Q171R (CAG→CGG)	dinD	DNA damage‑inducible protein D
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (9/11);  total (9/11)

CATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGTGAAGACGAAAAGCGTCTGTTTCTGAGAAATGAGTTGAAAGAACATAACAAACAATTGGTTGAGGCAGCACAGCAAGCTGGGGTAGCGACAGCTACTGATTTCGCCATATTTCAGAATCATGGTTACCAGGGGCTGTATGGTGGATTAGATCAGAAAGCTATTCATCAGCGGAAGGGGCTGAAAAAGAATCAGAAGATCCTGGATCATATGGGTTCAACAGAACTGGCGGCTAATCTCTTTCGAGCTACCCAAACAGAAGAAAAACTCAAGCGGGATGGCGTGAATTCAAAACAGCAAGCAAACACCACACACTTTGACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAAC  >  NC_000913/3818062‑3818495                                                                                                                                                                                                                  |                                                                                                                                                                                                                                 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)                                                                       aggCAACTTCGTGAAGACGAAAATCGTCTTTTTCTGAGAAATGAATTGAAAGAACATAACAAACAATTGGTTGAGGCAGCGCAGCAAGCTGGGGGAGCGACAGCTGTCGATTTTGCCATATTTCAGAATCATGGTTACCGGGGGCTGTATGGTGGGTGAGATCAGAAAGCCATCCATCAGCGGAGGGGGCTGAAAAAGAGTCAGAAGATCCTGGATCATATGGGTTCAACAGAACTGGCGGCGAATCTCtt                                                                                                                   <  1:290851/249‑1 (MQ=255)                                                                          cAACTTCGTGAAGACGAAAAGCGTCTTTTTCTGAGAAATGAATTGAAAGAACATAACAAACAATTGGTTGAGGCAGCGCAGCAAGCTGGGGGAGCGACAGCTATCGATTTCGCCATATTTCAGAATCATGGTTACCGGGGGCTGTATGGTGGATTAGATCAGAAAGCCATCCATCAGCGGAAGGGGCTGAAAAAGAGTCAGAAGATCCTGGATCATATGGGTTCAACAGAACTGGCGGCGAATCTCTTTCg                                                                                                                >  2:360097/1‑251 (MQ=255)                                                                                   gAAGACGAAAAGCGTCTTTTTCTGAGAAATGAATTGAAAGAACATAACAAACAAGTGGTTGAGGCAGCGCAGCAAGCTGGGGGAGCGACAGCTATCGATTTCGCCATATTTCAGAATCATGGTTACCGGGGGCTGTATGGTGGATTAGATCAGAAAGCCATCCATCAGCGGAAGGGGCTGAAAAAGAGTCAGAAGATCCTGGATCATATGGGTTCAACAGAACTGGCGGCGAATCTCTTTCGAGCTcccc                                                                                                        <  2:97317/250‑1 (MQ=255)                                                                                          aaaaGCGTCTTTTTCTGAGAAATGAATTGAAAGAACATAACAAACAATTGGTTGAGGCAGCGCAGCAAGCTGGGGGAGCGACAGCTATCGATTTCGCCATATTTCAGAATCATGGTTACCGGGGGCTGTATGGTGGATTAGATCAGAAAGCAATCCATCAGCGGAAGGGGCTGAAAAAGAGTCAGAAGATCCTGGATCATATGGGTTCAACAGCACTGGCGGCGAATCTCTTTCGAGCTACCCAAACaga                                                                                                 >  1:16846/1‑250 (MQ=255)                                                                                             aGCGTCTTTTTCTGAGAAATGAATTGAAAGAACATAACAAACAATTGGTTGAGGCAGCGCAGCAAGCTGGGGGAGCGACAGCTATCGATTTCGCCATATTTCAGAATCATGGTTACCGGGGGCTGTATGGTGGATTAGATCAGAAAGCCATCCATCAGCGGAAGGGGCTGAAAAAGAGTCAGAAGATCCTGGATCATATGGGTTCAACAGAACTGGCGGCGAATCTCTTTCGAGCTACCCAAACAgaagaa                                                                                             <  1:193388/251‑1 (MQ=255)                                                                                                gTCTTTTTCTGAGAAATGAATTGAAAGAACATAACAAACAATTGGTTGAGGCAGCGCAGCAAGCTGGGGGAGCGACAGCTATCGATTTCGCCATATTTCAGAATCATGGTTACCGGGGGCTGTATGGTGGATTAGATCAGAAAGCCATCCATCAGCGGAAGGGGCTGAAAAAGAGTCAGAAGATCCTGGATCATATGGGTTCAACAGAACTGGCGGCGAATCTCTTTCGAGCTACCCAAACAGAAGaaaa                                                                                           >  1:315286/1‑250 (MQ=255)                                                                                                                      gAAAGAACATAACATACACGTGGGTGGGGCAGCGCAGCAAGCGGGGGGAGCTACAGCGAGCGATTTCGCCATAGTTCAGAATCATGGTTACCGGGGGCTGTATGGTGGATTTGATCAGAAAGCCATCCATCAGCGGAAGGGGCTGAAAAAGAGTCAGAAGATCCTGGATCATATGGGTGCAACAGAACTGGCGGCTAATCGCTTTCGAGCGACCCAAACAGAAGAAAAACTCAAGCGGGATGGCgggtat                                                                     <  2:72087/250‑6 (MQ=255)                                                                                                                            aCATAACAAACAATTGGTTGAGGCAGCGCAGCAAGCTGGGGGAGCGACAGCTATCGATTTCGCCATATTTCAGAATCATGGTTACCGGGGGCTGTATGGTGGATTAGATCAGAAAGCCATCCATCAGCGGAAGGGGCTGAAAAAGAGTCAGAAGATCCTGGATCATATGGGTTCAACAGAACTGGCGGCGAATCTCTTTCGAGCTACCCAAACAGAAGAAAAACTCAAGCGGGATGGCGTGAAGTCAAAAc                                                              <  2:54891/251‑1 (MQ=255)                                                                                                                                                 gtcagcgcagcaAGCTGGGGGAGCGACAGCTATCGATTTCGCCATATTTCAGAATCATGGTTAGCGGGGGCTGTATGGTGGATTAGATCAGAAAGCCATCCATCAGCGGAAGGGGCTGAAAAAGAGTCAGAAGATCCTGGATCATATGGGTTCAACAGAACTGGCGGCGAATCTCTTTCGAGCTACCCAAACAGAAGAAAAACTCAAGCGGGATGGCGTGAATTCAAAACAGCAAGCTAACACCACACACt                                         <  2:241655/249‑1 (MQ=255)                                                                                                                                                      cgcagcaAGCTGGGGGAGCGACAGCTATCGATTTCGCCATATTGCAGAATCATGGTTACCGGGGGCTGTATGGTGGATTAGATCAGAAAGCCATCCATCAGCGGAAGGGGCTGAAAAAGAGTCAGAAGATCCTGGATCATATGGGTTCAACAGAACTGGCGGCGAATCTCTTTCGAGCTACCCAAACAGAAGAAAAACTCAAGCGGGATGGCGTGAATTCAAAACAGCAAGCTAATACCACACACTTTGAc                                    <  1:217436/249‑1 (MQ=255)                                                                                                                                                                                  cgATTTCGCCATATTTCAGAATCATGGTTACCGGGGGCTGTATGGTGGATTAGATCAGAAAGCCATCCATCAGCGGAAGGGGCTGAAAAAGAGTCAGAAGATCCTGGATCATATGGGTTCAACAGAACTGGCGGCGAATCTCTTTCGAGCTACCCAAACAGAAGAAAAACTCAAGCGGGATGGCGTGAATTCAAAACAGCAAGCTAATACCACACACTTTGACGTGGGTAGCAAGGTGAGGCAAACCATTc        >  2:74712/2‑251 (MQ=255)                                                                                                                                                                                        cGCCATATTTCAGAATCTTGGTTACCGGGGGCTGTATGGCGGATTAGATCAGAAAGCCATCCATCAGCGGAAGGGGCTGAAAAAGAGTCAGCAGATCCTGGATCATATGGGTTCAACAGAACTGGCGGCGAATCTCTTTCGAGCTACCCAAACAGAAGAAAAACTCAAGCGGGATGGCGTGAATTCAAAACAGCAAGCTCATACCACCCACTTTGACGTGGGTAGCAAGGTGAGGCAAACCATTCAGGAAc  >  2:392893/1‑251 (MQ=255)                                                                                                                                                                                                                  |                                                                                                                                                                                                                                 CATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGTGAAGACGAAAAGCGTCTGTTTCTGAGAAATGAGTTGAAAGAACATAACAAACAATTGGTTGAGGCAGCACAGCAAGCTGGGGTAGCGACAGCTACTGATTTCGCCATATTTCAGAATCATGGTTACCAGGGGCTGTATGGTGGATTAGATCAGAAAGCTATTCATCAGCGGAAGGGGCTGAAAAAGAATCAGAAGATCCTGGATCATATGGGTTCAACAGAACTGGCGGCTAATCTCTTTCGAGCTACCCAAACAGAAGAAAAACTCAAGCGGGATGGCGTGAATTCAAAACAGCAAGCAAACACCACACACTTTGACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAAC  >  NC_000913/3818062‑3818495

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: ‑