Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_hsa_tpiA 2,106,741 A→T S162S (TCT→TCA)  wbbH ← O‑antigen polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA2,106,7410AT100.0% 141.1 / NA 42S162S (TCT→TCA) wbbHO‑antigen polymerase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (20/22);  total (20/22)

ACTCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTAGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAAC  >  NC_000913_3_hsa_tpiA/2106595‑2106882
                                                                                                                                                  |                                                                                                                                             
aCTCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTCCCTTTGTATTAGTAAATTTTTTTGa                                                                                                                                             >  2:163755/1‑149 (MQ=255)
  tctATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTATGTATTAGTAAATTTTTTTGAcc                                                                                                                                           >  1:259886/1‑149 (MQ=255)
    tATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAAGGTTTTACTTACCTTTCTATTAGTAAATTTTTTTAACCaa                                                                                                                                         >  2:364454/1‑149 (MQ=255)
        aaCACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTGCTTACCTTTTTATGAGTAAATTTTTTTGACCaaataa                                                                                                                                     >  1:460082/1‑149 (MQ=255)
             cTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAAtaaataaa                                                                                                                                >  1:619254/1‑149 (MQ=255)
                      aTGCATAAGAGATGATAAACATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAagta                                                                                                                       <  2:11424/149‑1 (MQ=255)
                        gCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAagtagt                                                                                                                     <  2:164886/149‑1 (MQ=255)
                           gaaGTGATGAGAAACATGAACACAATTAGCTTACCAGTATTTAGTATCAGTGCAAAGAGTAATACAATAAACACAACTAAACTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGATCAAATAAATCAAGCAAAAGTAGCTAg                                                                                                                  <  2:388496/148‑1 (MQ=255)
                                 atgatAACCATACAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGAt                                                                                                            <  2:561231/149‑1 (MQ=255)
                                     tAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTgcc                                                                                                        >  2:19773/1‑147 (MQ=255)
                                          aTAAAGACAATTTGCTTACCAGTATTGAGTATAATTGCAAAGATGAATACAATAAAAACAAGGAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCACATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCAt                                                                                                   <  2:640264/149‑1 (MQ=255)
                                           tAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATg                                                                                                  >  2:613694/1‑149 (MQ=255)
                                               gACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTTATCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGTCTGCATGTATg                                                                                              >  2:484913/1‑149 (MQ=255)
                                                   aTTTGCTTACGCGTATTCAGTATAATTGTCAAGATTAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATCAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTga                                                                                          <  1:43861/149‑1 (MQ=255)
                                                   aTTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTga                                                                                          >  2:60675/1‑149 (MQ=255)
                                                         ttACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGAAAACGTAGTTAGAATGATGTACTCCATGTATGCTCAGAAAtt                                                                                    >  2:43377/1‑149 (MQ=255)
                                                          aaCCAGTATTCAGTATAATTGCAAAGATTAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAAttt                                                                                   <  1:625645/148‑1 (MQ=255)
                                                            ccAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCt                                                                                 >  1:536017/1‑149 (MQ=255)
                                                            caaGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCt                                                                                 <  2:447199/147‑1 (MQ=255)
                                                             cAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGACTAGCTGCATGTATGCTGAGAAATTACtt                                                                                >  1:354317/1‑149 (MQ=255)
                                                             aaGTAATCAGTATAATTGCAAAGAAGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCtt                                                                                <  2:58190/148‑1 (MQ=255)
                                                                          aaTTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCaa                                                                   <  1:706359/149‑1 (MQ=255)
                                                                           aTTGCAAAGATGAATACAATAAAACCAAGTAAAGTAAAGGTCTCATCGACGTTTTGATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAAc                                                                  <  1:168543/149‑1 (MQ=255)
                                                                            ttGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACa                                                                 >  1:224394/1‑149 (MQ=255)
                                                                              gCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATGAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATc                                                               <  1:420514/149‑1 (MQ=255)
                                                                               cAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTATTTACATTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCa                                                              <  2:352097/149‑1 (MQ=255)
                                                                               cAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTAGGCTGAGAAATTTCTTTATGTTTCTTCAACATCa                                                              >  2:587298/1‑149 (MQ=255)
                                                                                             aTAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCaa                                                >  1:81323/1‑149 (MQ=255)
                                                                                                      agtaaagtaaATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCAtata                                       <  1:163755/149‑1 (MQ=255)
                                                                                                       gtaaagtaaaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTTAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAg                                      <  1:484913/149‑1 (MQ=255)
                                                                                                           agtaaaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCtac                                  >  2:693470/1‑148 (MQ=255)
                                                                                                             taaaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCtaagt                                <  2:401940/149‑1 (MQ=255)
                                                                                                                 tGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGc                            <  2:619254/149‑1 (MQ=255)
                                                                                                                   ttttACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTa                          <  2:224394/149‑1 (MQ=255)
                                                                                                                            ccTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAAc                 >  1:526998/1‑149 (MQ=255)
                                                                                                                             cTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACt                >  1:507830/1‑149 (MQ=255)
                                                                                                                               ttGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTgg              <  1:613694/149‑1 (MQ=255)
                                                                                                                                  tATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCCTGTATGATGAGAAATTTATTGCTGTGTCTTTAACATCAGTATCTTTTATCAATTTCAAACAGCTAAGTAAGCTAGTCCCGAACTGGTat           >  2:250252/1‑148 (MQ=255)
                                                                                                                                     tAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTa        <  1:568757/149‑1 (MQ=255)
                                                                                                                                        tAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACtt     >  1:605727/1‑149 (MQ=255)
                                                                                                                                        tAAATTTTTTTGACCAAATAAAGAAAGCAAAAGTAGTTAGAAGGATCGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACtt     <  1:43377/149‑1 (MQ=255)
                                                                                                                                           aTTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGTATCTATCTTAAGTTATATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAAc  <  1:175040/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
ACTCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTAGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAAC  >  NC_000913_3_hsa_tpiA/2106595‑2106882

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: ‑