Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,818,617	C→T	intergenic (+33/‑257)	dinD → / → yicG	DNA damage‑inducible protein D/PF03458 family inner membrane protein YicG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,818,617	0	C	T	100.0%	44.2 / NA	16	intergenic (+33/‑257)	dinD/yicG	DNA damage‑inducible protein D/PF03458 family inner membrane protein YicG
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (4/12);  total (4/12)

GAGCTACCCAAACAGAAGAAAAACTCAAGCGGGATGGCGTGAATTCAAAACAGCAAGCAAACACCACACACTTTGACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAGCAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTC‑‑CCCCTCAGTTTTAACCTATTTTTTCTTATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCTCTGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCCCTCCCCTTAAAAAATCTCAATCGTGACAATGCGCACAAATCGCTACCCTGCCA  >  NC_000913/3818385‑3818850                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                            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)                                                      caagctaATACCACACACTTTGACGTGGGTAGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAGCAGAAGCTAAAAACTCTTTTTGTTGATTGAGATACCCGATGCGTAAGCCAAAGGCCATGTGCGGTGAACACAATGGCTAAACATTGCACCTTTCTTTTc‑‑caccc                                                                                                                                                                    <  1:392893/251‑5 (MQ=255)                                                                ccacACACTTTGACGTGGGTAGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAGCAGAAACTAAAAACTCTTTTTGTTGATTGAGATACCCGATGCGTAAGCCAAAGTCCAGGTGCGGTGAACACAATGGCTAAACATTGCACCTTGCTTTTCCACCCCTCAATTTTa                                                                                                                                                            <  1:318767/251‑1 (MQ=255)                                                                                    agCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAGCAGAAACTAAAAACTCTTTTTGTTGATTGAGATACCCGATGCGTAAGCCAAAGTCCAGGTGCGGTGAACACAATGGCTAAACATTGCACCTTTCTTTTCCACCCCTCAATTTTAACCTATTTTTTCTTATGCAt                                                                                                                                        >  1:150038/2‑251 (MQ=255)                                                                                     gCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAGCAGAAACTAAAAACTCTTTTTGTTGATTGAGATACCCGATGCGTAAGCCAAAGTCCAGGTGCGGTGAACACAATGGCTAAACATTGCACCTTTCTTTTCCACCCCTCAATTTTAACCTATTTTTTCTTATGCAtt                                                                                                                                       <  1:74712/251‑1 (MQ=255)                                                                                                                                                    aGGTCAGCATCAAGCAGTTGTAACACAGTGTAAAAATTACAGAGAAGAAATAGCAGAAACTAAAAACTCTTTTTGTTGATTGAGATACCCTATGCGTAAGCCAAAGTCCAGGTGCGGTGAACACAATGGCTAAACATTGCACCTTTCTTTTCCACCCCTCAATTTTAACCTATTTTTTCTTATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCTACGCTACAGAAAATAGCGATTTCACATAACTa                                                                        <  1:211951/251‑1 (MQ=255)                                                                                                                                                              cAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAGCAGAAACTAAAAACTCTTTTTGTTGATTGAGATACCCGATGCGTAAGCCAAAGTCCAGGTGCGGTGAACACAATGGCTAAACATTGCACCTTTCTTTTCCACCCCTCAATTTTAACCTATTTTTTCTTATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCTACGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCa                                                              <  1:100846/251‑1 (MQ=255)                                                                                                                                                                                                           gAAACTAAAAACTCTTTTTGTTGATTGAGATACCCGATGCGTAAGCCAAAGTCCAGGTGCGGTGAACACAATGGCTAAACATTGCACCTTTCTTTTCCACCCCTCAATTTTAACCTATTTTTTCTTATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCTACGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGcaca                  >  2:18705/1‑250 (MQ=255)                                                                                                                                                                                                                 aaaaaCTCTTTTTGTTGATTGAGATACCCGATGCGTAAGCCAAAGTCCAGGTGCGGTGAACACAATGGCTAAACATTGCACCTTTCTTTTCCACCCCTCAATTTTAACCTATTTTTTCTTATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCTACGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTaaaaa                                          >  1:324271/1‑220 (MQ=255)                                                                                                                                                                                                                 aaaaaCTCTTTTTGTTGATTGAGATACCCGATGCGTAAGCCAAAGTCCAGGTGCGGTGAACACAATGGCTAAACATTGCACCTTTCTTTTCCACCCCTCAATTTTAACCTATTTTTTCTTATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCTACGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTaaaaa                                          <  2:324271/220‑1 (MQ=255)                                                                                                                                                                                                                           tttGTTGATTGAGATACCCGATGCGTAAGCCAAAGTCCAGGTGCGGTGAACACAATGGCTAAACATTGCACCTTTCTTTTCCACCCCTCAATTTTAACCTATTTTTTCTTATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCTACGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCCCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCca  <  2:22960/250‑1 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                            GAGCTACCCAAACAGAAGAAAAACTCAAGCGGGATGGCGTGAATTCAAAACAGCAAGCAAACACCACACACTTTGACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAGCAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTC‑‑CCCCTCAGTTTTAACCTATTTTTTCTTATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCTCTGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCCCTCCCCTTAAAAAATCTCAATCGTGACAATGCGCACAAATCGCTACCCTGCCA  >  NC_000913/3818385‑3818850

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