Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,157,900 +A coding (32/1434 nt) ptsG → fused glucose‑specific PTS enzymes: IIB component/IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,157,8961.A93.8% 88.4 / ‑0.2 33coding (28/1434 nt)ptsGfused glucose‑specific PTS enzymes: IIB component/IIC component
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base A (11/19);  total (13/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.36e-01

TCTGTGTAATAAATAAAGGGCGCTTAGATGCCCTGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAG  >  NC_000913/1157758‑1158019
                                                                                                                                           |                                                                                                                            
tctGTGTAATAAATAAAGGGCGCTTAGATGCCCTGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCaaaa                                                                                                                         >  1:107025/1‑143 (MQ=255)
actGAGTAATAAATAAAACGCGCTTAGATGCCCTGTACACGGCGAGACTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCaaaa                                                                                                                         <  2:541945/142‑1 (MQ=255)
  tgtgTAATAAATAAAGGGCGCTTAGATGCCCTGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCaaaaag                                                                                                                       >  2:581976/1‑141 (MQ=255)
    ggtAATCAATAAATCGCTCTGAGGTGCCCTGTACGCGGCGAGGCTCTGCCCCTTCGTCAGACGTGAGAACGTAAAACAAGCATCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCaaaaaggt                                                                                                                     <  2:626729/142‑5 (MQ=255)
       aataaataAAGGGCGCTTAGATGCCCTGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCgg                                                                                                                   >  1:344949/1‑143 (MQ=255)
         taaataaaGGGCGCTTAGATGCCCTGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTa                                                                                                                 >  1:349304/1‑143 (MQ=255)
          aaataaaGGGCGCTTAGATGCCCTGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCACTTATGCTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTaa                                                                                                                >  2:332037/1‑143 (MQ=255)
               aaGGGCGCTTAGATGCCCTGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGc                                                                                                           <  1:662319/143‑1 (MQ=255)
                              cccTGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACGCAGGAGCACTCTCAATTATGTTTCAGAATGCATTTGCTAACCTGCGAAAAGGTCGGTAAATCGCTGATGCTGCCGGTAt                                                                                            >  1:426598/1‑143 (MQ=255)
                                 tGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGATAAATCGCTGATGCTGCCGGTATGCg                                                                                         >  2:158600/1‑143 (MQ=255)
                                   tACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTa                                                                                       <  1:222594/143‑1 (MQ=255)
                                    acacGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTACATCGCTGATGCTGCCGGTATCCGTAc                                                                                      >  2:652065/1‑143 (MQ=255)
                                    acacGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTAc                                                                                      <  2:543164/143‑1 (MQ=255)
                                           gAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCCCTGATGCTGCCCGTATGCGTACTGCCTAt                                                                               >  2:300062/1‑143 (MQ=255)
                                               ctctCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCa                                                                           >  1:190311/1‑143 (MQ=255)
                                                tctcCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATAGCTGATGCTGCCGGTATCCGTACCGCCTATCGCAg                                                                          >  1:109423/1‑143 (MQ=255)
                                                cctcCCCCCTTCCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAg                                                                          <  2:612620/142‑1 (MQ=255)
                                                         ttGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTctgc                                                                 <  2:344949/143‑1 (MQ=255)
                                                           gCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTctgctg                                                               <  2:711126/143‑1 (MQ=255)
                                                                 cgTGAGAACGTAAAAAAAGCACCCATCCTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTc                                                         <  2:441026/143‑1 (MQ=255)
                                                                   tGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCgg                                                       >  1:165115/1‑143 (MQ=255)
                                                                    gagaACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATCTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGt                                                      <  1:550600/143‑1 (MQ=255)
                                                                       taagTAAAAAATGCTGCCTTACACAAAAGCACGCTCAATTAAGTCTAAGAATGCATGTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTcc                                                   <  1:686470/140‑1 (MQ=255)
                                                                       aaCGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTcc                                                   <  1:581976/143‑1 (MQ=255)
                                                                           tAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGa                                               <  1:684906/143‑1 (MQ=255)
                                                                                     aCCCACACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAgctg                                     <  1:146062/143‑1 (MQ=255)
                                                                                        cATACTCGGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAGTCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAgctggct                                  <  2:429389/143‑1 (MQ=255)
                                                                                         aTACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAgctggctg                                 <  1:28939/143‑1 (MQ=255)
                                                                                          ggcgcAGGAGCGAGCTCAGTTAGGTTTAAGGAGGGAGTTGGTACCCTGCAAAGAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTGTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGc                                <  1:154892/139‑1 (MQ=255)
                                                                                              cAGGAGCACTCTCAATGATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCccgc                            <  2:694684/143‑1 (MQ=255)
                                                                                               aGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCccgcc                           <  2:349304/143‑1 (MQ=255)
                                                                                                  agccccctcAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGtt                        <  1:18170/137‑1 (MQ=255)
                                                                                                       tctcAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATc                   >  1:309615/1‑143 (MQ=255)
                                                                                                       tcgcAAGTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCAGCTAAGAGTCTGTTGCGGTAATTTCAGCTGGCTGCCCGCCGTTGTATc                   <  1:140227/140‑1 (MQ=255)
                                                                                                                  tttAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCGCGCCGTGGTATCTCATGTTATgc        >  1:607454/1‑142 (MQ=255)
                                                                                                                   ttAAGAATGCATTTGCTAACCTGCAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCCTGTTATGGCa      >  2:659965/1‑143 (MQ=255)
                                                                                                                        aaTGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAg  <  1:375880/143‑1 (MQ=255)
                                                                                                                                           |                                                                                                                            
TCTGTGTAATAAATAAAGGGCGCTTAGATGCCCTGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAG  >  NC_000913/1157758‑1158019

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: