Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,183,270 C→A A128S (GCG→TCG)  yegX ← putative hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,183,2700CA80.0% 38.4 / 8.3 20A128S (GCG→TCG) yegXputative hydrolase
Reads supporting (aligned to +/- strand):  ref base C (4/0);  new base A (8/8);  total (12/8)
Fisher's exact test for biased strand distribution p-value = 1.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

GTGAGAAATCCACCGTTTGCAGAAATAATCTCGCCTGAACTGAAGCGGATACCGACGGGGAGAAATAATGATACGCCCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGTCCACCAGCTTTTCGCCTTCCGTC  >  NC_000913/2183195‑2183349
                                                                           |                                                                               
gtgAGAAATCCACCGTTTGCAGAAATAATCTCGCCTGAACTGAAGCGGATACCGACGGGGAGAAATAATGATACGACCCa                                                                             <  1:357613/80‑1 (MQ=255)
 tgAGAAATCCACCGTTTGCAGAAATAATCTCGCCTGAACTGAAGCGGATACCGACGGGGAGAAATAATGATACGACCCAc                                                                            <  1:357614/80‑1 (MQ=255)
          cACCGTTTGCAGAAATAATCTCGCCTGAACTGAAGCGGATACCGACGGGGAGAAATAATGATACGACCCACGCAGCAGGc                                                                   >  1:357634/1‑80 (MQ=255)
                          aatCTCGCCTGAACTGAAGCGGATACCGACGGGGAGAAATAATGATACGACCCACGCAGCAGGCCATTTTCGCGGCTTAGt                                                  >  2:357626/1‑81 (MQ=255)
                               cGCCTGAACTGAAGCGGATACCGACGGGGAGAAATAATGATACGACCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGc                                               >  2:357639/1‑79 (MQ=255)
                                     aactgaaGCGGATACCGACGGGGAGAAATAATGATACGACCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTAcg                                       <  1:357615/81‑1 (MQ=255)
                                     aactgaaGCGGATACCGACGGGGAGAAATAATGATACGACCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTAc                                        <  2:357616/80‑1 (MQ=255)
                                       ctgaaGCGGATACCGACGGGGAGAAATAATGATACGACCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTAcgc                                      >  1:357640/1‑80 (MQ=255)
                                          aaGCGGATACCGACGGGGAGAAATAATGATACGACCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGaa                                   >  1:357647/1‑80 (MQ=255)
                                           aGCGGATACCGACGGGGAGAAATAATGATACGACCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGaaa                                  >  1:357649/1‑80 (MQ=255)
                                                   ccGACGGGGAGAAATAATGATACGACCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAgggg                          <  1:357617/80‑1 (MQ=255)
                                                    cGACGGGGAGAAATAATGATACGACCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGTc                        <  1:357618/81‑1 (MQ=255)
                                                      aCGGGGAGAAATAATGATACGACCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGTcca                      <  1:357619/81‑1 (MQ=255)
                                                         gggAGAAATAATGATACGACCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGTccacca                   <  1:357620/81‑1 (MQ=255)
                                                             gaAATAATGATACGCCCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGTCCACCAGCtt               >  1:357660/1‑81 (MQ=255)
                                                             gaAATAATGATACGCCCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGTCCACCAGCtt               >  2:357646/1‑81 (MQ=255)
                                                              aaaTAATGATACGCCCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGTCCACCAGCttt              >  2:357661/1‑81 (MQ=255)
                                                               aataatGATACGACCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGTCCACCAGCtttt             >  2:357648/1‑81 (MQ=255)
                                                                        aCGACCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATCGGGGTCCACCAGCTTTTCGCCTTCCg    >  1:357669/1‑81 (MQ=255)
                                                                          gCCCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGTCCACCAGCTTTTCGCCTTCCGTc  >  2:357654/1‑81 (MQ=255)
                                                                           |                                                                               
GTGAGAAATCCACCGTTTGCAGAAATAATCTCGCCTGAACTGAAGCGGATACCGACGGGGAGAAATAATGATACGCCCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGTCCACCAGCTTTTCGCCTTCCGTC  >  NC_000913/2183195‑2183349

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