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.9% 90.9 / ‑0.2 34coding (28/1434 nt)ptsGfused glucose‑specific PTS enzymes: IIB component/IIC component
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base A (12/19);  total (14/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.39e-01

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