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% 146.8 / NA 42T35M (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 (20/22);  total (20/22)

GCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAACCCGGCGAAAGAGCTGGCGAACGTGCA  >  NC_000913/3936247‑3936504
                                                                                                                                      |                                                                                                                           
gCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCtt                                                                                                                         >  1:286754/1‑139 (MQ=255)
    taataaCGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTaaca                                                                                                                     >  2:420944/1‑139 (MQ=255)
     aataaCGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTaacaa                                                                                                                    >  2:520976/1‑139 (MQ=255)
      ataaCGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTaacaac                                                                                                                   <  1:119523/139‑1 (MQ=255)
             cTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCt                                                                                                            >  1:54444/1‑139 (MQ=255)
              tACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCtt                                                                                                           >  1:677089/1‑139 (MQ=255)
              tACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCtt                                                                                                           >  1:622490/1‑139 (MQ=255)
                cAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTg                                                                                                         <  2:286754/139‑1 (MQ=255)
                 aGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATTGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGt                                                                                                        >  1:441075/1‑139 (MQ=255)
                  ggACATCTTGAATATGAACATGAAAAAACTAGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTa                                                                                                       <  2:570733/139‑1 (MQ=255)
                                 gAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTa                                                                                                       <  1:82445/124‑1 (MQ=255)
                                 gAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTa                                                                                                       >  2:82445/1‑124 (MQ=255)
                                 gAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGgc                                                                                        >  1:588308/1‑139 (MQ=255)
                                 gAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGgc                                                                                        >  2:143180/1‑139 (MQ=255)
                                   aCATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGgcgc                                                                                      >  2:678898/1‑139 (MQ=255)
                                    cATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGgcgcg                                                                                     <  1:215926/139‑1 (MQ=255)
                                      tGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCa                                                                                   >  2:466397/1‑139 (MQ=255)
                                           aaaCTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAag                                                                              <  1:520976/139‑1 (MQ=255)
                                             aCTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAagag                                                                            <  2:161179/139‑1 (MQ=255)
                                                  cTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATTGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGa                                                                       <  2:441075/139‑1 (MQ=255)
                                                       cTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAAc                                                                  <  1:678898/139‑1 (MQ=255)
                                                       cTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAAc                                                                  <  2:1101489/139‑1 (MQ=255)
                                                            ttCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGc                                                             <  2:588308/139‑1 (MQ=255)
                                                              ccGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCta                                                           <  1:143180/139‑1 (MQ=255)
                                                               cGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCtat                                                          >  1:346262/1‑139 (MQ=255)
                                                                gCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAagag                                                                            >  2:153915/1‑120 (MQ=255)
                                                                gCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAagag                                                                            <  1:153915/120‑1 (MQ=255)
                                                                    ttGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCta                                                           >  1:335020/1‑133 (MQ=255)
                                                                    ttGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCta                                                           <  2:335020/133‑1 (MQ=255)
                                                                    ttGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCta                                                           <  2:40894/133‑1 (MQ=255)
                                                                    ttGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCta                                                           >  1:40894/1‑133 (MQ=255)
                                                                              cgcCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTa                                                                                                       >  2:202760/1‑79 (MQ=255)
                                                                              cgcCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTa                                                                                                       <  1:202760/79‑1 (MQ=255)
                                                                               gcCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTgg                                          <  1:523940/139‑1 (MQ=255)
                                                                                 cACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCtgg                                                   <  1:551949/128‑1 (MQ=255)
                                                                                 cACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCtgg                                                   >  2:551949/1‑128 (MQ=255)
                                                                                                  cGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAACCCGGc                       >  1:602754/1‑139 (MQ=255)
                                                                                                                 ccATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCtgg                                                   >  2:262110/1‑96 (MQ=255)
                                                                                                                 ccATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCtgg                                                   <  1:262110/96‑1 (MQ=255)
                                                                                                                     cgcgCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAACCCGGCGAAAGAGCTGGCGAACGTg    <  2:836098/139‑1 (MQ=255)
                                                                                                                     cgcgCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAACCCGGCGAAAGAGCTGGCGAACGTg    <  2:956516/139‑1 (MQ=255)
                                                                                                                       cgcTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAACCCGGCGAAAGAGCTGGCGAACGTGCa  <  1:1107238/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                           
GCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAACCCGGCGAAAGAGCTGGCGAACGTGCA  >  NC_000913/3936247‑3936504

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: