Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,235,581 T→A 16.7% intergenic (+174/‑325) pgi → / → yjbE glucosephosphate isomerase/extracellular polysaccharide production threonine‑rich protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,235,5810TA16.7% 60.1 / 6.8 24intergenic (+174/‑325)pgi/yjbEglucosephosphate isomerase/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base T (10/10);  new base A (2/2);  total (12/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

GCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTT  >  NC_000913/4235462‑4235674
                                                                                                                       |                                                                                             
gCCTGATAAGACGCGACCGCGTCGCATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGt                                                                                           >  2:960413/1‑124 (MQ=38)
gCCTGATAAGACGCGACCGCGTCGCATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGt                                                                                           <  1:960413/124‑1 (MQ=38)
                  gcgTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACATAACGTCGGCAtt                                                                                    <  2:917631/113‑1 (MQ=255)
                  gcgTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACATAACGTCGGCAtt                                                                                    >  1:917631/1‑113 (MQ=255)
                          tCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCAt                                                                                     <  1:1407181/104‑1 (MQ=255)
                          tCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCAt                                                                                     >  2:1407181/1‑104 (MQ=255)
                          tCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACATAACGTCGGCAtt                                                                                    >  2:1086936/1‑105 (MQ=255)
                          tCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACATAACGTCGGCAtt                                                                                    <  1:1086936/105‑1 (MQ=255)
                                                       cATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGtt                                                                                               <  1:1306868/65‑1 (MQ=255)
                                                       cATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGtt                                                                                               >  2:1306868/1‑65 (MQ=255)
                                                            tATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGt                                                                                           <  3:1105399/64‑1 (MQ=255)
                                                            tATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGt                                                                                           >  4:1105399/1‑64 (MQ=255)
                                                                  gCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGaa                                            <  4:805521/105‑1 (MQ=255)
                                                                  gCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGaa                                            >  3:805521/1‑105 (MQ=255)
                                                                                 cACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGaaa                                           >  4:488043/1‑91 (MQ=255)
                                                                                 cACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGaaa                                           <  3:488043/91‑1 (MQ=255)
                                                                                    gattgattAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAg                               <  2:887651/100‑1 (MQ=255)
                                                                                    gattgattAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAg                               >  1:887651/1‑100 (MQ=255)
                                                                                             aaTAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGaaa                                           <  2:371300/79‑1 (MQ=255)
                                                                                             aaTAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGaaa                                           >  1:371300/1‑79 (MQ=255)
                                                                                              aTAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCtt                                                              >  2:191231/1‑59 (MQ=255)
                                                                                              aTAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCtt                                                              <  1:191231/59‑1 (MQ=255)
                                                                                                      aataaCAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGAt             >  3:1321656/1‑100 (MQ=255)
                                                                                                      aataaCAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGAt             <  4:1321656/100‑1 (MQ=255)
                                                                                                         aaCAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCgtt  <  2:1144377/108‑1 (MQ=255)
                                                                                                         aaCAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCgtt  >  1:1144377/1‑108 (MQ=255)
                                                                                                                       |                                                                                             
GCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTT  >  NC_000913/4235462‑4235674

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