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_0009133,304,4550GA100.0% 168.5 / NA 53*336* (TAG→TAAyhbWputative luciferase‑like monooxygenase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/33);  total (20/33)

GGCTTGCAGTCGATCCTGCGCGAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAGTGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGAT  >  NC_000913/3304327‑3304573
                                                                                                                                |                                                                                                                      
ggCTTGCAGTCGATCCTGCGCGAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTGGGGATGATGTGTCTTaa                                                                                                              <  1:382822/139‑1 (MQ=255)
      cAGTCGATCCTGCGCGAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCggg                                                                                                        <  2:1529553/139‑1 (MQ=255)
            aTCCTGCGCGAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCt                                                                                                  >  2:928692/1‑139 (MQ=255)
              ccTGCGCGAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTa                                                                                                >  1:791984/1‑139 (MQ=255)
                    cgAAACCGACGCCGATGAGATTTTGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAg                                                                                          <  1:162284/139‑1 (MQ=255)
                       aaCCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTg                                                                                       >  2:949762/1‑139 (MQ=255)
                       aaCCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTg                                                                                       >  2:146588/1‑139 (MQ=255)
                           gACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCaa                                                                                   <  2:1009791/139‑1 (MQ=255)
                             cGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                 <  1:1360646/139‑1 (MQ=255)
                             cGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                 <  1:845888/139‑1 (MQ=255)
                             cGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                 <  1:1502028/139‑1 (MQ=255)
                             cGCCGATGAGATTATGGTCAACGGGCACATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                 <  1:928692/139‑1 (MQ=255)
                                 gATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGt                                                                             <  2:791984/139‑1 (MQ=255)
                                  aTGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTc                                                                            <  1:1095829/139‑1 (MQ=255)
                                  aTGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTc                                                                            <  2:571107/139‑1 (MQ=255)
                                       agtATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACa                                                                       <  1:949762/137‑1 (MQ=255)
                                        ttATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGa                                                                                           >  2:1488060/1‑118 (MQ=255)
                                        ttATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGa                                                                                           <  1:1488060/118‑1 (MQ=255)
                                            ggTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGAc                                                                  <  1:146588/139‑1 (MQ=255)
                                            ggTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGAc                                                                  <  1:761061/139‑1 (MQ=255)
                                                 aCGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCggg                                                                                                        <  1:413689/96‑1 (MQ=255)
                                                 aCGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGGTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCggg                                                                                                        >  2:413689/1‑96 (MQ=255)
                                                        gATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                 <  1:769948/112‑1 (MQ=255)
                                                        gATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                 >  2:769948/1‑112 (MQ=255)
                                                                   accaGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                 >  1:438110/1‑101 (MQ=255)
                                                                   accaGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                 <  2:438110/101‑1 (MQ=255)
                                                                   accaGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGc                                           <  1:1449500/139‑1 (MQ=255)
                                                                    ccaGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATaagcgaagc                                                    >  2:1518675/1‑129 (MQ=255)
                                                                    ccaGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATaagcgaagc                                                    <  1:1518675/129‑1 (MQ=255)
                                                                    ccaGGCGCGGCTGCATCCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCa                                          >  2:1572595/1‑139 (MQ=255)
                                                                               tGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                 >  2:875090/1‑89 (MQ=255)
                                                                               tGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                 <  1:875090/89‑1 (MQ=255)
                                                                                  aTTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGgcagca                            <  2:1071582/139‑1 (MQ=255)
                                                                                           aGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCg                   >  2:767192/1‑139 (MQ=255)
                                                                                                gCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCt                                                                               >  1:534039/1‑74 (MQ=255)
                                                                                                gCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCt                                                                               <  2:534039/74‑1 (MQ=255)
                                                                                                gCGATGGATGTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATaa              >  2:854520/1‑139 (MQ=255)
                                                                                                       atgTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCg       <  1:1287659/139‑1 (MQ=255)
                                                                                                        tgTTAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGc      <  1:854520/139‑1 (MQ=255)
                                                                                                          ttAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                 <  1:422821/62‑1 (MQ=255)
                                                                                                          ttAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                 >  2:422821/1‑62 (MQ=255)
                                                                                                          ttAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGc                                           <  1:719727/100‑1 (MQ=255)
                                                                                                          ttAAGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGc                                           >  2:719727/1‑100 (MQ=255)
                                                                                                            aaGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGAt  <  1:1572595/139‑1 (MQ=255)
                                                                                                            aaGGAAGAGTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGAt  <  1:969234/139‑1 (MQ=255)
                                                                                                                   agTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                 >  2:792072/1‑53 (MQ=255)
                                                                                                                   agTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACg                                                                                 <  1:792072/53‑1 (MQ=255)
                                                                                                                   agTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGc      <  1:1288564/128‑1 (MQ=255)
                                                                                                                   agTTGTTGGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGc      >  2:1288564/1‑128 (MQ=255)
                                                                                                                         tgGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACa                                                                       >  2:1021069/1‑57 (MQ=255)
                                                                                                                         tgGGATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACa                                                                       <  1:1021069/57‑1 (MQ=255)
                                                                                                                          gggATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTgg                                                                                      <  2:559279/41‑1 (MQ=255)
                                                                                                                          gggATAATGTGTCTTAACGCGGGAAGCCTTATCCGAGCTgg                                                                                      >  1:559279/1‑41 (MQ=255)
                                                                                                                                |                                                                                                                      
GGCTTGCAGTCGATCCTGCGCGAAACCGACGCCGATGAGATTATGGTCAACGGGCAGATTTTCGACCACCAGGCGCGGCTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAGTGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGAT  >  NC_000913/3304327‑3304573

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: