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

GCAGTCGATCCTGCGCGAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAGTGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATA  >  NC_000913/3304332‑3304581                                                                                                                            |                                                                                                                               gCAGTCGATCCTGCGCGAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCgg                                                                                                                 <  2:148285/139‑1 (MQ=255)                   aaCCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGa                                                                                                   <  1:660774/135‑1 (MQ=255)                   aaCCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGa                                                                                                   >  2:660774/1‑135 (MQ=255)                   aaCCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTg                                                                                               >  2:1018094/1‑139 (MQ=255)                         cGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                         <  2:258809/139‑1 (MQ=255)                         cGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                         <  2:178834/139‑1 (MQ=255)                            cGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTg                                                                                      <  1:860795/139‑1 (MQ=255)                              aTGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTc                                                                                    <  2:243857/139‑1 (MQ=255)                               tGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTcc                                                                                   <  1:436013/139‑1 (MQ=255)                                   aTTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACa                                                                               >  1:869546/1‑139 (MQ=255)                                      aTGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAg                                                                            >  2:1355906/1‑139 (MQ=255)                                        ggTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGAc                                                                          >  2:72484/1‑139 (MQ=255)                                            aaCGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                         >  1:1114120/1‑120 (MQ=255)                                            aaCGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                         <  2:1114120/120‑1 (MQ=255)                                            aaCGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGa                                                                      >  2:836971/1‑139 (MQ=255)                                             aCGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACa                                                                               >  2:661618/1‑129 (MQ=255)                                             aCGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACa                                                                               <  1:661618/129‑1 (MQ=255)                                             aCGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGAt                                                                     >  2:566929/1‑139 (MQ=255)                                                 gCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATaagc                                                                 >  1:376003/1‑139 (MQ=255)                                                      ttttCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATaagcgaagc                                                            <  1:566929/139‑1 (MQ=255)                                                           gACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATc                                                       <  2:429007/139‑1 (MQ=255)                                                           gACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATc                                                       <  1:1018094/139‑1 (MQ=255)                                                           gACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATc                                                       >  1:485335/1‑139 (MQ=255)                                                               accaGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGc                                                   <  1:1355906/139‑1 (MQ=255)                                                               accaGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCCGGAAGCCTTATCCGATCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGc                                                   <  2:869546/139‑1 (MQ=255)                                                                ccaGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                         >  2:1289956/1‑100 (MQ=255)                                                                ccaGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                         <  1:1289956/100‑1 (MQ=255)                                                                       cgGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTa                                           <  2:376003/139‑1 (MQ=255)                                                                           tGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGgca                                       >  2:665490/1‑139 (MQ=255)                                                                                       aGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATg                              <  1:942223/136‑1 (MQ=255)                                                                                       aGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATg                              >  2:942223/1‑136 (MQ=255)                                                                                            gCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGc                                                   >  1:549266/1‑110 (MQ=255)                                                                                            gCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGc                                                   <  2:549266/110‑1 (MQ=255)                                                                                                   atgTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCg               <  1:665490/139‑1 (MQ=255)                                                                                                            aaGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCTCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCg               <  1:966285/130‑1 (MQ=255)                                                                                                            aaGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCCTCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGAt          <  1:578040/135‑1 (MQ=255)                                                                                                            aaGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCg               >  2:966285/1‑130 (MQ=255)                                                                                                            aaGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGAt          >  2:578040/1‑135 (MQ=255)                                                                                                            aaGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACt      <  2:697771/139‑1 (MQ=255)                                                                                                            aaGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAACCTTATCCGAGCTGGCAACg                                                                                         <  2:743908/56‑1 (MQ=255)                                                                                                            aaGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAACCTTATCCGAGCTGGCAACg                                                                                         >  1:743908/1‑56 (MQ=255)                                                                                                              gagTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGc                                                                                                 >  1:87925/1‑46 (MQ=255)                                                                                                              gagTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGc                                                                                                 <  2:87925/46‑1 (MQ=255)                                                                                                               agTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGt                                                                                     <  2:43142/57‑1 (MQ=255)                                                                                                               agTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGt                                                                                     >  1:43142/1‑57 (MQ=255)                                                                                                               agTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGa                                  <  2:1064249/108‑1 (MQ=255)                                                                                                               agTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGa                                  >  1:1064249/1‑108 (MQ=255)                                                                                                                gttgttGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATa  >  2:683093/1‑139 (MQ=255)                                                                                                                gttgttGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATa  >  2:893222/1‑139 (MQ=255)                                                                                                                gttgttGGGATAATGTGTCTTAACGCGGGAAGCCTTAACCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATa  >  2:402282/1‑139 (MQ=255)                                                                                                                 ttgttgGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                         >  2:757821/1‑51 (MQ=255)                                                                                                                 ttgttgGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                         <  1:757821/51‑1 (MQ=255)                                                                                                                   gttgGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCt             <  1:905966/125‑1 (MQ=255)                                                                                                                   gttgGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCt             >  2:905966/1‑125 (MQ=255)                                                                                                                            |                                                                                                                               GCAGTCGATCCTGCGCGAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAGTGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATA  >  NC_000913/3304332‑3304581

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