Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,827,351 G→A *320* (TAG→TAA srlE → glucitol/sorbitol‑specific enzyme IIB component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,827,3510GA92.3% 63.6 / ‑1.0 26*320* (TAG→TAAsrlEglucitol/sorbitol‑specific enzyme IIB component of PTS
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (10/14);  total (11/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.95e-01

GTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAGAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTATTGCTTC  >  NC_000913/2827217‑2827484
                                                                                                                                      |                                                                                                                                     
gtcggtTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGc                                                                                                                                   <  1:1194094/139‑1 (MQ=255)
        gTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACAt                                                                                                                           <  1:196892/139‑1 (MQ=255)
           gCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGAc                                                                                                                        <  2:83292/139‑1 (MQ=255)
               gCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGtt                                                                                                                    >  2:162697/1‑139 (MQ=255)
                     gCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGttatttat                                                                                                              <  2:474471/139‑1 (MQ=255)
                                                         gTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTTGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCAtt                                                                          <  1:162697/139‑1 (MQ=255)
                                                         gTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCAtt                                                                          <  1:641662/139‑1 (MQ=255)
                                                               gTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAGAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGcc                                                                    >  2:521966/1‑139 (MQ=255)
                                                                           ttAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAGAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAg                                                        <  1:521966/139‑1 (MQ=255)
                                                                                  gcgcgcCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCg                                                                                                                      >  1:670296/1‑70 (MQ=255)
                                                                                  gcgcgcCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCg                                                                                                                      <  2:670296/70‑1 (MQ=255)
                                                                                  gcgcgcCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCAc                                                                                                <  2:1040560/92‑1 (MQ=255)
                                                                                  gcgcgcCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCAc                                                                                                >  1:1040560/1‑92 (MQ=255)
                                                                                   cgcgcCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCtcat                                                <  1:313740/139‑1 (MQ=255)
                                                                                              tACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTg                                     >  1:786301/1‑139 (MQ=255)
                                                                                                     cGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGgcgc                              >  2:962758/1‑139 (MQ=255)
                                                                                                      gCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGc                                                                              <  1:310677/90‑1 (MQ=255)
                                                                                                      gCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGc                                                                              >  2:310677/1‑90 (MQ=255)
                                                                                                        cTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGa                                                                                                          >  1:1662/1‑60 (MQ=255)
                                                                                                        cTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGa                                                                                                          <  2:1662/60‑1 (MQ=255)
                                                                                                        cTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGa                                                       <  1:1133251/111‑1 (MQ=255)
                                                                                                        cTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGa                                                       >  2:1133251/1‑111 (MQ=255)
                                                                                                                tctcCGGTTTTATCTATCAATAAAGGCTGAAACATGAc                                                                                                                        <  2:293365/38‑1 (MQ=255)
                                                                                                                tctcCGGTTTTATCTATCAATAAAGGCTGAAACATGAc                                                                                                                        >  1:293365/1‑38 (MQ=38)
                                                                                                                tctcCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGAcc                   <  1:905549/139‑1 (MQ=255)
                                                                                                                                 cAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTATTGCttc  >  2:386293/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                     
GTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAGAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTATTGCTTC  >  NC_000913/2827217‑2827484

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