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

ATTTGTTGACCCTGAATACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAA‑CCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTC  >  NC_000913/252951‑253215                                                                                                                                      |                                                                                                                                     aTTTGTTGACCCTGAATACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGtt                                                                                                                                 >  2:948474/1‑139 (MQ=255)     gttgACCCTGAATACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAagc                                                                                                                             >  2:647929/1‑139 (MQ=255)        gACCCTGAATACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAagcagc                                                                                                                          <  2:498526/139‑1 (MQ=255)         aCCCTGAATACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAagcagca                                                                                                                         <  1:647929/139‑1 (MQ=255)           ccTGAATACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAagcagcagc                                                                                                                       >  2:340590/1‑139 (MQ=255)                aTACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGtt                                                                                                                  <  1:340590/139‑1 (MQ=255)                    acCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAg                                                                                                                                   >  2:225725/1‑118 (MQ=255)                    acCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAg                                                                                                                                   <  1:225725/118‑1 (MQ=255)                    acCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAAt                                                                                                              >  1:319625/1‑139 (MQ=255)                         ccgTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCgg                                                                                                         >  1:828372/1‑139 (MQ=255)                         ccgTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCgg                                                                                                         >  1:1287715/1‑139 (MQ=255)                         ccgTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCgg                                                                                                         >  2:264377/1‑139 (MQ=255)                          cgTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCggg                                                                                                        <  2:319625/139‑1 (MQ=255)                                      aGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTa                                                                                            <  2:1287715/139‑1 (MQ=255)                                        cgcTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACt                                                                                          <  1:1158960/139‑1 (MQ=255)                                                      cTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCg                                                                            <  1:264377/139‑1 (MQ=255)                                                      cTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCg                                                                            <  1:1137998/139‑1 (MQ=255)                                                         tGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATg                                                                              >  1:1303767/1‑134 (MQ=255)                                                         tGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATg                                                                              <  2:1303767/134‑1 (MQ=255)                                                         tGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGG‑TTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTatat                                                                               <  1:97039/132‑1 (MQ=255)                                                         tGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGG‑TTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTatat                                                                               >  2:97039/1‑132 (MQ=255)                                                                 tCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAA‑CCGc                                                                <  1:1076076/139‑1 (MQ=255)                                                                           aaCTTTCGGTGGACGCAAGCATAACCGCAATACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAA‑CCGAAACATTAAAt                                                      >  2:587402/1‑139 (MQ=255)                                                                           aaCTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAA‑CCg                                                                 <  1:1113488/128‑1 (MQ=255)                                                                           aaCTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAA‑CCg                                                                 >  2:1113488/1‑128 (MQ=255)                                                                           aaCTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAA‑CCGCAACATTAAAt                                                      >  2:128997/1‑139 (MQ=255)                                                                                      gACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAg                                                                                                     <  1:800435/82‑1 (MQ=255)                                                                                      gACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAg                                                                                                     >  2:800435/1‑82 (MQ=255)                                                                                        cGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAA‑CCGCAACATTAAATCCAGCTTGCAATg                                         >  1:463612/1‑139 (MQ=255)                                                                                           aaGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCg                                                                                                   >  1:398847/1‑79 (MQ=255)                                                                                           aaGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCg                                                                                                   <  2:398847/79‑1 (MQ=255)                                                                                              cATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAA‑CCGCAACATTAAATCCAGCTTGCAATGAAAATa                                   >  2:592848/1‑139 (MQ=255)                                                                                                 aaCCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAA‑CCGCAACATTAAATCCAGCTTGCAATGAAAATAAcg                                <  2:463612/139‑1 (MQ=255)                                                                                                 aaCCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAA‑CCGCAACATTAAATCCAGCTTGCAATGAAAATAAcg                                <  1:128997/139‑1 (MQ=255)                                                                                                   ccGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAA‑CCGCa                                                               <  2:558875/106‑1 (MQ=255)                                                                                                   ccGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAA‑CCGCa                                                               >  1:558875/1‑106 (MQ=255)                                                                                                   ccGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAA‑CCGCAACATTAAATCCAGCTTGCAATGAAAATAAcgcc                              <  1:592848/139‑1 (MQ=255)                                                                                                    cGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCCGCAACATTAAATCCAGCTTGCAATGAAAATAAcgcc                              <  1:1014072/139‑1 (MQ=255)                                                                                                       aaaaCCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGtt                                                                                                                  >  2:117890/1‑52 (MQ=255)                                                                                                       aaaaCCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGtt                                                                                                                  <  1:117890/52‑1 (MQ=255)                                                                                                       aaaaCCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAA‑CCGCAACATTAAATCCAGCTTGCAATGAAAATAAcgcccgcc                          >  2:632107/1‑139 (MQ=255)                                                                                                       aaaaCCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAA‑CCGCAACATTAAATCCAGCTTGCAATGAAAATAAcgcccgcc                          >  2:993199/1‑139 (MQ=255)                                                                                                       aaaaCCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAA‑CCGCAACATTAAATCCAGCTTGCAATGAAAATAAcgcccgcc                          <  1:632107/139‑1 (MQ=255)                                                                                                         aaCCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGTCGGGGAGCg                                                                                                   >  2:1331733/1‑65 (MQ=255)                                                                                                         aaCCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGTCGGGGAGCg                                                                                                   <  1:1331733/65‑1 (MQ=255)                                                                                                            ccTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAA‑CCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTa                     <  1:993199/139‑1 (MQ=255)                                                                                                                   tGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAA‑CCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTgcgc              >  2:148278/1‑139 (MQ=255)                                                                                                                   tGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAA‑CCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTgcgc              >  2:848284/1‑139 (MQ=255)                                                                                                                    gAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAA‑CCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGcgcg             <  1:848284/139‑1 (MQ=255)                                                                                                                       cGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAA‑CCGCAACATTAAATCCAGCTTGCAATGAAAATAAcgcc                              >  2:381493/1‑119 (MQ=255)                                                                                                                       cGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAA‑CCGCAACATTAAATCCAGCTTGCAATGAAAATAAcgcc                              <  1:381493/119‑1 (MQ=255)                                                                                                                          tATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAA‑CCg                                                                 <  2:1235838/81‑1 (MQ=255)                                                                                                                          tATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAA‑CCg                                                                 >  1:1235838/1‑81 (MQ=255)                                                                                                                          tATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAA‑CCGCAACATTAAATCCAGCt                                                <  1:1077831/98‑1 (MQ=255)                                                                                                                          tATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAA‑CCGCAACATTAAATCCAGCt                                                >  2:1077831/1‑98 (MQ=255)                                                                                                                          tATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAA‑CCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCt       >  2:310411/1‑139 (MQ=255)                                                                                                                          tATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAA‑CCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCt       >  1:305408/1‑139 (MQ=255)                                                                                                                           aTGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGaa                                                                                                               <  1:721228/35‑1 (MQ=255)                                                                                                                           aTGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGaa                                                                                                               >  2:721228/1‑35 (MQ=255)                                                                                                                               ttttCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAA‑CCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTc  >  2:272673/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                     ATTTGTTGACCCTGAATACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAA‑CCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTC  >  NC_000913/252951‑253215

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