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%	173.8 / NA	50	T126A (ACA→GCA)	yafP	GNAT family putative N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/24);  total (26/24)

ACCCTGAATACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTT  >  NC_000913/252959‑253206                                                                                                                              |                                                                                                                           aCCCTGAATACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAagcagca                                                                                                               >  2:873803/1‑139 (MQ=255)            acCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAAt                                                                                                    >  2:71200/1‑139 (MQ=255)                gccgTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCg                                                                                                <  1:1179066/139‑1 (MQ=255)                 ccgTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAATACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAagcag                                                                                                                 <  2:1229328/121‑1 (MQ=255)                 ccgTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAagcag                                                                                                                 >  1:1229328/1‑121 (MQ=255)                 ccgTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCgg                                                                                               >  1:658652/1‑139 (MQ=255)                       ggTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCg                                                                                         <  2:658652/139‑1 (MQ=255)                       ggTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCg                                                                                         <  2:887557/139‑1 (MQ=255)                            ccAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGtt                                                                                                        <  1:997469/119‑1 (MQ=255)                            ccAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGtt                                                                                                        >  2:997469/1‑119 (MQ=255)                                           atccTGTGATTAAGTCTGAATCCCATCTTACTGTTGACGTAATCATAACCGCAAATCTCTTTTTTGAACGTTATGGTTGTCAGGCATTTATGCAGTAGCGCGTAGAATGCCGGGGATCGTGGCTTACTAATTTTTatat                                                                     <  1:7392/137‑1 (MQ=255)                                               tttGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGa                                                                 >  2:716396/1‑139 (MQ=255)                                               tttGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGa                                                                 >  2:2601/1‑139 (MQ=255)                                                 tGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTatat                                                                     >  2:597342/1‑133 (MQ=255)                                                 tGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTatat                                                                     <  1:597342/133‑1 (MQ=255)                                                       aGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGtt                                                                                                        <  1:1099494/92‑1 (MQ=255)                                                       aGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGtt                                                                                                        >  2:1099494/1‑92 (MQ=255)                                                           tGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCg                                                                                         >  1:1009961/1‑103 (MQ=255)                                                           tGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCg                                                                                         <  2:1009961/103‑1 (MQ=255)                                                                         cGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGat                                                                <  2:964628/114‑1 (MQ=255)                                                                         cGGTGGACGCAAGCATAACCGCAAAACACTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGat                                                                >  1:964628/1‑114 (MQ=255)                                                                               aCGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTaa                                                                              >  1:1361865/1‑94 (MQ=255)                                                                               aCGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTaa                                                                              <  2:1361865/94‑1 (MQ=255)                                                                               aCGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAAt                                 >  2:1545402/1‑139 (MQ=255)                                                                                cGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATg                                >  1:571838/1‑139 (MQ=255)                                                                                     gCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAAt                           <  1:873803/139‑1 (MQ=255)                                                                                      cATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAAttttt                                                                         >  2:1588308/1‑92 (MQ=255)                                                                                      cATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAAttttt                                                                         <  1:1588308/92‑1 (MQ=255)                                                                                        tAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAAc                        >  1:1580980/1‑139 (MQ=255)                                                                                        tAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAAc                        >  2:54697/1‑139 (MQ=255)                                                                                         aaCCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCg                                                                                         <  2:1211384/73‑1 (MQ=255)                                                                                         aaCCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCg                                                                                         >  1:1211384/1‑73 (MQ=255)                                                                                         aaCCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAAcg                       >  1:486659/1‑139 (MQ=255)                                                                                           ccGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAAcgcc                     <  1:1545402/139‑1 (MQ=255)                                                                                               aaaaCCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAagca                                                                                                                  >  1:813728/1‑42 (MQ=255)                                                                                               aaaaCCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAagca                                                                                                                  <  2:813728/42‑1 (MQ=255)                                                                                               aaaaCCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAAcgcccgcc                 <  1:54697/139‑1 (MQ=255)                                                                                               aaaaCCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAAcg                       <  1:504682/133‑1 (MQ=255)                                                                                               aaaaCCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAAcg                       >  2:504682/1‑133 (MQ=255)                                                                                                   cccTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATg                                                                    >  2:73735/1‑84 (MQ=255)                                                                                                   cccTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATg                                                                    <  1:73735/84‑1 (MQ=255)                                                                                                             aaCGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGaa                                                                                                     >  2:869433/1‑41 (MQ=255)                                                                                                             aaCGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGaa                                                                                                     <  1:869433/41‑1 (MQ=255)                                                                                                              aCGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAAttttt                                                                         <  1:189060/68‑1 (MQ=255)                                                                                                              aCGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAAttttt                                                                         >  2:189060/1‑68 (MQ=255)                                                                                                              aCGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATg                                                                    <  1:1555526/73‑1 (MQ=255)                                                                                                              aCGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATg                                                                    >  2:1555526/1‑73 (MQ=255)                                                                                                              aCGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGtt  <  2:486659/139‑1 (MQ=255)                                                                                                                          tCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATg                                                                    <  2:887667/61‑1 (MQ=255)                                                                                                                          tCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATg                                                                    >  1:887667/1‑61 (MQ=255)                                                                                                                              |                                                                                                                           ACCCTGAATACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTT  >  NC_000913/252959‑253206

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