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%	97.9 / NA	28	intergenic (+224/‑13)	yhhY/yhhZ	N‑acetyltransferase YhhY/putative endonuclease YhhZ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (11/17);  total (11/17)

CGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAGGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTAT  >  NC_000913/3581613‑3582079                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                      agcgGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATGCAAAATCGACATTCCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACTCCTACAATTACCTGTAGTAAATAAATAACAATTAACAAGCATAATAGTGCCTTTGAGAAAATAGCATGCCATAAAAAGACTTTTCAGGGa                                                                                                                                                                                                                          <  2:315496/250‑1 (MQ=255)      gTGAAGAAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGt                                                                                                                                                                                                                     >  2:148893/1‑251 (MQ=255)               gATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAATTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAGCAATTAACAAGCATAATATTGCCATTGATAAAATGGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGtt                                                                                                                                                                                                            <  1:260664/251‑1 (MQ=255)                     gCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCt                                                                                                                                                                                                       <  2:356288/250‑1 (MQ=255)                     gCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTg                                                                                                                                                                                                      <  1:414897/251‑1 (MQ=255)                                gagaTGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACgg                                                                                                                                                                                            <  2:360121/250‑1 (MQ=255)                                   aTGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGgaga                                                                                                                                                                                        <  2:129027/251‑1 (MQ=255)                                                       aaaaTCGACATACCTTCCTTTAATGTATTTATTTGCCCAATACATATATTGAGATTTAACGCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGATTAAATAAATAACAATTAACAAGCATAATATTGCCATGTATAAAATAGCAGGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGTAGAACAACAAGGAAGCATCTCCg                                                                                                                                                                    <  2:417986/251‑1 (MQ=255)                                                                            aaGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGt                                                                                                                                               <  1:237763/251‑1 (MQ=255)                                                                              ggTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTcc                                                                                                                                                                     >  2:98003/1‑227 (MQ=255)                                                                              ggTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTcc                                                                                                                                                                     <  1:98003/227‑1 (MQ=255)                                                                                             aaTACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGAGAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGTAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTgg                                                                                                                              <  1:279644/251‑1 (MQ=255)                                                                                               tACATATATTGAGATTTAACTCATAATTACTTCCTCAAGTGTAATATTTTATTTTTTAATATATACGCTTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGTTGCGTAATATTGTTTACCTGACAGTAACGGGCGAACAACAAGGAAGCATCTCCGAAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGc                                                                                                                             >  2:49794/1‑250 (MQ=255)                                                                                                                                                  ttttttAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGATGAAATATTTACATTCTCACTCTTAAATAATATta                                                                          >  1:260331/1‑250 (MQ=255)                                                                                                                                                     tttAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGATGAAATATTTACATTCTCACTCTTAAATAATATtaata                                                                       >  1:330304/1‑250 (MQ=255)                                                                                                                                                        aaTATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATAc                                                                   <  2:79877/251‑1 (MQ=255)                                                                                                                                                             ataCGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGc                                                              >  1:260891/1‑251 (MQ=255)                                                                                                                                                                                 aGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGAt                                                                                                            >  1:282523/1‑185 (MQ=255)                                                                                                                                                                                 aGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGAt                                                                                                            <  2:282523/185‑1 (MQ=255)                                                                                                                                                                                                   aaCAAGCATAATATTGCCATTGATACAATAGAATGCCATCAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGAATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATAGTTGGCAGAGCGGGCATGAgg                                                                                                              >  2:60440/1‑165 (MQ=255)                                                                                                                                                                                                   aaCAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAAGCGTTGGCAGAGCGGGCATGAgg                                                                                                              <  1:60440/165‑1 (MQ=255)                                                                                                                                                                                                       aGCATAAAATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAAGAAGCATCTCCGCAGGTTGTGGGGCTTCTGAGTCTACAGGTAATCGTTGTCAGAGCGGACATGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGTTGCTTGGTTCACATTTCCATGGTATAACTTTTTCTTTATTAATTg                    <  2:241869/251‑1 (MQ=255)                                                                                                                                                                                                         cATAATATTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGGTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGTGAATCGTTGGCAGAGCGGGCATGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAAGACTGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGAt                  <  1:282276/251‑1 (MQ=255)                                                                                                                                                                                                            aataTTGCCATTGATAAAATAGCATGCCATAAAAAGACTTTTCAAGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATaaa               <  1:250520/251‑1 (MQ=255)                                                                                                                                                                                                                     aTTGATAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACACGGAAGCATCTCCGCAGGTTGTGGGACTTCTGCGTCTACAGGTAATCGTTGGCAGAGCGGGCATGCGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTcc       >  1:35831/1‑250 (MQ=255)                                                                                                                                                                                                                         aTAAAATAGCATGCCATAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCattat  >  1:158610/1‑251 (MQ=255)                                                                                                                                                                                                                                        aTAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGCTGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGtt                                                 >  2:398294/1‑189 (MQ=255)                                                                                                                                                                                                                                        aTAAAAAGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGtt                                                 <  1:398294/189‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                      CGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAGGACTTTTCAGGGATGAGTAATATTGTTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTAT  >  NC_000913/3581613‑3582079

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