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

AAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAGCTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACGTGTCTTGCTACATCA  >  NC_000913/4078306‑4078570                                                                                                                                     |                                                                                                                                     aaaaGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATg                                                                                                                                <  1:375501/139‑1 (MQ=255)                gATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtgct                                                                                                                 <  1:106948/139‑1 (MQ=255)                  tGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtgctgc                                                                                                               <  2:535786/139‑1 (MQ=255)                               cTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAg                                                                                                  <  1:306022/139‑1 (MQ=255)                                                       aGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTc                                                                                               <  1:59539/118‑1 (MQ=255)                                                       aGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTc                                                                                               >  2:59539/1‑118 (MQ=255)                                                              gAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGt                                                                   <  1:278721/139‑1 (MQ=255)                                                                     cGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGgag                                                                                                                             >  2:46933/1‑74 (MQ=255)                                                                     cGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGgag                                                                                                                             <  1:46933/74‑1 (MQ=255)                                                                      gTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAgg                                                                                                 >  2:1117944/1‑101 (MQ=255)                                                                      gTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAgg                                                                                                 <  1:1117944/101‑1 (MQ=255)                                                                            cGTTCTGCCCGCGAAGCCGTTTGGCCCGTCTGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGCGCGGTGCTGCCTCTCTCGTTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGa                                                     >  1:714896/1‑139 (MQ=255)                                                                              ttCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCctct                                                                                                           <  1:1120186/83‑1 (MQ=255)                                                                              ttCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCctct                                                                                                           >  2:1120186/1‑83 (MQ=255)                                                                              ttCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCAt                                                                                             >  1:1218571/1‑97 (MQ=255)                                                                              ttCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCAt                                                                                             <  2:1218571/97‑1 (MQ=255)                                                                                cTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtgctgc                                                                                                               >  1:438732/1‑77 (MQ=255)                                                                                cTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtgctgc                                                                                                               <  2:438732/77‑1 (MQ=255)                                                                                cTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTAttt                                                                       <  2:919919/117‑1 (MQ=255)                                                                                cTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTAttt                                                                       >  1:919919/1‑117 (MQ=255)                                                                                           gCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCctc                                                                                                            >  2:607679/1‑69 (MQ=255)                                                                                           gCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCctc                                                                                                            <  1:607679/69‑1 (MQ=255)                                                                                                                 ggCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCt                        >  1:1084877/1‑131 (MQ=255)                                                                                                                 ggCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCt                        <  2:1084877/131‑1 (MQ=255)                                                                                                                 ggCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACgt                >  2:1031262/1‑139 (MQ=255)                                                                                                                 ggCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACgt                >  1:651248/1‑139 (MQ=255)                                                                                                                  gCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACgtg               >  2:1102971/1‑139 (MQ=255)                                                                                                                  gCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACgtg               >  1:1111147/1‑139 (MQ=255)                                                                                                                   cGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCtata                                             >  1:585559/1‑108 (MQ=255)                                                                                                                   cGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCtata                                             <  2:585559/108‑1 (MQ=255)                                                                                                                   cGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACgtgt              >  1:816476/1‑139 (MQ=255)                                                                                                                     ggAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACGTGTCt            <  2:651248/139‑1 (MQ=255)                                                                                                                        aCGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACGTGTCTTGc         >  2:517035/1‑139 (MQ=255)                                                                                                                           aagaagAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCAt                                                                                             <  1:1185156/52‑1 (MQ=255)                                                                                                                           aagaagAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCAt                                                                                             >  2:1185156/1‑52 (MQ=255)                                                                                                                           aagaagAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTatat                                                                                       <  1:928129/58‑1 (MQ=255)                                                                                                                           aagaagAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTatat                                                                                       >  2:928129/1‑58 (MQ=255)                                                                                                                            agaagaGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACGTGTCTTGCTACa     >  2:975109/1‑139 (MQ=255)                                                                                                                             gaagaGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACGTGTCTTGCTACAt    <  2:994219/139‑1 (MQ=255)                                                                                                                               agagTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACGTGTCTTGCTACATCa  <  2:1111147/139‑1 (MQ=255)                                                                                                                               agagTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACGTGTCTTGCTACATCa  <  2:816476/139‑1 (MQ=255)                                                                                                                                 agTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCAtt                                                                                            >  1:806508/1‑47 (MQ=255)                                                                                                                                 agTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCAtt                                                                                            <  2:806508/47‑1 (MQ=255)                                                                                                                                     |                                                                                                                                     AAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAGCTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACGTGTCTTGCTACATCA  >  NC_000913/4078306‑4078570

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