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%	150.8 / NA	47	*330* (TAG→TAA)	yiiD	GNAT family putative N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/25);  total (22/25)

ACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAGCTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACGTGTCTTGCTAC  >  NC_000913/4078303‑4078566                                                                                                                                        |                                                                                                                                 aCGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTg                                                                                                                               >  2:1159348/1‑139 (MQ=255)      aaGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGga                                                                                                                          <  1:1323693/139‑1 (MQ=255)                   gATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGAAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtgct                                                                                                             >  1:676175/1‑139 (MQ=255)                     tGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtgctgc                                                                                                           <  1:1159348/139‑1 (MQ=255)                     tGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtgctgc                                                                                                           <  1:518256/139‑1 (MQ=255)                     tGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtgctgc                                                                                                           <  2:676175/139‑1 (MQ=255)                         ggTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCctct                                                                                                       >  1:1278146/1‑139 (MQ=255)                         ggTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCctct                                                                                                       >  1:192956/1‑139 (MQ=255)                          gTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCctctc                                                                                                      >  2:1191124/1‑139 (MQ=255)                              aaaTCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCAt                                                                                                  >  1:485164/1‑139 (MQ=255)                              aaaTCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCAt                                                                                                  >  1:459302/1‑139 (MQ=255)                               aaTCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCAtt                                                                                                 <  2:459302/139‑1 (MQ=255)                               aaTCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCAtt                                                                                                 <  2:485164/139‑1 (MQ=255)                               aaTCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTGCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCAtt                                                                                                 <  1:207905/139‑1 (MQ=255)                                   tttGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAgg                                                                                             <  1:175402/139‑1 (MQ=255)                                   tttGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAgg                                                                                             <  1:722871/139‑1 (MQ=255)                                                  ccGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTc                                                                              <  1:1318442/139‑1 (MQ=255)                                                        gCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGc                                                                        >  2:144044/1‑139 (MQ=255)                                                             gtTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTAttt                                                                   <  2:573668/139‑1 (MQ=255)                                                                tGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCAt                                                                                         >  2:758778/1‑114 (MQ=255)                                                                tGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCAt                                                                                         <  1:758778/114‑1 (MQ=255)                                                                tGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCgg                                                                <  1:1191124/139‑1 (MQ=255)                                                                     gCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTtata                                                           <  2:1278146/139‑1 (MQ=255)                                                                      cACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTTTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTatat                                                          <  1:219287/139‑1 (MQ=255)                                                                      cACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTatat                                                          >  1:546652/1‑139 (MQ=255)                                                                      cACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCCCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTatat                                                          <  2:762543/139‑1 (MQ=255)                                                                        cGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCg                                                                                                                  <  2:642781/81‑1 (MQ=255)                                                                        cGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCg                                                                                                                  >  1:642781/1‑81 (MQ=255)                                                                        cGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTa                                                        >  1:1148448/1‑139 (MQ=255)                                                                                 ttCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCAt                                                                                         >  1:1029479/1‑97 (MQ=255)                                                                                 ttCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCAt                                                                                         <  2:1029479/97‑1 (MQ=255)                                                                                    tGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAgggg                                                                                                                    <  2:865089/67‑1 (MQ=255)                                                                                    tGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAgggg                                                                                                                    >  1:865089/1‑67 (MQ=255)                                                                                      cccGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtgctgc                                                                                                           <  2:449875/74‑1 (MQ=255)                                                                                      cccGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtgctgc                                                                                                           >  1:449875/1‑74 (MQ=255)                                                                                      cccGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCtat                                          >  2:1348570/1‑139 (MQ=255)                                                                                         gcgAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTa                                                                                 >  1:649107/1‑97 (MQ=255)                                                                                         gcgAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTa                                                                                 <  2:649107/97‑1 (MQ=255)                                                                                                                     gCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTgg                                                                         <  1:582247/77‑1 (MQ=255)                                                                                                                     gCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTgg                                                                         >  2:582247/1‑77 (MQ=255)                                                                                                                     gCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACgtg           >  2:505917/1‑139 (MQ=255)                                                                                                                      cGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCAtt                                                                                                 <  1:1366079/52‑1 (MQ=255)                                                                                                                      cGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCAtt                                                                                                 >  2:1366079/1‑52 (MQ=255)                                                                                                                      cGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACgtgt          <  1:505917/139‑1 (MQ=255)                                                                                                                              aagaagAGTAACTGATGGAGAGGGGCGGTGCTGCctct                                                                                                       >  1:992354/1‑38 (MQ=255)                                                                                                                              aagaagAGTAACTGATGGAGAGGGGCGGTGCTGCctct                                                                                                       <  2:992354/38‑1 (MQ=255)                                                                                                                              aagaagAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACGTGTCTTGCTAc  >  2:507384/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                 ACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAGCTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACGTGTCTTGCTAC  >  NC_000913/4078303‑4078566

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