Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,936,333 T→C 8.0% V19A (GTC→GCC)  rbsB → D‑ribose ABC transporter periplasmic binding protein; ribose chemotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,936,3330TC 8.0% 143.3 / 3.1 50V19A (GTC→GCC) rbsBD‑ribose ABC transporter periplasmic binding protein; ribose chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (24/22);  new base C (2/2);  total (26/24)
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.43e-01

CAGATGATCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAAC  >  NC_000913/3936191‑3936475
                                                                                                                                              |                                                                                                                                              
cAGATGATCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGc                                                                                                                                          >  3:138333/1‑149 (MQ=255)
 aGATGATCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGc                                                                                                                                          <  2:94579/148‑1 (MQ=255)
 aGATGATCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGc                                                                                                                                          >  1:94579/1‑148 (MQ=255)
 aGATGATCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGCCAGTGc                                                                                                                                          >  7:135244/1‑148 (MQ=255)
 aGATGATCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGCCAGTGc                                                                                                                                          >  7:131820/1‑148 (MQ=255)
     gatCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATg                                                                                                                                     >  7:107398/1‑149 (MQ=255)
       tCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGc                                                                                                                                    >  5:194583/1‑148 (MQ=255)
       tCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGc                                                                                                                                    <  6:194583/148‑1 (MQ=255)
         gTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGAt                                                                                                                                 >  1:128945/1‑149 (MQ=255)
             aaGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGc                                                                                                                              >  7:29067/1‑148 (MQ=255)
                cGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCaaaa                                                                                                                          <  8:29067/149‑1 (MQ=255)
                        aTTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCAt                                                                                                                  >  4:211321/1‑149 (MQ=255)
                          tttGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATcg                                                                                                                >  8:244686/1‑149 (MQ=255)
                                  cGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCtggtg                                                                                                        <  5:148560/149‑1 (MQ=255)
                                        tGGTACACAACAAAAAGCAGTAATATCGTCTACAGGAGATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCAGCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTcc                                                                                                  <  1:58628/149‑1 (MQ=255)
                                          gTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCAc                                                                                                >  1:221306/1‑149 (MQ=255)
                                          gTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCAc                                                                                                >  8:125699/1‑149 (MQ=255)
                                          gTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCAc                                                                                                >  5:175629/1‑149 (MQ=255)
                                                aacaaAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTaa                                                                                          <  8:107398/149‑1 (MQ=255)
                                                    aaaaGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTaacaac                                                                                      <  8:211251/149‑1 (MQ=255)
                                                    aaaaGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTaacaac                                                                                      <  3:215755/149‑1 (MQ=255)
                                                          aGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTc                                                                                <  7:125699/149‑1 (MQ=255)
                                                            taataaCGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCtt                                                                              <  5:8073/149‑1 (MQ=255)
                                                                aaCGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGTGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTa                                                                          >  6:48074/1‑149 (MQ=255)
                                                                 aCGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTggtggt                                                                                                      >  7:237738/1‑120 (MQ=255)
                                                                 aCGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTggtggt                                                                                                      <  8:237738/120‑1 (MQ=255)
                                                                  cGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATc                                                                        <  4:20062/149‑1 (MQ=255)
                                                                    aCTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTCGTGGTCTCCACGCTTAACAACCCGTTCATTGTATCg                                                                       >  4:228995/1‑148 (MQ=255)
                                                                         aGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGa                                                                   >  1:191934/1‑147 (MQ=255)
                                                                         aGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCGTAACAACCCGTGCTTTGTCTCGCTGa                                                                   >  5:176574/1‑147 (MQ=255)
                                                                                cTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGTGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGgcg                                                          >  2:138042/1‑149 (MQ=255)
                                                                                 ttGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGgcgc                                                         >  2:25897/1‑149 (MQ=255)
                                                                                        tGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGaaa                                                  <  3:137095/149‑1 (MQ=255)
                                                                                        tGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGaaa                                                  <  2:128945/149‑1 (MQ=255)
                                                                                           aCATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAagag                                               <  7:244686/149‑1 (MQ=255)
                                                                                           aCATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAagag                                               >  7:198084/1‑149 (MQ=255)
                                                                                             aTGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGc                                             >  8:69829/1‑149 (MQ=255)
                                                                                                         gCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGTGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACtt                                   <  1:138042/147‑1 (MQ=255)
                                                                                                         gCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACtt                                   <  7:69829/147‑1 (MQ=255)
                                                                                                         gCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGCCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTgg                                 <  8:131820/149‑1 (MQ=255)
                                                                                                         gCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGCCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTgg                                 <  8:135244/149‑1 (MQ=255)
                                                                                                             cccTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCtat                             <  2:221306/149‑1 (MQ=255)
                                                                                                                  gTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCt                        <  3:211321/149‑1 (MQ=255)
                                                                                                                   tttCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCtg                       <  1:25897/149‑1 (MQ=255)
                                                                                                                              gcgcTAAGTGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTgg             <  5:48074/148‑1 (MQ=255)
                                                                                                                               cgcTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGAc           >  7:77770/1‑149 (MQ=255)
                                                                                                                                gcTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACt          <  2:191934/149‑1 (MQ=255)
                                                                                                                                     gcgcCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTcc        >  5:101904/1‑146 (MQ=255)
                                                                                                                                     gcgcCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTcc        >  3:24317/1‑146 (MQ=255)
                                                                                                                                         cACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGaac  >  5:137886/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                              
CAGATGATCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAAC  >  NC_000913/3936191‑3936475

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