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

TCTGTGTAATAAATAAAGGGCGCTTAGATGCCCTGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGC‑AAA‑AGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAG  >  NC_000913/1157758‑1158019
                                                                                                                                           |                                                                                                                            
tctGTGTAATAAATAAAGGGCGCTTAGATGCCCTGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGC‑aaa‑a                                                                                                                         >  1:107025/1‑143 (MQ=255)
actGAGTAATAAATAAAACGCGCTTAGATGCCCTGTACACGGCGAGACTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGC‑aaa‑a                                                                                                                         <  2:541945/142‑1 (MQ=255)
  tgtgTAATAAATAAAGGGCGCTTAGATGCCCTGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGC‑aaa‑aag                                                                                                                       >  2:581976/1‑141 (MQ=255)
    ggtAATCAATAAATCGCTCTGAGGTGCCCTGTACGCGGCGAGGCTCTGCCCCTTCGTCAGACGTGAGAACGTAAAACAAGCATCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGC‑aaa‑aaggt                                                                                                                     <  2:626729/142‑5 (MQ=255)
       aataaataAAGGGCGCTTAGATGCCCTGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAA‑AGGTCgg                                                                                                                   >  1:344949/1‑143 (MQ=255)
         taaataaaGGGCGCTTAGATGCCCTGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAA‑AGGTCGGTa                                                                                                                 >  1:349304/1‑143 (MQ=255)
          aaataaaGGGCGCTTAGATGCCCTGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCACTTATGCTTAAGAATGCATTTGCTAACCTGCAAAA‑AGGTCGGTaa                                                                                                                >  2:332037/1‑143 (MQ=255)
               aaGGGCGCTTAGATGCCCTGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAA‑AGGTCGGTAAATCGc                                                                                                           <  1:662319/143‑1 (MQ=255)
                              cccTGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACGCAGGAGCACTCTCAATTATGTTTCAGAATGCATTTGCTAACCTGCGAAA‑AGGTCGGTAAATCGCTGATGCTGCCGGTAt                                                                                            >  1:426598/1‑143 (MQ=255)
                                 tGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAA‑AGGTCGATAAATCGCTGATGCTGCCGGTATGCg                                                                                         >  2:158600/1‑143 (MQ=255)
                                   tACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAA‑AGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTa                                                                                       <  1:222594/143‑1 (MQ=255)
                                    acacGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAA‑AGGTCGGTACATCGCTGATGCTGCCGGTATCCGTAc                                                                                      >  2:652065/1‑143 (MQ=255)
                                    acacGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAA‑AGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTAc                                                                                      <  2:543164/143‑1 (MQ=255)
                                           gAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAA‑AGGTCGGTAAATCCCTGATGCTGCCCGTATGCGTACTGCCTAt                                                                               >  2:300062/1‑143 (MQ=255)
                                               ctctCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAA‑AGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCa                                                                           >  1:190311/1‑143 (MQ=255)
                                                tctcCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAA‑AGGTCGGTAAATAGCTGATGCTGCCGGTATCCGTACCGCCTATCGCAg                                                                          >  1:109423/1‑143 (MQ=255)
                                                cctcCCCCCTTCCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAA‑AGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAg                                                                          <  2:612620/142‑1 (MQ=255)
                                                         ttGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAA‑AGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTctgc                                                                 <  2:344949/143‑1 (MQ=255)
                                                           gCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAA‑AGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTctgctg                                                               <  2:711126/143‑1 (MQ=255)
                                                                 cgTGAGAACGTAAAAAAAGCACCCATCCTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAA‑AGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTc                                                         <  2:441026/143‑1 (MQ=255)
                                                                   tGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAA‑AGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCgg                                                       >  1:165115/1‑143 (MQ=255)
                                                                    gagaACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATCTGCTAACCTGCAAAA‑AGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGt                                                      <  1:550600/143‑1 (MQ=255)
                                                                       taagTAAAAAATGCTGCCTTACACAAAAGCACGCTCAATTAAGTCTAAGAATGCATGTGCTAACCTGCAAAA‑AGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTcc                                                   <  1:686470/140‑1 (MQ=255)
                                                                       aaCGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAA‑AGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTcc                                                   <  1:581976/143‑1 (MQ=255)
                                                                           tAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAA‑AGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGa                                               <  1:684906/143‑1 (MQ=255)
                                                                                     aCCCACACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAA‑AGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAgctg                                     <  1:146062/143‑1 (MQ=255)
                                                                                        cATACTCGGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAA‑AGGTCGGTAAGTCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAgctggct                                  <  2:429389/143‑1 (MQ=255)
                                                                                         aTACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAA‑AGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAgctggctg                                 <  1:28939/143‑1 (MQ=255)
                                                                                          ggcgcAGGAGCGAGCTCAGTTAGGTTTAAGGAGGGAGTTGGTACCCTGC‑AAAGAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTGTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGc                                <  1:154892/139‑1 (MQ=255)
                                                                                              cAGGAGCACTCTCAATGATGTTTAAGAATGCATTTGCTAACCTGCAAAA‑AGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCccgc                            <  2:694684/143‑1 (MQ=255)
                                                                                               aGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAA‑AGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCccgcc                           <  2:349304/143‑1 (MQ=255)
                                                                                                  agccccctcAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAA‑AGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGtt                        <  1:18170/137‑1 (MQ=255)
                                                                                                       tctcAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAA‑AGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATc                   >  1:309615/1‑143 (MQ=255)
                                                                                                       tcgcAAGTATGTTTAAGAATGCATTTGCTAACCTGCAAAA‑AGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCAGCTAAGAGTCTGTTGCGGTAATTTCAGCTGGCTGCCCGCCGTTGTATc                   <  1:140227/140‑1 (MQ=255)
                                                                                                                  tttAAGAATGCATTTGCTAACCTGCAAAA‑AGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCGCGCCGTGGTATCTCATGTTATgc        >  1:607454/1‑142 (MQ=255)
                                                                                                                   ttAAGAATGCATTTGCTAACCTGC‑AAA‑AGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCCTGTTATGGCa      >  2:659965/1‑143 (MQ=255)
                                                                                                                        aaTGCATTTGCTAACCTGCAAAA‑AGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAg  <  1:375880/143‑1 (MQ=255)
                                                                                                                                           |                                                                                                                            
TCTGTGTAATAAATAAAGGGCGCTTAGATGCCCTGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGC‑AAA‑AGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAG  >  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: ‑