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,3510GA93.8% 86.0 / ‑2.3 32*320* (TAG→TAAsrlEglucitol/sorbitol‑specific enzyme IIB component of PTS
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (16/14);  total (17/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 = 9.83e-01

GTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAGAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTATTGCTTC  >  NC_000913/2827217‑2827484
                                                                                                                                      |                                                                                                                                     
gtcggtTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGc                                                                                                                                   >  1:1409893/1‑139 (MQ=255)
           gCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGAc                                                                                                                        <  1:1377615/139‑1 (MQ=255)
           gCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGAc                                                                                                                        <  1:1569495/139‑1 (MQ=255)
           gCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGAc                                                                                                                        <  1:430940/139‑1 (MQ=255)
           gCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGAc                                                                                                                        <  2:604753/139‑1 (MQ=255)
             tGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCg                                                                                                                      <  2:1409893/139‑1 (MQ=255)
               gCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGtt                                                                                                                    >  1:967958/1‑139 (MQ=255)
                                     ttCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGacc                                                                                                        <  1:580010/129‑1 (MQ=255)
                                     ttCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGacc                                                                                                        >  2:580010/1‑129 (MQ=255)
                                         cgTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTAt                                                                                          >  2:879264/1‑139 (MQ=255)
                                                 tCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGa                                                                                  <  2:55447/139‑1 (MQ=255)
                                                           aCTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGtt                                                                                                                    >  2:662653/1‑95 (MQ=255)
                                                           aCTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGtt                                                                                                                    <  1:662653/95‑1 (MQ=255)
                                                              ggTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGc                                                                     <  2:47249/139‑1 (MQ=255)
                                                               gTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGcc                                                                    >  1:655732/1‑139 (MQ=255)
                                                                  aGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTc                                                                 >  1:1521604/1‑139 (MQ=255)
                                                                       cTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGaccacc                                                                                                     >  1:863027/1‑98 (MQ=255)
                                                                       cTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGaccacc                                                                                                     <  2:863027/98‑1 (MQ=255)
                                                                                 ggCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCtc                                                  >  2:1191911/1‑139 (MQ=255)
                                                                                                    tCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACg                                                                      >  2:644801/1‑100 (MQ=255)
                                                                                                    tCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACg                                                                      <  1:644801/100‑1 (MQ=255)
                                                                                                    tCGCCTGGTTTGTCTCCGATTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGATGgcg                               >  2:1060147/1‑139 (MQ=255)
                                                                                                      gCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGgcgcg                             >  1:656030/1‑139 (MQ=255)
                                                                                                        cTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGt                                                                                >  2:1366461/1‑86 (MQ=255)
                                                                                                        cTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGt                                                                                <  1:1366461/86‑1 (MQ=255)
                                                                                                        cTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTAGTTATCAGACCACCATCACCCGTATGGACGCGAGTGCCATTGACGGCCTCAGTGACCAGATGGGCATCAACTGTGGTGAAGgggcgcg                           >  1:1129169/1‑138 (MQ=255)
                                                                                                                   ccGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCg                <  1:1191911/139‑1 (MQ=255)
                                                                                                                          tatctatcAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTa         <  2:656030/139‑1 (MQ=255)
                                                                                                                          tatctatcAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTa         >  2:79016/1‑139 (MQ=255)
                                                                                                                                tcAATAGAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACTCCCTCAGTGACCAGATGCTCATCACCTTTCTTGAAgg                                 >  2:365305/1‑109 (MQ=255)
                                                                                                                                tcAATAGAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAgg                                 <  1:365305/109‑1 (MQ=255)
                                                                                                                                 cAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTATTGCttc  >  2:1502924/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 < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: