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%	123.9 / NA	40	*336* (TAG→TAA)	yhbW	putative luciferase‑like monooxygenase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/24);  total (16/24)

TCATGGCTTGCAGTCGATCCTGCGCGAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAGTGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCGG  >  NC_000913/3304323‑3304587                                                                                                                                     |                                                                                                                                     tCATGGCTTGCAGTCGATCCTGCGCGAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTc                                                                                                                                >  1:449777/1‑139 (MQ=255)       cTTGCAGTCGATCCTGCGCGAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAAcg                                                                                                                          >  1:405425/1‑139 (MQ=255)         tGCAGTCGATCCTGCGCGAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAAcgcg                                                                                                                        >  2:231747/1‑139 (MQ=255)           cAGTCGATCCTGCGCGAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCggg                                                                                                                      >  1:646835/1‑139 (MQ=255)                     tGCGCGAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATc                                                                                                            <  1:46650/139‑1 (MQ=255)                        gcgAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGa                                                                                                         <  1:423524/139‑1 (MQ=255)                        gcgAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGa                                                                                                         >  2:423524/1‑139 (MQ=255)                            aaCCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTg                                                                                                     >  1:97447/1‑139 (MQ=255)                                  cGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                               <  2:449777/139‑1 (MQ=255)                                  cGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                               <  1:136176/139‑1 (MQ=255)                                  cGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                               <  2:408257/139‑1 (MQ=255)                                       aTGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTc                                                                                          <  2:452422/139‑1 (MQ=255)                                       aTGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTc                                                                                          <  2:222321/139‑1 (MQ=255)                                                     aaCGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGa                                                                            >  2:649807/1‑139 (MQ=255)                                                      aCGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTCTCCTACATAGACCTGAt                                                                           >  1:474631/1‑139 (MQ=255)                                                                        accaGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGc                                                         <  2:754692/139‑1 (MQ=255)                                                                         ccaGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCa                                                        <  2:461277/139‑1 (MQ=255)                                                                         ccaGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCa                                                        <  2:97447/139‑1 (MQ=255)                                                                            ggCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATtg                                                     >  2:43346/1‑139 (MQ=255)                                                                                  gCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAgg                                               >  2:457893/1‑139 (MQ=255)                                                                                             ttGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATaagcgaa                                                                    <  2:607708/107‑1 (MQ=255)                                                                                             ttGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATaagcgaa                                                                    >  1:607708/1‑107 (MQ=255)                                                                                                  cTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCgg                                                                                                                       <  2:448323/51‑1 (MQ=255)                                                                                                  cTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCgg                                                                                                                       >  1:448323/1‑51 (MQ=255)                                                                                                     gCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGa                                                                                                         <  1:678942/62‑1 (MQ=255)                                                                                                     gCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGa                                                                                                         >  2:678942/1‑62 (MQ=255)                                                                                                         tggatgTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCa                        <  2:191797/139‑1 (MQ=255)                                                                                                         tggatgTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCa                        <  2:474631/139‑1 (MQ=255)                                                                                                            atgTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCg                     <  2:283727/139‑1 (MQ=255)                                                                                                             tgTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGc                    <  2:722617/139‑1 (MQ=255)                                                                                                             tgTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGc                    <  1:43346/139‑1 (MQ=255)                                                                                                               ttAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGa                                                                                                         >  1:670868/1‑52 (MQ=255)                                                                                                               ttAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGa                                                                                                         <  2:670868/52‑1 (MQ=255)                                                                                                                 aaGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGAt                <  1:457893/139‑1 (MQ=255)                                                                                                                     aaGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACt            <  1:575375/139‑1 (MQ=255)                                                                                                                         gttgttGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                               >  2:36998/1‑52 (MQ=255)                                                                                                                         gttgttGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                               <  1:36998/52‑1 (MQ=255)                                                                                                                             ttgGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTgg                                                                                                    >  1:652948/1‑43 (MQ=255)                                                                                                                             ttgGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTgg                                                                                                    <  2:652948/43‑1 (MQ=255)                                                                                                                               gggATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCgg  <  2:405425/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                     TCATGGCTTGCAGTCGATCCTGCGCGAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAGTGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCGG  >  NC_000913/3304323‑3304587

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