Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,818,446	C→T	N229N (AAC→AAT)	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,446	0	C	T	96.3%	82.2 / ‑5.8	27	N229N (AAC→AAT)	dinD	DNA damage‑inducible protein D
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (9/17);  total (9/18)
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.97e-01

GGCAGCACAGCAAGCTGGGGTAGCGACAGCTACTGATTTCGCCATATTTCAGAATCATGGTTACCAGGGGCTGTATGGTGGATTAGATCAGAAAGCTATTCATCAGCGGAAGGGGCTGAAAAAGAATCAGAAGATCCTGGATCATATGGGTTCAACAGAACTGGCGGCTAATCTCTTTCGAGCTACCCAAACAGAAGAAAAACTCAAGCGGGATGGCGTGAATTCAAAACAGCAAGCAAACACCACACACTTTGACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAGCAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCT  >  NC_000913/3818206‑3818688                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   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)                                                                   ggggCTGTATGGTGGATTAGATCAGAAAGCCATCCATCAGCGGAAGGGGCTGAAAAAGAGTCAGAAGATCCTGGATCATATGGGTTCAACAGAACTGGCGGCGAATCTCTTTCGAGCTACCCAAACAGAAGAAAAACTCAAGCGGGATGGCGTGAATTCAAAACAGCAAGCTAATACCACACACTTTGACGTGGGTAGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGc                                                                                                                                                                        <  1:360097/251‑1 (MQ=255)                                                                                     aGATCAGAAAGCCATCCATCAGCGGAAGGGGCTGAACAAGAGTCAGAAGATCCTGGATCATATGGGTTCAACAGAACTGGCGGCGAATCTCTTTCGAGCTACCCAAACAGAAGAAAAACTCAAGCGGGCTGGCGTGAATTCAAAACAGCAAGCTAATACCACACACTTTGAc                                                                                                                                                                                                                                     >  2:296326/1‑172 (MQ=255)                                                                                     aGATCAGAAAGCCATCCATCAGCGGAAGGGGCTGAAAAAGAGTCAGAAGATCCTGGATCATATGGGTTCAACAGAACTGGCGGCGAATCTCTTTCGAGCTACCCAAACAGAAGAAAAACTCAAGCGGGATGGCGTGAATTCAAAACAGCAAGCTAATACCACACACTTTGAc                                                                                                                                                                                                                                     <  1:296326/172‑1 (MQ=255)                                                                                       aTCAGAAAGCGATCCATCAGCGGAAGGTGCTGAAAAAGAGTCATAAGATCCTGGATCATATGGGTTCAACAGAACTGGCGGCGAATCTCTGTCGAGCTACCCAAACATAACAAAAACTCAAGCGGGATGGCGTGAATGCAAAACAGCAAGCTAATACCACACACGGTGACGTGTGTAGCAAGGTGAGGCAAACCAATCAGGAACTTGGCGGAATGATGCCTTAGGAATTGCCGACCCCGGAGGgcagcagc                                                                                                                                                    <  2:249945/251‑3 (MQ=255)                                                                                               gccatccatcaGCGGAAGGGGCTGAAACAGAGTCAGAAGATCCTGGATCATATGGGTTCAACAGAACTGGCGGCGAATCTCTTTCGAGCTACCCAAACAGAAGAAAAACTCAAGCGGGATGGCGTGAATTCAAAACAGCAAGCTAATACCACACACTTTGACGTGGGTAGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATg                                                                                                                                                                                     >  2:21979/7‑210 (MQ=255)                                                                                               gccatccatcaGCGGAAGGGGCTGAAAAAGAGTCAGAAGATCCTGGATCATATGGGTTCAACAGAACTGGCGGCGAATCTCTTTCGAGCTACCCAAACAGAAGAAAAACTCAAGCGGGATGGCGTGAATTCAAAACAGCAAGCTAATACCACACACTTTGACGTGGGTAGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATg                                                                                                                                                                                     <  1:21979/204‑1 (MQ=255)                                                                                                  atccatcaGCGGAAGGGGCTGAAAAAGAGTCAGAAGATCCTGGATCATATGGGTTCAACAGAACTGGCGGCGAATCTCTTTCGAGCTACCCAAACAGAAGAAAAACTCAAGCGGGATGGCGTGAATTCAAAACAGCAAGCTAATACCACACACTTTGACGTGGGTAGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCATGAGGAGTTGCCGACCCCGCAGGTCAGCGTCAAGCAGTTgg                                                                                                                                          >  1:22960/4‑250 (MQ=255)                                                                                                            ggAAGGGGCTGAAAAAGAGTCAGAAGATCCTGGATCATATGGGTTCAACAGAACTGGCGGCGAATCTCTTTCGAGCTACCCAAAAAGAAGAAAAACTCAAGCGGGATGGCGTGAATTCAAAACAGCAAGCTAATACCACACACTTTGACGTGGGTAGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAgtgt                                                                                                                                >  1:126550/1‑250 (MQ=255)                                                                                                             gAAGGGGCTGAAAAAGAGTCAGAAGATCCTGGATCATATGGGTTCAACAGAACTGGCGGCGAATCTCTTTCGAGCTACCCAAACAGAAGAAAAACTCAAGCGGGATGGCGTGAATTCAAAACAGCAAGCTAATACCACACACTTTGACGTGGGTAGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTaa                                                                                                                              <  1:432086/251‑1 (MQ=255)                                                                                                                    cTGAAAAAGAGTCAGAAGATCCTGGATCATATGGGTTCAACAGAACTGGCGGCGCATCTCTTTCGAGCTACCCAAACAGAAGAAAAACTCAAGCGGGATGGCGTGAATTCAAAACAGCAAGCTAATACCACACACTTTGACGTGGGTAGCAAGGTGAGGc                                                                                                                                                                                                                  >  2:304305/1‑160 (MQ=255)                                                                                                                    cTGAAAAAGAGTCAGAAGATCCTGGATCATATGGGTTCAACAGAACTGGCGGCGAATCTCTTTCGAGCTACCCAAACAGAAGAAAAACTCAAGCGGGATGGCGTGAATTCAAAACAGCAAGCTAATACCACACACTTTGACGTGGGTAGCAAGGTGAGGc                                                                                                                                                                                                                  <  1:304305/160‑1 (MQ=255)                                                                                                                       aaaaaGAGTCAGAAGATCCTGGATCATATGGGTTCAACAGAACTGGCGGCGAATCTCTTTCGAGCTACCCAAACAGAAGAAAAACTCAAGCGGGATGGCGTGAATTCAAAACAGCGAGCTAATACCACACACTTTGACGTGGGTAGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACaga                                                                                                                    <  1:108899/251‑1 (MQ=255)                                                                                                                                cAGAAGATCCTGGATCATATGGGTTCAACAGAACTGGCGGCGAATCGCTTGCGAGCTACCCAAACAGAAGAAAATCTCAAGCGGGATGGCGTGAATTCAAAACAGCAAGCTAATACCACACACTTTGACGTGGGTAGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATa                                                                                                           <  1:68932/251‑1 (MQ=255)                                                                                                                                           ggATCATATGGGTTCAACAGAACTGGCGGCGAATCTCTTTCGAGCTACCCAAACAGAAGAAAAACTCAAGCGGGATGGCGTGAATTCAAAACAGCAAGCTAATACCACACACTTTGACGTGGGTAGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAGCAGAAACt                                                                                                  <  2:108379/249‑1 (MQ=255)                                                                                                                                                    gggTTCAACAGAACTGGCGGCGAATCTCTTTCGAGCTACCCAAACAGAAGAAAAACTCAAGCGGGATGGCGTGAATTCAAAACAGCAAGCTAATACCACACACTTTGACGTGGGTAGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAGCAGAAACTAAAAACTCttt                                                                                       <  1:137511/251‑1 (MQ=255)                                                                                                                                                                  tGGCGGCGAATCTCTTTCGAGCTACCCAAACAGAAGAAAAACTCAAGCGGGATGGCGTGAATTCAAAACAGCAAGCTAATACCACACACTTTGACGTGGGTAGCAAGGTGAGGCAACCCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAGCAGAAACTAAAAACTCTTTTTGttgattg                                                                             >  2:100846/1‑247 (MQ=255)                                                                                                                                                                                    gAGCTACCCAAACAGAAGAAAAACTCAAGCGGGATGGCGTGAATTCAAAACAGCAAGCTAATACCACACACTTTGACGTGGGTAGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAGCAGAAACTAAAAACTCTTTTTGTTGATTGAGATACCCGATGCGTAAGCCaa                                                       <  1:107762/251‑1 (MQ=255)                                                                                                                                                                                                   aagaaAAACTCAAGCGGGATGGCGTGAATTCAAAACAGCAAGCTAATACCACACACGTTGACGTGGGTAGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAGCAGAAACTAAAAACTCTTTTTGTTGATTGAGATACCCGATGCGTAAGCCAAAGTCCAGGTGCGGTg                                        <  1:233434/251‑1 (MQ=255)                                                                                                                                                                                                             cAAGCGGGATGGCGTGAATTCAAAACAGCAAGCTAATACCACACACTTTGACGTGGGTAGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGTAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAGCAGAAACTAAAAACTCTTTTTGTGGATTGAGATACCCGATGCGTAAGCCAAAGTCCAGGTGCGGTGAACACAATg                               <  2:16846/250‑1 (MQ=255)                                                                                                                                                                                                                gCGGGATGGCGTGAATTCAAAACAGCAAGCTAATACCACACACTTTGACGTGGGTAGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAGCAGAAACTAAAAACTCTTTTTGTTGATTGAGATACCCGATGCGTAAGCCAAAGTCCAGGTGCGGTGAACACAATGGCTa                           <  2:315286/251‑1 (MQ=255)                                                                                                                                                                                                                                 cAAAACAGCAAGCTAATACCACACACTTTGACGTGGGTAGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCAGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAGCAGAAACTAAAAACTCTTTTTGTTGATTGAGATACCCGATGCGTAAGCCAAAGTCCAGGTGCGGTGAACACAATGGCTAAACATTGCACctgtct           <  2:244652/250‑1 (MQ=255)                                                                                                                                                                                                                                 cAAAACAGCAAGCTAATACCACACACTTTGACGTGGGTAGCAAGGTGAGGCAAACCATTCAGGAACTTAGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAGCAGAAACTAAAAACTCTTTTTGTTGATTGAGATACCCGATGCGTAAGCCAAAGTCCAGGTGCGGTGAACACAATGGCTAAACATTGCACctttctt          >  1:415977/1‑251 (MQ=255)                                                                                                                                                                                                                                         caagctaATACCACACACTTTGACGTGGGTAGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAGCAGAAGCTAAAAACTCTTTTTGTTGATTGAGATACCCGATGCGTAAGCCAAAGGCCATGTGCGGTGAACACAATGGCTAAACATTGCACCTTTCTTTTccaccc  <  1:392893/251‑5 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   GGCAGCACAGCAAGCTGGGGTAGCGACAGCTACTGATTTCGCCATATTTCAGAATCATGGTTACCAGGGGCTGTATGGTGGATTAGATCAGAAAGCTATTCATCAGCGGAAGGGGCTGAAAAAGAATCAGAAGATCCTGGATCATATGGGTTCAACAGAACTGGCGGCTAATCTCTTTCGAGCTACCCAAACAGAAGAAAAACTCAAGCGGGATGGCGTGAATTCAAAACAGCAAGCAAACACCACACACTTTGACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAGCAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCT  >  NC_000913/3818206‑3818688

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