Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,936,381 C→T T35M (ACG→ATG)  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,3810CT100.0% 62.3 / NA 19T35M (ACG→ATG) rbsBD‑ribose ABC transporter periplasmic binding protein, ribose chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (8/11);  total (8/11)

AGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAACCCG  >  NC_000913/3936246‑3936481
                                                                                                                                       |                                                                                                    
aGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCt                                                                                                   >  1:126062/1‑139 (MQ=255)
               tACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGtt                                                                                       <  1:477574/136‑1 (MQ=255)
               tACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGtt                                                                                       >  2:477574/1‑136 (MQ=255)
                   ggACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTa                                                                                <  1:212154/139‑1 (MQ=255)
                                   aaCATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGgcg                                                                >  1:257487/1‑139 (MQ=255)
                                        caaaaaaCTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAg                                                           >  1:411728/2‑139 (MQ=255)
                                         aaaaaaCTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTaacaac                                                                                            <  1:138757/105‑1 (MQ=255)
                                         aaaaaaCTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTaacaac                                                                                            >  2:138757/1‑105 (MQ=255)
                                                   cTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGa                                                <  2:550571/139‑1 (MQ=255)
                                                         tGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACt                                          <  2:126062/139‑1 (MQ=255)
                                                            tttCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACGCCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTgg                                       <  2:257487/139‑1 (MQ=255)
                                                            tttCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTaacaac                                                                                            >  1:284116/1‑86 (MQ=255)
                                                            tttCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTaacaac                                                                                            <  2:284116/86‑1 (MQ=255)
                                                               ccGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCta                                    <  2:76325/139‑1 (MQ=255)
                                                                  cTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCTGATAAACTTGGCTATaa                                 <  2:411728/139‑1 (MQ=255)
                                                                                       tCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGc                                      <  1:58741/113‑1 (MQ=255)
                                                                                       tCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGc                                      >  2:58741/1‑113 (MQ=255)
                                                                                             cGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGaacaa      <  1:190590/139‑1 (MQ=255)
                                                                                                 tGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAACCCg  >  1:53879/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                    
AGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAACCCG  >  NC_000913/3936246‑3936481

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: