Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,304,455	G→A	*336* (TAG→TAA)	yhbW →	putative luciferase‑like monooxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,304,455	0	G	A	100.0%	159.3 / NA	50	*336* (TAG→TAA)	yhbW	putative luciferase‑like monooxygenase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/31);  total (19/31)

TGGCTTGCAGTCGATCCTGCGCGAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAGTGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGAT  >  NC_000913/3304326‑3304573                                                                                                                                  |                                                                                                                       tGGCTTGCAGTCGATCCTGCGCGAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTa                                                                                                               >  1:119603/1‑139 (MQ=255)  ggCTTGCAGTCGATCCTGCGCGAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTaa                                                                                                              <  1:1582743/139‑1 (MQ=255)      tGCAGTCGATCCTGCGCGAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAAcgcg                                                                                                          >  2:1499584/1‑139 (MQ=255)                   gcgcGAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATcc                                                                                             <  2:846879/139‑1 (MQ=255)                     gcgAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGa                                                                                           <  1:791598/139‑1 (MQ=255)                     gcgAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGa                                                                                           >  2:364274/1‑139 (MQ=255)                     gcgAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGa                                                                                           >  1:587878/1‑139 (MQ=255)                         aaCCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCgg                                                                                                         >  2:349893/1‑121 (MQ=255)                         aaCCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCgg                                                                                                         <  1:349893/121‑1 (MQ=255)                         aaCCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTAt                                                                                               >  2:495269/1‑131 (MQ=255)                         aaCCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTAt                                                                                               <  1:495269/131‑1 (MQ=255)                             gACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCaa                                                                                   <  1:746545/139‑1 (MQ=255)                               cGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                 <  2:551608/139‑1 (MQ=255)                               cGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                 <  2:1390563/139‑1 (MQ=255)                               cGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                 <  2:119603/139‑1 (MQ=255)                               cGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                 <  1:364274/139‑1 (MQ=255)                               cGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                 <  1:1594822/139‑1 (MQ=255)                                                            ttttCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATaagcgaagc                                                    <  1:1011174/139‑1 (MQ=255)                                                                cGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCAt                                                <  1:1463000/139‑1 (MQ=255)                                                                     accaGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGc                                           <  2:383552/139‑1 (MQ=255)                                                                             cgGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTa                                   <  1:1499584/139‑1 (MQ=255)                                                                               gCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAGGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAgg                                 <  1:413105/139‑1 (MQ=255)                                                                                cTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGgc                                >  2:602768/1‑139 (MQ=255)                                                                                   cATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGc                                                                                         >  2:593153/1‑79 (MQ=255)                                                                                   cATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGc                                                                                         <  1:593153/79‑1 (MQ=255)                                                                                           tGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGt                     >  2:756665/1‑139 (MQ=255)                                                                                             aGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCt                                                                               <  1:622678/79‑1 (MQ=255)                                                                                             aGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCt                                                                               >  2:622678/1‑79 (MQ=255)                                                                                                  gCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGc                                                                                                    >  2:584888/1‑53 (MQ=255)                                                                                                  gCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGc                                                                                                    <  1:584888/53‑1 (MQ=255)                                                                                                  gCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATaa              >  2:1115899/1‑139 (MQ=255)                                                                                                      tggatgTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGGATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAggga          <  1:352756/139‑3 (MQ=255)                                                                                                      tggatgTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCa          <  1:701082/139‑1 (MQ=255)                                                                                                         atgTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCg       <  1:1115899/139‑1 (MQ=255)                                                                                                         atgTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCg       <  1:602768/139‑1 (MQ=255)                                                                                                            ttAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCa                                          <  2:1212036/101‑1 (MQ=255)                                                                                                            ttAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCa                                          >  1:1212036/1‑101 (MQ=255)                                                                                                              aaGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGAt  <  1:223824/139‑1 (MQ=255)                                                                                                                    gagTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTc                                                                            <  2:265157/59‑1 (MQ=255)                                                                                                                    gagTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTc                                                                            >  1:265157/1‑59 (MQ=255)                                                                                                                    gagTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGc                                           <  1:385617/92‑1 (MQ=255)                                                                                                                    gagTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGc                                           >  2:385617/1‑92 (MQ=255)                                                                                                                     agTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTc                                                                            <  1:65677/58‑1 (MQ=255)                                                                                                                     agTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTc                                                                            >  1:1487893/1‑58 (MQ=255)                                                                                                                     agTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTc                                                                            >  2:65677/1‑58 (MQ=255)                                                                                                                     agTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTc                                                                            <  2:1487893/58‑1 (MQ=255)                                                                                                                       ttgttgGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGc      <  2:1155082/126‑1 (MQ=255)                                                                                                                       ttgttgGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGc      >  1:1155082/1‑126 (MQ=255)                                                                                                                            gggATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATtgt                                      <  2:1294360/89‑1 (MQ=255)                                                                                                                            gggATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATtgt                                      >  1:1294360/1‑89 (MQ=255)                                                                                                                                  |                                                                                                                       TGGCTTGCAGTCGATCCTGCGCGAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAGTGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGAT  >  NC_000913/3304326‑3304573

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: ‑