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,3510GA94.6% 110.5 / ‑0.2 37coding (960/960 nt)srlEglucitol/sorbitol‑specific enzyme IIB component of PTS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (14/21);  minor base . (1/1);  total (15/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

TCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAGAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAG  >  NC_000913/2827226‑2827473
                                                                                                                             |                                                                                                                          
tCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATg                                                                                                               <  1:398979/139‑1 (MQ=255)
  gCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGAc                                                                                                             <  1:129039/139‑1 (MQ=255)
  gCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGAc                                                                                                             <  2:537393/139‑1 (MQ=255)
  gCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGAc                                                                                                             <  2:953012/139‑1 (MQ=255)
  gCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGAc                                                                                                             <  2:872902/139‑1 (MQ=255)
  gCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGAc                                                                                                             <  2:799369/139‑1 (MQ=255)
             cccGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATc                                                                                                  <  2:872670/139‑1 (MQ=255)
                gTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGCTGAAACATGACCGTTATTTATCAGa                                                                                               <  2:199609/138‑1 (MQ=255)
                gTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGCTGAAACATGACCGTTATTTATCAGa                                                                                               >  1:199609/1‑138 (MQ=255)
                      acacGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGAccacca                                                                                         <  2:130/139‑1 (MQ=255)
                                    ggTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGAccacca                                                                                         <  2:1006588/125‑1 (MQ=255)
                                    ggTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGAccacca                                                                                         >  1:1006588/1‑125 (MQ=255)
                                    ggTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGgc                                                                           >  1:778203/1‑139 (MQ=255)
                                            ttCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGc                                                                   <  2:778203/139‑1 (MQ=255)
                                             tCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGAc                                                                                                             >  1:963593/1‑96 (MQ=255)
                                             tCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGAc                                                                                                             <  2:963593/96‑1 (MQ=255)
                                             tCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCAcc                                                                                    <  1:610282/121‑1 (MQ=255)
                                             tCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCAcc                                                                                    >  2:610282/1‑121 (MQ=255)
                                                 tACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACAt                                                                                                                <  1:701244/89‑1 (MQ=255)
                                                 tACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACAt                                                                                                                >  2:701244/1‑89 (MQ=255)
                                                  aCTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCAc                                                                                     <  1:107673/115‑1 (MQ=255)
                                                  aCTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCAc                                                                                     >  2:107673/1‑115 (MQ=255)
                                                         aGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGc                                                          >  2:298633/1‑135 (MQ=255)
                                                         aGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGc                                                          <  1:298633/135‑1 (MQ=255)
                                                         aGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTACCAATAAAGGCTGAAACa                                                                                                                 >  2:611536/1‑80 (MQ=255)
                                                         aGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTACCAATAAAGGCTGAAACa                                                                                                                 <  1:611536/80‑1 (MQ=255)
                                                                ttttAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGAcc                                               >  1:141308/1‑139 (MQ=255)
                                                                         gcgcgcCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGAc                                                                                                             <  1:80124/68‑1 (MQ=255)
                                                                         gcgcgcCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGAc                                                                                                             >  2:80124/1‑68 (MQ=255)
                                                                          cgcgcCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCtcat                                     <  2:1053107/139‑1 (MQ=255)
                                                                                           tCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGgcg                    >  1:435243/1‑139 (MQ=255)
                                                                                             gCCTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGgcgcg                  >  2:678802/1‑139 (MQ=255)
                                                                                               cTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCAtt                                                               >  1:878656/1‑92 (MQ=255)
                                                                                               cTGGTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCAtt                                                               <  2:878656/92‑1 (MQ=255)
                                                                                                 ggTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCt                                                       <  2:633123/98‑1 (MQ=255)
                                                                                                 ggTTTGTCTCCGGTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCt                                                       >  1:633123/1‑98 (MQ=255)
                                                                                                             gTTTTATCTATCAATAAAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAag  <  1:462955/139‑1 (MQ=255)
                                                                                                                             |                                                                                                                          
TCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAGAGGCTGAAACATGACCGTTATTTATCAGACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAG  >  NC_000913/2827226‑2827473

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