Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,228,724 A→G T216A (ACT→GCT)  ycgM → putative isomerase/hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,228,7240AG100.0% 84.2 / NA 24T216A (ACT→GCT) ycgMputative isomerase/hydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (8/16);  total (8/16)

GATCGCTTATATGAGCAAGTTTTTTACCCTCAAGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAACTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCTGATGAGCGATGTACCTTTCTGGCA  >  NC_000913/1228594‑1228852
                                                                                                                                  |                                                                                                                                
gATCGCTTATATGAGCAAGTTTTTTACCCTCAAGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGtt                                                                                                                          <  2:394253/139‑1 (MQ=255)
                               aaGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAgcggcg                                                                                           >  2:14900/1‑139 (MQ=255)
                                              gtgtTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCg                                                                            <  2:2765/139‑1 (MQ=255)
                                                   gACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCtgtg                                                                       <  2:1000717/139‑1 (MQ=255)
                                                       ggCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAg                                                                   <  2:586292/139‑1 (MQ=255)
                                                        gCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGa                                                                  <  1:14900/139‑1 (MQ=255)
                                                        gCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGa                                                                  <  2:64162/139‑1 (MQ=255)
                                                         cACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGAc                                                                 <  1:912997/139‑1 (MQ=255)
                                                              cTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGtt                                                            >  1:519065/1‑139 (MQ=255)
                                                                     cGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGt                                                     <  1:329771/139‑1 (MQ=255)
                                                                      gTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTg                                                    <  1:920141/139‑1 (MQ=255)
                                                                                         gTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCgg                                 >  2:738727/1‑139 (MQ=255)
                                                                                           gatgaGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTCTAACGTAATGGCGGaa                               <  2:148678/139‑1 (MQ=255)
                                                                                                 cTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCtg                         >  2:998080/1‑139 (MQ=255)
                                                                                                    aCAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCg                                                                            >  2:634921/1‑85 (MQ=255)
                                                                                                    aCAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCg                                                                            <  1:634921/85‑1 (MQ=255)
                                                                                                    aCAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCtgatg                      <  1:738727/139‑1 (MQ=255)
                                                                                                      aGTCACTTTCGATGGGCATTCTTTGACAGCTcgcg                                                                                                                            <  1:872533/35‑1 (MQ=255)
                                                                                                      aGTCACTTTCGATGGGCATTCTTTGACAGCTcgcg                                                                                                                            >  2:872533/1‑35 (MQ=255)
                                                                                                                    ggCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCTGATGAGCGATGTACCTTTCt      >  1:127432/1‑139 (MQ=255)
                                                                                                                        ttCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCTGATGAGCGATGTACCTTTCTGGCa  <  1:998080/139‑1 (MQ=255)
                                                                                                                        ttCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCTGATGAGCGATGTACCTTTCTGGCa  <  2:519065/139‑1 (MQ=255)
                                                                                                                         tCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGc                                   >  1:85303/1‑105 (MQ=255)
                                                                                                                         tCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGc                                   <  2:85303/105‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                
GATCGCTTATATGAGCAAGTTTTTTACCCTCAAGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAACTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCTGATGAGCGATGTACCTTTCTGGCA  >  NC_000913/1228594‑1228852

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: