Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 860,157 C→T 20.0% intergenic (‑129/+17) ybiV ← / ← ybiW sugar phosphatase; fructose‑1‑P/ribose‑5‑P/glucose‑6‑P phosphatase/putative pyruvate formate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913860,1570CT20.0% 40.7 / 6.3 20intergenic (‑129/+17)ybiV/ybiWsugar phosphatase; fructose‑1‑P/ribose‑5‑P/glucose‑6‑P phosphatase/putative pyruvate formate lyase
Reads supporting (aligned to +/- strand):  ref base C (8/8);  new base T (2/2);  total (10/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-01

AATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTT  >  NC_000913/860068‑860296
                                                                                         |                                                                                                                                           
aataaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCata                                                                                                                   <  1:985141/116‑1 (MQ=255)
aataaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCata                                                                                                                   >  2:985141/1‑116 (MQ=255)
             cGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAGATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCgaggag                                                                                                                              <  4:672197/92‑1 (MQ=255)
             cGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCgaggag                                                                                                                              >  3:672197/1‑92 (MQ=255)
                         aTTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGtcgccg                                                                                                                                        <  4:583763/70‑1 (MQ=38)
                         aTTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGtcgccg                                                                                                                                        >  3:583763/1‑70 (MQ=39)
                                    aaTGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGt                                                                                          <  1:563521/105‑1 (MQ=255)
                                    aaTGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGt                                                                                          >  2:563521/1‑105 (MQ=255)
                                                 aTGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGt                                                                                                                >  3:508323/1‑70 (MQ=255)
                                                 aTGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGt                                                                                                                <  4:508323/70‑1 (MQ=255)
                                                       ttGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCg                                                                                                    <  3:132663/76‑1 (MQ=255)
                                                       ttGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCg                                                                                                    >  4:132663/1‑76 (MQ=255)
                                                       ttGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGGTATCGTCCTGAg                                                                                    >  3:369560/1‑92 (MQ=255)
                                                       ttGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGGTATCGTCCTGAg                                                                                    <  4:369560/92‑1 (MQ=255)
                                                                     gATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCgatgat                                                                                               >  3:516338/1‑67 (MQ=255)
                                                                     gATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCgatgat                                                                                               <  4:516338/67‑1 (MQ=255)
                                                                               ggcggcGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTc                                                                                                              >  1:528882/1‑42 (MQ=25)
                                                                               ggcggcGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTc                                                                                                              <  2:528882/42‑1 (MQ=25)
                                                                               ggcggcGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGa                            <  3:194859/124‑1 (MQ=255)
                                                                               ggcggcGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGa                            >  4:194859/1‑124 (MQ=255)
                                                                                 cgacgaCGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATa                    <  4:1020296/130‑1 (MQ=255)
                                                                                 cgacgaCGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATa                    >  3:1020296/1‑130 (MQ=255)
                                                                                       cgccgcCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATAATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGtttttt  <  4:232024/142‑1 (MQ=255)
                                                                                       cgccgcCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATAATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGtttttt  >  3:232024/1‑142 (MQ=255)
                                                                                         |                                                                                                                                           
AATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTT  >  NC_000913/860068‑860296

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