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

TGGCTTGCAGTCGATCCTGCGCGAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAGTGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACAC  >  NC_000913/3304326‑3304584                                                                                                                                  |                                                                                                                                  tGGCTTGCAGTCGATCCTGCGCGAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTa                                                                                                                          >  1:1389972/1‑139 (MQ=255)  ggCTTGCAGTCGATCCTGCGCGAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTaa                                                                                                                         <  2:148622/139‑1 (MQ=255)  ggCTTGCAGTCGATCCTGCGCGAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTaa                                                                                                                         <  1:243648/139‑1 (MQ=255)            cGATCCTGCGCGAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGc                                                                                                               >  2:396149/1‑139 (MQ=255)                ccTGCGCGAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCt                                                                                                                            >  1:725053/1‑122 (MQ=255)                ccTGCGCGAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCt                                                                                                                            <  2:725053/122‑1 (MQ=255)                   gcgcGAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATcc                                                                                                        >  2:1156693/1‑139 (MQ=255)                     gcgAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGa                                                                                                      >  2:1262503/1‑139 (MQ=255)                         aaCCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGcc                                                                                                              >  1:770539/1‑127 (MQ=255)                         aaCCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGcc                                                                                                              <  2:770539/127‑1 (MQ=255)                             gACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCaa                                                                                              <  1:1232129/139‑1 (MQ=255)                              aCGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAAc                                                                                             >  1:1475567/1‑139 (MQ=255)                               cGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                            <  1:1156693/139‑1 (MQ=255)                               cGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                            <  1:883226/139‑1 (MQ=255)                               cGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCAACAGGCGCGGCTGCATTCGTTTGAGCTAGCGATGGATGTTAAGGAAGAGTTGTTGGGATGATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                            <  1:1006996/139‑1 (MQ=255)                                   gATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGt                                                                                        <  1:3578/139‑1 (MQ=255)                                    aTGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTc                                                                                       <  1:675872/139‑1 (MQ=255)                                         aTTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACa                                                                                  <  2:1452453/139‑1 (MQ=255)                                              ggTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGa                                                                                                      >  2:1153431/1‑114 (MQ=255)                                              ggTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGa                                                                                                      <  1:1153431/114‑1 (MQ=255)                                              ggTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                            <  2:487984/124‑1 (MQ=255)                                              ggTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                            >  1:487984/1‑124 (MQ=255)                                                  aaCGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                            <  2:351318/120‑1 (MQ=255)                                                  aaCGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                            >  1:351318/1‑120 (MQ=255)                                                   aCGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                            >  2:1478511/1‑119 (MQ=255)                                                   aCGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                            <  1:1478511/119‑1 (MQ=255)                                                   aCGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGAt                                                                        >  2:83042/1‑139 (MQ=255)                                                        cAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATaagcg                                                                   >  1:718331/1‑139 (MQ=255)                                                           aTTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATaagcgaag                                                                <  2:1475567/139‑1 (MQ=255)                                                              ttCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAgcgc                                                             <  1:1579205/139‑1 (MQ=255)                                                                cGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCAt                                                           <  1:396149/139‑1 (MQ=255)                                                                cGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCAt                                                           <  2:1389972/139‑1 (MQ=255)                                                                   tcacctGCCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTACCGCGGGAAGACTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGGGAAGCGCATCAg                                                        <  2:632036/138‑1 (MQ=255)                                                                     accaGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACa                                                                                  >  1:1360031/1‑111 (MQ=255)                                                                     accaGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACa                                                                                  <  2:1360031/111‑1 (MQ=255)                                                                     accaGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGc                                                      <  1:290323/139‑1 (MQ=255)                                                                     accaGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTGGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGc                                                      <  2:1603138/139‑1 (MQ=255)                                                                        aGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCAtt                                                   <  1:83042/139‑1 (MQ=255)                                                                                cTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                            >  2:797083/1‑90 (MQ=255)                                                                                cTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                            <  1:797083/90‑1 (MQ=255)                                                                                            gAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTc                               <  2:1428804/139‑1 (MQ=255)                                                                                                  gCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCGGAGTTGGCAACGCTTCCTTACATAGCCCTGATAAGCGAAGCGCACCAGGCATTGGGTAGCCAGAAGAAAGGTCGGaaaa                         >  2:1319687/1‑139 (MQ=255)                                                                                                            ttAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCaa                                                                                              <  2:837542/60‑1 (MQ=255)                                                                                                            ttAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCaa                                                                                              >  1:837542/1‑60 (MQ=255)                                                                                                            ttAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                            >  1:1374802/1‑62 (MQ=255)                                                                                                            ttAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                            <  2:1374802/62‑1 (MQ=255)                                                                                                              aaGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGAt             <  2:1040686/139‑1 (MQ=255)                                                                                                              aaGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGAt             >  1:1040686/1‑139 (MQ=255)                                                                                                                   agagTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGc                                                      <  1:395833/93‑1 (MQ=255)                                                                                                                   agagTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGc                                                      >  2:395833/1‑93 (MQ=255)                                                                                                                    gagTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                            <  2:1016730/54‑1 (MQ=255)                                                                                                                    gagTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACa                                                                                  <  2:1444225/64‑1 (MQ=255)                                                                                                                    gagTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACa                                                                                  >  1:1444225/1‑64 (MQ=255)                                                                                                                    aagTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                            >  1:1016730/2‑54 (MQ=255)                                                                                                                     agTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACa                                                                                  >  2:1484105/1‑63 (MQ=255)                                                                                                                     agTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACa                                                                                  <  1:1484105/63‑1 (MQ=255)                                                                                                                     agTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGAt      >  2:867101/1‑139 (MQ=255)                                                                                                                      gttgttGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                            <  1:1104382/52‑1 (MQ=255)                                                                                                                      gttgttGGGATAATGAGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                            >  2:1104382/1‑52 (MQ=255)                                                                                                                         gttgGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATacac  >  2:1634524/1‑139 (MQ=255)                                                                                                                               aTAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                            >  2:1193152/1‑43 (MQ=255)                                                                                                                                  |                                                                                                                                  TGGCTTGCAGTCGATCCTGCGCGAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAGTGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACAC  >  NC_000913/3304326‑3304584

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