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,7240AG95.7% 74.0 / ‑5.9 23T216A (ACT→GCT) ycgMputative isomerase/hydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (10/12);  minor base C (0/1);  total (10/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGATCGCTTATATGAGCAAGTTTTTTACCCTCAAGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAACTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCTGATGAGCGATGTACCTTTCTGGCA  >  NC_000913/1228592‑1228852
                                                                                                                                    |                                                                                                                                
cTGATCGCTTATATGAGCAAGTTTTTTACCCTCAAGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTcgcg                                                                                                                            <  2:650595/139‑1 (MQ=255)
  gATCGCTTATATGAGCAAGTTTTTTACCCTCAAGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGtt                                                                                                                          <  2:169258/139‑1 (MQ=255)
   aTCGCTTATATGAGCAAGTTTTTTACCCTCAAGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGttt                                                                                                                         >  1:58916/1‑139 (MQ=255)
     cGCTTATATGAGCAAGTTTTTTACCCTCAAGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTg                                                                                                                       <  1:1203661/139‑1 (MQ=255)
               aGCAAGTTTTTTACCCTCAAGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACtttt                                                                                                             >  2:238155/1‑139 (MQ=255)
                 cAAGTTTTTTACCCTCAAGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTg                                                                                                           <  1:238155/139‑1 (MQ=255)
                         ttACCCTCAAGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGa                                                                                                   >  1:1130433/1‑139 (MQ=255)
                                 aaGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAgcggcg                                                                                           >  2:391536/1‑139 (MQ=255)
                                       ggTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCaaaa                                                                                     >  2:430144/1‑139 (MQ=255)
                                            cGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGc                                                                                <  2:1130433/139‑1 (MQ=255)
                                                      aCAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCAGTGc                                                                      >  2:728044/1‑139 (MQ=255)
                                                         ggCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAg                                                                   <  1:928292/139‑1 (MQ=255)
                                                         ggCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAg                                                                   <  1:430144/139‑1 (MQ=255)
                                                          gCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGa                                                                  >  2:635676/1‑139 (MQ=255)
                                                          gCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGa                                                                  <  1:1034217/139‑1 (MQ=255)
                                                           ctcGCCTGATTGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCTCTTTCGATGCGCATTCTTTGACCCCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGAc                                                                 <  2:762942/137‑1 (MQ=255)
                                                           cACGCCTGATGGGGTCGTCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGCCAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGAc                                                                 <  2:687833/139‑1 (MQ=255)
                                                           cACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGAc                                                                 <  2:58916/139‑1 (MQ=255)
                                                                                                         gTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCTGATGAGc                   <  1:635676/139‑1 (MQ=255)
                                                                                                           cACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTGTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATAGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCTGATgagaga                 >  2:951316/1‑139 (MQ=255)
                                                                                                                tCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCTGATGAGCGATGTAc            >  1:1107289/1‑139 (MQ=255)
                                                                                                                tCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTACGTTTTAACGTAATGGCGGAACACCTGATGAGCGATGTAc            >  2:147409/1‑139 (MQ=255)
                                                                                                                          ttCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCTGATGAGCGATGTACCTTTCTGGCa  <  1:391536/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                
CTGATCGCTTATATGAGCAAGTTTTTTACCCTCAAGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAACTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCTGATGAGCGATGTACCTTTCTGGCA  >  NC_000913/1228592‑1228852

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: