Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,078,438	G→A	*330* (TAG→TAA)	yiiD →	GNAT family putative N‑acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,078,438	0	G	A	100.0%	147.0 / NA	46	*330* (TAG→TAA)	yiiD	GNAT family putative N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/30);  total (16/30)

GACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAGCTGATGGAGAGGGGCGGTGCTGCCT‑CTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACGTGTCTTGCTACATCA  >  NC_000913/4078302‑4078570                                                                                                                                         |                                                                                                                                      gACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACt                                                                                                                                     >  1:601537/1‑139 (MQ=255)   cGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGa                                                                                                                                   <  1:324000/139‑1 (MQ=255)            cGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAggg                                                                                                                          <  2:160695/139‑1 (MQ=255)                tGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCgg                                                                                                                      >  2:71019/1‑139 (MQ=255)                tGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCgg                                                                                                                      >  2:634189/1‑139 (MQ=255)                  cAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtg                                                                                                                    <  2:818106/139‑1 (MQ=255)                  cAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtg                                                                                                                    <  2:54776/139‑1 (MQ=255)                     aTGCAGGTCGAATTCTTTGGCGACGGGACGCCGGGTGCAGTGTTTGAAGGCCCTTTTATCTTTCTGCCCCCGAGGCCATTTGGCCCGTATGAAGGGGGCGGGAACGAAGAA‑AGGAACTGATGGAGAGGGGCGGtgctgc                                                                                                                <  2:446291/139‑1 (MQ=255)                      tGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtgctgc                                                                                                                <  1:829278/139‑1 (MQ=255)                          ggTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCct‑ct                                                                                                           <  2:896749/139‑1 (MQ=255)                          ggTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCct‑ct                                                                                                           >  2:380969/1‑139 (MQ=255)                                aaTCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCT‑CTCTCAtt                                                                                                     <  1:380969/139‑1 (MQ=255)                                    tttGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCT‑CTCTCATTCAgg                                                                                                 <  2:828201/139‑1 (MQ=255)                                    tttGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCT‑CTCTCATTCAgg                                                                                                 <  1:278653/139‑1 (MQ=255)                                    tttGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCT‑CTCTCATTCAgg                                                                                                 <  2:1043294/139‑1 (MQ=255)                                      tGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTCCTGCCCGCGAAGCCATTGGGCCCGTATGAAGAGGGCGGGAACGAAGAAGGGTGACTGATGTAGAGGGGCGGTGCTGACT‑CTCTCATTCAGGTc                                                                                               <  2:271054/139‑1 (MQ=255)                                        gCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCT‑CTCTCATTCAGGTCAt                                                                                             <  2:341879/139‑1 (MQ=255)                                        gCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCT‑CTCTCATTCAGGTCAt                                                                                             <  1:805480/139‑1 (MQ=255)                                        gCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCT‑CTCTCATTCAGGTCAt                                                                                             <  1:834107/139‑1 (MQ=255)                                        gCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCT‑CTCTCATTCAGGTCAt                                                                                             <  2:187998/139‑1 (MQ=255)                                        gCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCT‑CTCTCATTCAGGTCAt                                                                                             <  1:414027/139‑1 (MQ=255)                                         cgacgaGACGCCGGGTGCAGTGTTTGAATGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCT‑CTCTCATTCAGGTCAtt                                                                                            <  1:634189/139‑1 (MQ=255)                                             gagaCGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCT‑CTCTCATTCAGGTCATTtata                                                                                        <  1:581883/139‑1 (MQ=255)                                                             tgtTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCT‑CTCTCATTCAGGTCATTTATATTAGTCCTTGGCTAtt                                                                        <  2:601537/139‑1 (MQ=255)                                                                  gAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCT‑CTCTCATTCAGGTCATTTATATTAGTcc                                                                                 >  1:544920/1‑125 (MQ=255)                                                                  gAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCT‑CTCTCATTCAGGTCATTTATATTAGTcc                                                                                 <  2:544920/125‑1 (MQ=255)                                                                         cGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCT‑CTCTCATTCAGGTCAt                                                                                             >  1:417236/1‑106 (MQ=255)                                                                         cGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCT‑CTCTCATTCAGGTCAt                                                                                             <  2:417236/106‑1 (MQ=255)                                                                         cGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCT‑CTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTa                                                            <  2:794426/139‑1 (MQ=255)                                                                                    cTGCCCGCGAAGCCATTTGGCCCGTATGGAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCT‑CTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAActc                                                 >  2:322639/1‑139 (MQ=255)                                                                                    cTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCT‑CTCTCATTCAGGTCAt                                                                                             >  1:1012496/1‑95 (MQ=255)                                                                                    cTGCCCGCGAAGCCATTTGGCCCGGATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCT‑CTCTCATTCAGGTCAt                                                                                             <  2:1012496/95‑1 (MQ=255)                                                                                     tGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCT‑CTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAActct                                                >  2:1016705/1‑139 (MQ=255)                                                                                       cccGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCT‑CTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCtat                                              >  2:391581/1‑139 (MQ=255)                                                                                           cgAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCT‑CTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAAt                                          >  1:877669/1‑139 (MQ=255)                                                                                                        cccGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtgctgc                                                                                                                >  2:709411/1‑57 (MQ=255)                                                                                                        cccGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtgctgc                                                                                                                <  1:709411/57‑1 (MQ=255)                                                                                                             aTGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCT‑CTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCt                        <  2:877669/139‑1 (MQ=255)                                                                                                              tGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCT‑CTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTc                       >  1:95051/1‑139 (MQ=255)                                                                                                              tGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCT‑CTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTc                       <  1:71019/139‑1 (MQ=255)                                                                                                                 ggaggcggggAACAAAGAAGAGTAATTATTGGGAAGGGGC‑GTTCTGCCTCCTTTATTTAAGGTCATTT‑TATTAGTCCTTTGGTTATTCGGTTAATTTTAAAAAGAAATCTTAAAATTTTTATTCTTCGTTGGCTCCATa                   <  1:974730/132‑1 (MQ=255)                                                                                                                     ggCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCgg                                                                                                                      <  2:439504/38‑1 (MQ=255)                                                                                                                     ggCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCgg                                                                                                                      >  1:439504/1‑38 (MQ=255)                                                                                                                      gCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCT‑CTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACgtg               >  2:340512/1‑139 (MQ=255)                                                                                                                      gCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCT‑CTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACgtg               >  2:815777/1‑139 (MQ=255)                                                                                                                                   agagTAACTGATGGAGAGGGGCGGTGCTGCCT‑CTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACGTGTCTTGCTACATCa  <  1:1016705/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                      GACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAGCTGATGGAGAGGGGCGGTGCTGCCT‑CTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACGTGTCTTGCTACATCA  >  NC_000913/4078302‑4078570

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