Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,818,828	G→A	intergenic (+244/‑46)	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,828	0	G	A	100.0%	94.7 / NA	32	intergenic (+244/‑46)	dinD/yicG	DNA damage‑inducible protein D/PF03458 family inner membrane protein YicG
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/11);  total (21/11)

GAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTC‑‑CCCCTCAGTTTTAACCTATTTTTTCTTATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCT‑CTGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCCCTCCCCTTAAAAAATCTCAATCGTGACAATGCGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATT  >  NC_000913/3818587‑3819071                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    gAAACTAAAAACTCTTTTTGTTGATTGAGATACCCGATGCGTAAGCCAAAGTCCAGGTGCGGTGAACACAATGGCTAAACATTGCACCTTTCTTTTCCACCCCTCAATTTTAACCTATTTTTTCTTATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCT‑ACGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGcaca                                                                                                                                                                                                                                               >  2:18705/1‑250 (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)                                                                                                                                                                                                                        cctcctTAAAAAATCTCAACCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTCCTGCCGAAGCCATGACGGGGGCGCTGGCGGCCGGGCGTCGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTGGGt                        <  2:126550/248‑1 (MQ=255)                                                                                                                                                                                                                         cacctTAAAAAATCTCAAGCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTgg                          <  1:219269/246‑1 (MQ=255)                                                                                                                                                                                                                          accttccaaaaTCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCCCCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTGGGTc                       >  2:182965/8‑250 (MQ=255)                                                                                                                                                                                                                               aaaaaaTCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTGGGTCAAacac                 <  1:232032/251‑1 (MQ=255)                                                                                                                                                                                                                                            cGTGACAATACGCACAAATCGCTACCCTGCCAGACCGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTGGGTCAAACACCCGGAATAtgtgc    >  1:102221/1‑250 (MQ=255)                                                                                                                                                                                                                                               gACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTGGGTCAAACACCCGGAATATGTGAtt  <  1:18705/250‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    GAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTC‑‑CCCCTCAGTTTTAACCTATTTTTTCTTATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCT‑CTGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCCCTCCCCTTAAAAAATCTCAATCGTGACAATGCGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATT  >  NC_000913/3818587‑3819071

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