Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 164,205 A→G S701G (AGT→GGT)  hrpB → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913164,2050AG100.0% 141.1 / NA 41S701G (AGT→GGT) hrpBputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/20);  total (21/20)

GGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCG  >  NC_000913/164078‑164334
                                                                                                                               |                                                                                                                                 
ggCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATGGTGAATTGCCt                                                                                                                        <  1:405932/139‑1 (MQ=255)
   aGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATGGTGAATTGCCTgcg                                                                                                                     >  1:818302/1‑139 (MQ=255)
        tGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATGGTGAATTGCCTGCGCATTa                                                                                                                <  1:8030/139‑1 (MQ=255)
          gAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATGGTGAATTGCCTGCGCATTaca                                                                                                              <  2:573116/139‑1 (MQ=255)
                       tgcCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATGGTGAATTGCCTGCGCATTACACTGTGCCGACggg                                                                                                 >  1:263864/1‑139 (MQ=255)
                           acaTATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATGGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGc                                                                                             >  1:302643/1‑139 (MQ=255)
                                     ggCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATGGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCa                                                                                   >  2:940514/1‑139 (MQ=255)
                                      gCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATGGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCAt                                                                                  >  2:643840/1‑139 (MQ=255)
                                       cGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATGGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCAtt                                                                                 >  2:513083/1‑139 (MQ=255)
                                       cGTACATTCACTACGCGGCCTGAAATCACTCGACATTTAACAGGCACTACGCGGATTACTTGATTGGGGAATGCAGAAACGACTGGATGGTGAATTGCCTGCGCATTACACTGTGCCGAAGGGAAGCCGGATCGCCAtt                                                                                 >  1:459433/1‑139 (MQ=255)
                                        gTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATGGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTc                                                                                >  1:1057534/1‑139 (MQ=255)
                                         tACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATGGTGAATTGCCTGCGCATTACACTGTGCGGACGGGAAGCCgg                                                                                          >  2:625712/1‑128 (MQ=255)
                                         tACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATGGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCgg                                                                                          <  1:625712/128‑1 (MQ=255)
                                           cATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATGGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGtt                                                                             <  1:118520/139‑1 (MQ=255)
                                              tCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATGGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATc                                                                          <  1:503503/139‑1 (MQ=255)
                                               cACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATGGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCa                                                                         <  2:818302/139‑1 (MQ=255)
                                                      cgGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATGGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGAt                                                                  >  2:862409/1‑139 (MQ=255)
                                                        gCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGAATGCAGCAACGTCTGGATGGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAAc                                                               >  1:1097981/1‑139 (MQ=255)
                                                                   cTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATGGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCgcgc                                                     >  2:335990/1‑139 (MQ=255)
                                                                   cTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATGGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAAACCGCCCgcgc                                                     >  1:227737/1‑139 (MQ=255)
                                                                          tttATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATGGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGt                                              <  1:862409/139‑1 (MQ=255)
                                                                                 ggCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATGGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATg                                       <  2:263864/139‑1 (MQ=255)
                                                                                  gCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATGGTGAATTGc                                                                                                                          >  1:615387/1‑55 (MQ=255)
                                                                                  gCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATGGTGAATTGc                                                                                                                          <  2:615387/55‑1 (MQ=255)
                                                                                  gCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATGGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGc                                                                                     >  1:482503/1‑92 (MQ=255)
                                                                                  gCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATGGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGc                                                                                     <  2:482503/92‑1 (MQ=255)
                                                                                               tACTTGATTGGGAATGCAGCAACGTCTGGATGGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGc                        <  2:1097981/139‑1 (MQ=255)
                                                                                                aCTTGATTGGGGAATGCAGCAACGTCTGGATGGTGAATTGCCTGCGCATTACACTGTGCCGACGGGa                                                                                                <  1:914907/67‑1 (MQ=255)
                                                                                                aCTTGATTGGGGAATGCAGCAACGTCTGGATGGTGAATTGCCTGCGCATTACACTGTGCCGACGGGa                                                                                                >  2:914907/1‑67 (MQ=255)
                                                                                                aCTTGATTGGGGAATGCAGCAACGTCTGGATGGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTga                                            >  2:1082441/1‑119 (MQ=255)
                                                                                                aCTTGATTGGGGAATGCAGCAACGTCTGGATGGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTga                                            <  1:1082441/119‑1 (MQ=255)
                                                                                                     aTTGGGGAATGCAGCAACGTCTGGATGGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGc                   <  1:643840/139‑1 (MQ=255)
                                                                                                         gggAATGCAGCAACGTCTGGATGGTGAATTGCCTgc                                                                                                                      >  2:1056669/1‑36 (MQ=255)
                                                                                                         gggAATGCAGCAACGTCTGGATGGTGAATTGCCTgc                                                                                                                      <  1:1056669/36‑1 (MQ=37)
                                                                                                             aTGCAGCAACGTCTGGATGGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATc           <  1:301076/139‑1 (MQ=255)
                                                                                                                   caACGTCTGGATGGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCcgacga     <  1:335990/139‑1 (MQ=255)
                                                                                                                      cGTCTGGATGGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCg  >  2:599951/1‑139 (MQ=255)
                                                                                                                        tCTGGATGGTGAATTGCCTGCGCATTACACTGTGcc                                                                                                       <  1:847106/36‑1 (MQ=255)
                                                                                                                        tCTGGATGGTGAATTGCCTGCGCATTACACTGTGcc                                                                                                       >  2:847106/1‑36 (MQ=255)
                                                                                                                        tCTGGATGGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGccac                >  1:79455/1‑123 (MQ=255)
                                                                                                                        tCTGGATGGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGccac                <  2:79455/123‑1 (MQ=255)
                                                                                                                               |                                                                                                                                 
GGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCG  >  NC_000913/164078‑164334

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: