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,3510GA97.2% 109.5 / ‑5.1 36*320* (TAG→TAAsrlEglucitol/sorbitol‑specific enzyme IIB component of PTS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (16/19);  minor base C (0/1);  total (16/20)
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.94e-01

GTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAGAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTATTGCTT  >  NC_000913/2827217‑2827483
                                                                                                                                      |                                                                                                                                    
gtcggtTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGc                                                                                                                                  >  2:65899/1‑139 (MQ=255)
gtcggtTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGc                                                                                                                                  >  2:1068932/1‑139 (MQ=255)
     tttGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGaaa                                                                                                                             >  1:786143/1‑139 (MQ=255)
           gCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGAc                                                                                                                       <  1:65899/139‑1 (MQ=255)
              ggCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCga                                                                                                                    >  2:332967/1‑138 (MQ=255)
                  gAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGttatt                                                                                                                <  1:527649/139‑1 (MQ=255)
                     gCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGttatttat                                                                                                             <  2:911618/139‑1 (MQ=255)
                      cccGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATc                                                                                                            <  1:40937/139‑1 (MQ=255)
                                     ttCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGac                                                                                                        >  2:98106/1‑128 (MQ=255)
                                     ttCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGac                                                                                                        <  1:98106/128‑1 (MQ=255)
                                                      tCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGAc                                                                                                                       <  2:431654/96‑1 (MQ=255)
                                                      tCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGAc                                                                                                                       >  1:431654/1‑96 (MQ=255)
                                                         gTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCAtt                                                                         <  2:159718/139‑1 (MQ=255)
                                                           aCTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAACAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGa                                                                       <  2:758497/139‑1 (MQ=255)
                                                              ggTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTGTCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGc                                                                    <  2:195973/139‑1 (MQ=255)
                                                                 gAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTACCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCt                                                                 <  2:786143/139‑1 (MQ=255)
                                                                    ccGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATc                                                                                                 >  1:585532/1‑104 (MQ=255)
                                                                    ccGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATc                                                                                                 <  2:585532/104‑1 (MQ=255)
                                                                    ccGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGATATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAg                                                              <  1:332967/139‑1 (MQ=255)
                                                                      gCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTg                                                            >  2:215248/1‑139 (MQ=255)
                                                                         ttttAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGac                                                                                                        >  2:758103/1‑92 (MQ=255)
                                                                         ttttAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGac                                                                                                        <  1:758103/92‑1 (MQ=255)
                                                                                 ggCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGaccacc                                                                                                    >  2:163601/1‑88 (MQ=255)
                                                                                 ggCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGaccacc                                                                                                    <  1:163601/88‑1 (MQ=255)
                                                                                                    tCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGgcgc                                                                                   <  2:1045558/86‑1 (MQ=255)
                                                                                                    tCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGgcgc                                                                                   >  1:1045558/1‑86 (MQ=255)
                                                                                                     cGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGttatttat                                                                                                             <  2:1114460/59‑1 (MQ=255)
                                                                                                     cGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGttatttat                                                                                                             >  1:1114460/1‑59 (MQ=255)
                                                                                                      gCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGgcgcg                            >  1:441838/1‑139 (MQ=255)
                                                                                                        cTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGt                                     <  2:702718/128‑1 (MQ=255)
                                                                                                        cTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGt                                     >  1:702718/1‑128 (MQ=255)
                                                                                                        cTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGcc                          >  2:671820/1‑139 (MQ=255)
                                                                                                                tctcCGGTTTTATCTATCAATACAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGAcc                  <  1:984643/139‑1 (MQ=255)
                                                                                                                  tcCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACATTTCGTGAAGGCGCGCCTGCGGACCTc                >  1:1008109/1‑139 (MQ=255)
                                                                                                                         ttATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGt         <  2:81194/139‑1 (MQ=255)
                                                                                                                                tcAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTATTGCtt  <  1:215248/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                    
GTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAGAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTATTGCTT  >  NC_000913/2827217‑2827483

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: