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

AGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAGCTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACGTGTCTTGCTACATCA  >  NC_000913/4078309‑4078570                                                                                                                                  |                                                                                                                                     aGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGgag                                                                                                                             <  1:768127/139‑1 (MQ=255)           cAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtg                                                                                                                   >  1:1098623/1‑139 (MQ=255)               tGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtgctgc                                                                                                               <  1:1032670/139‑1 (MQ=255)               tGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtgctgc                                                                                                               <  1:164662/139‑1 (MQ=255)                     tCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtgctgc                                                                                                               >  1:552475/1‑133 (MQ=255)                     tCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATATGGCCCGGATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtgctgc                                                                                                               <  2:552475/133‑1 (MQ=255)                        aaaTCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAgggg                                                                                                                        <  2:694783/121‑1 (MQ=255)                        aaaTCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAgggg                                                                                                                        >  1:694783/1‑121 (MQ=255)                        aaaTCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCAt                                                                                                      >  1:575025/1‑139 (MQ=255)                                                 tGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtg                                                                                                                   <  2:71719/101‑1 (MQ=255)                                                 tGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtg                                                                                                                   >  1:71719/1‑101 (MQ=255)                                                      tgtTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTAtt                                                                        >  2:626441/1‑139 (MQ=255)                                                                cACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTatat                                                                                       >  2:471016/1‑114 (MQ=255)                                                                cACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTatat                                                                                       <  1:471016/114‑1 (MQ=255)                                                                  cGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCAt                                                                                             <  2:615565/106‑1 (MQ=255)                                                                  cGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCAt                                                                                             >  1:615565/1‑106 (MQ=255)                                                                  cGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTatat                                                                                       >  2:453662/1‑112 (MQ=255)                                                                  cGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTatat                                                                                       <  1:453662/112‑1 (MQ=255)                                                                  cGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTcc                                                                                 <  1:998919/118‑1 (MQ=255)                                                                  cGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTcc                                                                                 >  2:998919/1‑118 (MQ=255)                                                                         cGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGagag                                                                                                                           <  2:1046325/69‑1 (MQ=255)                                                                         cGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGagag                                                                                                                           >  1:1046325/1‑69 (MQ=255)                                                                           ttCTGCCCGCGAAGCCATTTGGTCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTAATTCAGGTCATTTATATTAGTCCTTCGCTATTTCGGTTATATTAGATAAGAAc                                                   <  2:924075/139‑1 (MQ=255)                                                                           ttCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGgag                                                                                                                             <  2:701951/65‑1 (MQ=255)                                                                           ttCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGgag                                                                                                                             >  1:701951/1‑65 (MQ=255)                                                                            tCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAgg                                                                                                 >  2:454892/1‑92 (MQ=255)                                                                            tCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAgg                                                                                                 <  1:454892/92‑1 (MQ=255)                                                                             cTGGCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGt                                                                   <  2:665073/121‑1 (MQ=255)                                                                             cTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtgctgc                                                                                                               <  2:1094862/77‑1 (MQ=255)                                                                             cTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtgctgc                                                                                                               >  1:1094862/1‑77 (MQ=255)                                                                             cTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGt                                                                   >  1:665073/1‑121 (MQ=255)                                                                              tGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTg                                                                                                                                   >  1:197693/1‑56 (MQ=255)                                                                              tGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTg                                                                                                                                   <  2:197693/56‑1 (MQ=255)                                                                              tGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAgggg                                                                                                                        >  2:537744/1‑67 (MQ=255)                                                                              tGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAgggg                                                                                                                        <  1:537744/67‑1 (MQ=255)                                                                               gCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtgctgc                                                                                                               >  1:457173/1‑75 (MQ=255)                                                                               gCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtgctgc                                                                                                               <  2:457173/75‑1 (MQ=255)                                                                                cccGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtg                                                                                                                   <  1:293394/70‑1 (MQ=255)                                                                                cccGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtg                                                                                                                   >  2:293394/1‑70 (MQ=255)                                                                                cccGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCAt                                                                                             <  2:1066122/92‑1 (MQ=255)                                                                                cccGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCAt                                                                                             >  1:1066122/1‑92 (MQ=255)                                                                                     gAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCAt                                                                                             <  1:497386/87‑1 (MQ=255)                                                                                     gAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCAt                                                                                             >  2:497386/1‑87 (MQ=255)                                                                                                     tATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGc                         <  1:60890/139‑1 (MQ=255)                                                                                                          agagGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCAt                    >  2:452483/1‑139 (MQ=255)                                                                                                           gagGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAgg                                                                                                 <  1:780732/61‑1 (MQ=255)                                                                                                           gagGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAgg                                                                                                 >  2:780732/1‑61 (MQ=255)                                                                                                               gCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACgtg               >  2:822043/1‑139 (MQ=255)                                                                                                               gCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACgtg               >  1:992407/1‑139 (MQ=255)                                                                                                                cGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtg                                                                                                                   <  1:152608/38‑1 (MQ=255)                                                                                                                cGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtg                                                                                                                   >  2:152608/1‑38 (MQ=255)                                                                                                                cGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTa                                                            >  1:130887/1‑93 (MQ=255)                                                                                                                cGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTa                                                            <  2:130887/93‑1 (MQ=255)                                                                                                                cGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATa                   <  1:25631/134‑1 (MQ=255)                                                                                                                cGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATa                   >  2:25631/1‑134 (MQ=255)                                                                                                                 gggAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACGTGTc             >  2:374877/1‑139 (MQ=255)                                                                                                                   gAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTC‑‑TAAATTTTTATTCTTCGtt                            <  1:334272/120‑1 (MQ=255)                                                                                                                   gAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTC‑‑TAAATTTTTATTCTTCGtt                            >  2:334272/1‑120 (MQ=255)                                                                                                                       gaagaaGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACGTGTCTTGCTa       <  1:822043/139‑1 (MQ=255)                                                                                                                            agagTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACGTGTCTTGCTACATCa  <  2:634871/139‑1 (MQ=255)                                                                                                                              agTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACGTGTCTTGCTACATCa  >  2:235699/1‑137 (MQ=255)                                                                                                                              agTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACGTGTCTTGCTACATCa  <  1:235699/137‑1 (MQ=255)                                                                                                                                  |                                                                                                                                     AGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAGCTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACGTGTCTTGCTACATCA  >  NC_000913/4078309‑4078570

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