Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,581,850	G→A	intergenic (+224/‑13)	yhhY → / → yhhZ	N‑acetyltransferase YhhY/putative endonuclease YhhZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,581,850	0	G	A	100.0%	107.2 / NA	32	intergenic (+224/‑13)	yhhY/yhhZ	N‑acetyltransferase YhhY/putative endonuclease YhhZ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/15);  total (17/15)

TATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTA‑CTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAGGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAA  >  NC_000913/3581603‑3582086                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                             tattatATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTA‑CTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGAc                                                                                                                                                                                                                                           >  1:232923/1‑251 (MQ=255)       atGGCGCGGGAGAAGTAAGATAGGGCCCTTTTTCTGAGCTGGAAAGAGGGGGTCATTCAAAATCGACATACCTTCCTTTAAGGAATTTAGTTGCCCAATACATATATTGAGATTTAATTAATAATTA‑CTTCTTACAGAGTAATATTTTATTTTTTAATATATACATCTACAATTTCCTGGAGTAATTTAATATTAATTAACAAGCATAATATTGTACTTGATAAAATAGCATGACATAAAAAAACtttta                                                                                                                                                                                                                                      >  2:17807/1‑249 (MQ=255)         ggCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTA‑CTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAg                                                                                                                                                                                                                                    <  2:13217/250‑1 (MQ=255)            gcgGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTA‑CTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATACTATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGAt                                                                                                                                                                                                                                >  2:427706/1‑251 (MQ=255)               ggTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTA‑CTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGGCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAg                                                                                                                                                                                                                             <  1:389994/251‑1 (MQ=255)                                         tGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTA‑CTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACgg                                                                                                                                                                                                   >  1:275301/1‑251 (MQ=255)                                          gagaTGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTA‑CTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACggg                                                                                                                                                                                                  >  1:436873/1‑251 (MQ=255)                                                   aaaGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTA‑CTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGaacaac                                                                                                                                                                                          >  2:53515/1‑250 (MQ=255)                                                                cAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTA‑CTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTcc                                                                                                                                                                            >  1:354149/1‑251 (MQ=255)                                                                           tcctttccttAAAGGT‑TTTATTTGGCCAATACAATTAATGAGAATTAACTCCAAATTAATTTCTTAAATTGTA‑TATTTTAGTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTtgtg                                                                                                                                                                  <  1:158981/245‑1 (MQ=255)                                                                                 ccttTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTA‑CTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTc                                                                                                                                                           >  1:450061/1‑251 (MQ=255)                                                                                  cttTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTA‑CTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCt                                                                                                                                                          >  2:499967/1‑251 (MQ=255)                                                                                                          aCATATATTGAGATTTAGCTCATAATTA‑CTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCTCAGATTGGGGGACTTCTGATTCTACAGGTAAGCGTTGGCag                                                                                                                                  <  2:81412/251‑1 (MQ=255)                                                                                                                atTGAGATTTAACTCATAATTA‑CTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGCGTCTACAGGTAATCGTTGGCAGAGCggg                                                                                                                            >  1:378604/1‑251 (MQ=255)                                                                                                                           aCTCATAATTA‑CTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGATg                                                                                                                  >  1:511133/1‑250 (MQ=255)                                                                                                                            cTCATAATAA‑CTTCCTTAAGTGTAATCTTTTATTTTTTGATATATACGCCTACAATTTCCTGGAGTAAATCACTAACAATTATCAAGCATAATATTGTCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGCGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAgcctga                                                                                                                 >  2:265892/1‑245 (MQ=255)                                                                                                                               aTAATTA‑CTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAGTCGTTGGCAGAGCGGGCATGAGGATGAAata                                                                                                             <  1:558477/251‑1 (MQ=255)                                                                                                                                             aaaGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGATGAAATATTTACATTCTCAc                                                                                                <  1:122544/251‑1 (MQ=255)                                                                                                                                                          ttattttTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGATGAAATATTTACATTCTCACTCTTAAATAatat                                                                                   >  1:175017/1‑251 (MQ=255)                                                                                                                                                               ttttAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCAATGATAAAATAGCATGCCATAAAAATACTTTTCAGGTATGAGTAATAGTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGTTAATCGTTGGCAGAGCGGGCATGAGGATGGAATATTGACATTCTCTCTCTTAAATAATATtaat                                                                               <  2:364334/250‑1 (MQ=255)                                                                                                                                                                           cGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTACTAATACGGGGCTTg                                                                  >  2:493681/1‑251 (MQ=255)                                                                                                                                                                              cTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGtt                                                               >  2:79550/1‑251 (MQ=255)                                                                                                                                                                                  aaTTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATAGTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACTGGGCTTGGTTcaca                                                           <  2:232923/251‑1 (MQ=255)                                                                                                                                                                                      tCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCAATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGtt                                                        <  1:17807/250‑1 (MQ=255)                                                                                                                                                                                              taaataaaTAACAATTAACAAGCCTAATATTGCCATTGCTAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATACTACGGGGCTTGGTTCACAGTTCCATGGtat                                               >  2:365865/1‑251 (MQ=255)                                                                                                                                                                                                 ataaataaCAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAAc                                            >  1:76606/1‑251 (MQ=255)                                                                                                                                                                                                                             gCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACt                <  1:437683/251‑1 (MQ=255)                                                                                                                                                                                                                             gCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACt                <  1:499967/251‑1 (MQ=255)                                                                                                                                                                                                                             gCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTGCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACt                <  2:275301/251‑1 (MQ=255)                                                                                                                                                                                                                                aTTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCa             <  2:436873/251‑1 (MQ=255)                                                                                                                                                                                                                                   gATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCatt           <  2:470543/250‑1 (MQ=255)                                                                                                                                                                                                                                           aGCATGCCATAAAAAGACTTTTCAGTGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGTACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTaa  <  1:313540/251‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                             TATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTA‑CTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAGGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAA  >  NC_000913/3581603‑3582086

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