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

TTATTTGTTGACCCTGAATACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCT  >  NC_000913/252949‑253210                                                                                                                                        |                                                                                                                               ttATTTGTTGACCCTGAATACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAg                                                                                                                             <  2:289448/139‑1 (MQ=255)      tgttgACCCTGAATACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAag                                                                                                                        >  1:603846/1‑139 (MQ=255)           aCCCTGAATACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAagcagca                                                                                                                   >  2:713284/1‑139 (MQ=255)           aCCCTGAATACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAagcagca                                                                                                                   >  1:255757/1‑139 (MQ=255)           aCCCTGAATACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGCTATGGTTTTCAGGCAGTTAagcagca                                                                                                                   >  2:333017/1‑139 (MQ=255)                  aTACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGtt                                                                                                            <  2:603846/139‑1 (MQ=255)                      acCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAAt                                                                                                        >  1:980232/1‑139 (MQ=255)                      acCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAAt                                                                                                        >  2:629004/1‑139 (MQ=255)                           ccgTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCgg                                                                                                   >  2:638140/1‑139 (MQ=255)                           ccgTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTATTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCgg                                                                                                   >  1:947140/1‑139 (MQ=255)                            cgTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCggg                                                                                                  <  2:255757/139‑1 (MQ=255)                                        aGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTa                                                                                      >  1:177345/1‑139 (MQ=255)                                             tttGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGcc                                                                                                     <  2:966024/119‑1 (MQ=255)                                             tttGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGcc                                                                                                     >  1:966024/1‑119 (MQ=255)                                                   aaaaCCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTtat                                                                           <  1:638140/139‑1 (MQ=255)                                                      aCCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATg                                                                        <  1:629004/139‑1 (MQ=255)                                                                  gTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCg                                                            <  2:980232/139‑1 (MQ=255)                                                                     tGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATg                                                                        <  2:112621/124‑1 (MQ=255)                                                                     tGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATg                                                                        >  1:112621/1‑124 (MQ=255)                                                                       aaTCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAagca                                                                                                                      >  1:205617/1‑76 (MQ=255)                                                                       aaTCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAagca                                                                                                                      <  2:205617/76‑1 (MQ=255)                                                                                 tACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAg                                             >  2:65073/1‑139 (MQ=255)                                                                                        gACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCaa                                      >  1:317426/1‑139 (MQ=255)                                                                                          cGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAgcg                                                                                                                <  2:387169/63‑1 (MQ=255)                                                                                          cGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAgcg                                                                                                                >  1:387169/1‑63 (MQ=255)                                                                                                   aaCCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAAcg                           <  2:633371/139‑1 (MQ=255)                                                                                                     ccGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAAcgcc                         <  2:317426/139‑1 (MQ=255)                                                                                                     ccGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAAcgcc                         <  2:630139/139‑1 (MQ=255)                                                                                                     ccGAAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGAGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAAcgcc                         <  2:323649/139‑1 (MQ=255)                                                                                                        cAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAagcagca                                                                                                                   <  2:180810/46‑1 (MQ=255)                                                                                                        cAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAagcagca                                                                                                                   >  1:180810/1‑46 (MQ=255)                                                                                                         aaaaCCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATCAACCGCAACATTAAATCCAGCTTGCAATGAAAATAAcgcccgcc                     <  1:713284/139‑1 (MQ=255)                                                                                                         aaaaCCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAAcgcccgcc                     >  1:633838/1‑139 (MQ=255)                                                                                                          aaaCCCTTTTTTGAACGTTATGGTTTTCCGGCAGTTAAGCAGCAGCGCGTTGAATGcc                                                                                                     >  2:710798/1‑58 (MQ=39)                                                                                                          aaaCCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGcc                                                                                                     <  1:710798/58‑1 (MQ=255)                                                                                                             cccTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGtt                                                                                                            <  2:895878/48‑1 (MQ=255)                                                                                                             cccTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGtt                                                                                                            >  1:895878/1‑48 (MQ=255)                                                                                                                ttttttGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTAtg              >  1:885265/1‑139 (MQ=255)                                                                                                                    ttGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTgcg          <  2:633838/139‑1 (MQ=255)                                                                                                                            tATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCt  >  1:456694/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                               TTATTTGTTGACCCTGAATACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCT  >  NC_000913/252949‑253210

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