Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	253,084	A→G	T126A (ACA→GCA)	yafP →	GNAT family putative N‑acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	253,084	0	A	G	100.0%	141.3 / NA	39	T126A (ACA→GCA)	yafP	GNAT family putative N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/17);  total (22/17)

ACCCTGAATACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTG  >  NC_000913/252959‑253211                                                                                                                              |                                                                                                                                aCCCTGAATACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAagcagca                                                                                                                    >  1:631279/1‑139 (MQ=255)          acacCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGa                                                                                                           <  2:631279/139‑1 (MQ=255)            acCCGCCGTGGGGTTGCCAGCGCTGTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAAt                                                                                                         >  2:35159/1‑139 (MQ=255)                 ccgTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCgg                                                                                                    >  1:1098910/1‑139 (MQ=255)                       ggTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCg                                                                                              <  1:657394/139‑1 (MQ=255)                         ttGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTg                                                                                            <  1:514763/139‑1 (MQ=255)                            ccAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGtt                                                                                         >  2:194095/1‑139 (MQ=255)                                                 tGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCCGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGata                                                                    >  2:673517/1‑139 (MQ=255)                                                 tGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGata                                                                    >  1:1067120/1‑139 (MQ=255)                                                       aGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAAcc                                                              >  1:979447/1‑139 (MQ=255)                                                             aaTCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGtt                                                                                         >  2:623989/1‑106 (MQ=255)                                                             aaTCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGtt                                                                                         <  1:623989/106‑1 (MQ=255)                                                                  gAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAagcagca                                                                                                                    >  1:929276/1‑74 (MQ=255)                                                                  gAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAagcagca                                                                                                                    <  2:929276/74‑1 (MQ=255)                                                                  gAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTaaa                                                   >  2:1110277/1‑139 (MQ=255)                                                                   aaCTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATg                                                                         >  1:34060/1‑116 (MQ=255)                                                                   aaCTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATg                                                                         <  2:34060/116‑1 (MQ=255)                                                                   aaCTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAAt                                                  >  2:235033/1‑139 (MQ=255)                                                                            tGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGc                                         <  1:235033/139‑1 (MQ=255)                                                                              gACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCaa                                       >  1:42622/1‑139 (MQ=255)                                                                                cGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCg                                                             <  2:1011825/115‑1 (MQ=255)                                                                                cGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCg                                                             >  1:1011825/1‑115 (MQ=255)                                                                                      cATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCg                                                             <  2:731605/109‑1 (MQ=255)                                                                                      cATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCg                                                             >  1:731605/1‑109 (MQ=255)                                                                                           ccGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATcc                                                >  2:124505/1‑117 (MQ=255)                                                                                           ccGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATcc                                                <  1:124505/117‑1 (MQ=255)                                                                                               aaaaCCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAAcgcccgcc                      <  1:673517/139‑1 (MQ=255)                                                                                                 aaCCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTAc                                                                                      >  2:335070/1‑73 (MQ=255)                                                                                                 aaCCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTAc                                                                                      <  1:335070/73‑1 (MQ=255)                                                                                                   cccTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGt                  >  1:75281/1‑139 (MQ=255)                                                                                                               cGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAAttttt                                                                              >  1:352733/1‑67 (MQ=255)                                                                                                               cGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAAttttt                                                                              <  2:352733/67‑1 (MQ=255)                                                                                                                 ttATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTc    <  2:1067120/139‑1 (MQ=255)                                                                                                                  tATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCt   >  2:494785/1‑139 (MQ=255)                                                                                                                   aTGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTg  <  1:194095/139‑1 (MQ=255)                                                                                                                    tGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACa                                                        <  2:949198/84‑1 (MQ=255)                                                                                                                    tGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACa                                                        >  1:949198/1‑84 (MQ=255)                                                                                                                          tCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCa                                                           >  2:1085905/1‑75 (MQ=255)                                                                                                                          tCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCa                                                           <  1:1085905/75‑1 (MQ=255)                                                                                                                              |                                                                                                                                ACCCTGAATACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTG  >  NC_000913/252959‑253211

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