Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,818,801	C→A	intergenic (+217/‑73)	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,801	0	C	A	100.0%	82.7 / NA	27	intergenic (+217/‑73)	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 A (20/7);  total (20/7)

GAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTC‑‑CCCCTCAGTTTTAACCTATTTTTTCTTATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCT‑CTGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCCCTCCCCTTAAAAAATCTCAATCGTGACAATGCGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTGG  >  NC_000913/3818587‑3819047                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                       gAAACTAAAAACTCTTTTTGTTGATTGAGATACCCGATGCGTAAGCCAAAGTCCAGGTGCGGTGAACACAATGGCTAAACATTGCACCTTTCTTTTCCACCCCTCAATTTTAACCTATTTTTTCTTATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCT‑ACGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGcaca                                                                                                                                                                                                                       >  2:18705/1‑250 (MQ=255)       aaaaaCTCTTTTTGTTGATTGAGATACCCGATGCGTAAGCCAAAGTCCAGGTGCGGTGAACACAATGGCTAAACATTGCACCTTTCTTTTCCACCCCTCAATTTTAACCTATTTTTTCTTATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCT‑ACGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTaaaaa                                                                                                                                                                                                                                               <  2:324271/220‑1 (MQ=255)       aaaaaCTCTTTTTGTTGATTGAGATACCCGATGCGTAAGCCAAAGTCCAGGTGCGGTGAACACAATGGCTAAACATTGCACCTTTCTTTTCCACCCCTCAATTTTAACCTATTTTTTCTTATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCT‑ACGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTaaaaa                                                                                                                                                                                                                                               >  1:324271/1‑220 (MQ=255)                 tttGTTGATTGAGATACCCGATGCGTAAGCCAAAGTCCAGGTGCGGTGAACACAATGGCTAAACATTGCACCTTTCTTTTCCACCCCTCAATTTTAACCTATTTTTTCTTATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCT‑ACGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCCCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCca                                                                                                                                                                                                       <  2:22960/250‑1 (MQ=255)                                   cGATGCGTAAGCCAAAGTCCAGGTGCGGTGAACACAATGGCTAAACATTGCACCTTTCTTTTCCACCCCTCAATTTTAACCTATTTTTTCTTATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCT‑ACGCTACAGAAAATAGCGATTTCACATAACTACATGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGgag                                                                                                                                                                                      >  2:52327/1‑249 (MQ=255)                                    gATGCGTAAGCCAAAGTCCAGGTGCGGTGAACACAATGGCTAAACATTGCACCTTTCTTTTCCACCCCTCAATTTTAACCTATTTTTTCTTATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCT‑ACGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGgag                                                                                                                                                                                      >  1:177315/1‑248 (MQ=255)                                                                           cTAAACATTGCACCTTTCTTTTCCACCCCTCAATTTTAACCTATTTTTTCTTATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCT‑ACGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACt                                                                                                                                            >  1:238330/1‑251 (MQ=255)                                                                                     cACCTTTCTTTTCCACCCCTCAATTTTAACCTATTTTTTCTTATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCT‑ACGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACACATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGcc                                                                                                                                   >  1:109437/1‑250 (MQ=255)                                                                                                             tttAACCTATTTTTTCTTATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCT‑ACGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCgggc                                                                                                          >  1:120856/1‑249 (MQ=255)                                                                                                              ttAACCTATTTTTTCTTATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCT‑ACGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCgggcg                                                                                                         >  2:379568/1‑248 (MQ=255)                                                                                                                aaCCTATTTTTTCTTATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCT‑ACGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAggg                                                                                                                                                                                        >  1:44181/1‑170 (MQ=255)                                                                                                                aaCCTATTTTTTCTTATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCT‑ACGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGAGTTTTAggg                                                                                                                                                                                        <  2:44181/170‑1 (MQ=255)                                                                                                                        ttttCTTATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCT‑ACGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATg                                                                                               >  2:232032/1‑251 (MQ=255)                                                                                                                           tCTTATGCATTTTCTCAGAAAGGAAGACCGAAAAATGCATCTACT‑CTACCGAACATAGCGATTTCCCCTAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACACATCGCTACCCTGCCAGACAGATTTTTAGGTAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGTCGGCCGGTCGTCGCCGCATGGAt                                                                                            >  2:424412/1‑251 (MQ=255)                                                                                                                             ttATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCT‑ACGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATAc                                                                                          >  2:219269/1‑251 (MQ=255)                                                                                                                                 gCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCT‑ACGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACAttt                                                                                      >  1:45289/1‑251 (MQ=255)                                                                                                                                                   agTCAGAAGAATGCATCT‑ACGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCg                                                                    >  2:126295/1‑251 (MQ=255)                                                                                                                                                       agaagaATGCATCT‑ACGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCg                                                                >  1:88481/1‑251 (MQ=255)                                                                                                                                                           gaaTGCATCT‑ACGCAACAGAAACTAGCGATTTCACATAACTACAAGATATCAATTTCCTCTCACCTTAAAAAATCTCAATCGAGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCAGGGCGACGCCGCATGGATACCTTTGGAGTGATTATTATTGCGACAGCCAc                                                            >  2:383755/1‑251 (MQ=255)                                                                                                                                                                           aCAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAg                                              <  2:45289/250‑1 (MQ=255)                                                                                                                                                                               aaaaTAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCa                                         <  1:379568/251‑1 (MQ=255)                                                                                                                                                                                           tCACATATCTACAAGTTATCAATTTCCTCTCACCTGAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTGGGCATTACTGCCGAAGCCAGGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCaa                                                       <  2:140608/225‑1 (MQ=255)                                                                                                                                                                                           tCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCaa                                                       >  1:140608/1‑225 (MQ=255)                                                                                                                                                                                                 aaCTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTctgct                        <  2:120856/250‑1 (MQ=255)                                                                                                                                                                                                            aTCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATcc            >  1:243762/1‑251 (MQ=255)                                                                                                                                                                                                                aTTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTc        >  1:436836/1‑251 (MQ=255)                                                                                                                                                                                                                  ttCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAATTCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGATGTTACACATTTTGTATATGGTTGGCATTACTGCCGAAGCCATGACGGGAGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGACACCGCAATTGTCGGCGGGTCGGTGCTGGTTATTCTTCTGGGCACCTATACGCTCgg      >  2:141833/1‑251 (MQ=255)                                                                                                                                                                                                                      cctcacctTAAAAAATCGCAACCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTTCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCTCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGAGATTCTGCTGGGCCACTATCCGCTCGGCTgg  <  2:415977/251‑1 (MQ=255)                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                       GAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTC‑‑CCCCTCAGTTTTAACCTATTTTTTCTTATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCT‑CTGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCCCTCCCCTTAAAAAATCTCAATCGTGACAATGCGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTGG  >  NC_000913/3818587‑3819047

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