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_tpiA	2,106,741	0	A	T	100.0%	152.8 / NA	46	S162S (TCT→TCA)	wbbH	O‑antigen polymerase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (24/22);  total (24/22)

TCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTAGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAAC  >  NC_000913_3_hsa_tpiA/2106597‑2106882                                                                                                                                                 |                                                                                                                                              tctATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAACTAGGGTAGATGGGTAACTTACCTTTGTATTAGTAAATTTTTTTGAcc                                                                                                                                           >  1:643750/1‑149 (MQ=255)   tATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCaa                                                                                                                                         >  1:441145/1‑149 (MQ=255)           ccTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAataaataa                                                                                                                                 <  1:652499/149‑1 (MQ=255)                  tGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCaaaa                                                                                                                          >  1:258389/1‑149 (MQ=255)                     aTGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGATCAAATAATTAAAGCAAAagta                                                                                                                       >  2:366310/1‑149 (MQ=255)                       gCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAagtagt                                                                                                                     <  1:499418/149‑1 (MQ=255)                         aTAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTACAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTa                                                                                                                   >  1:498175/1‑149 (MQ=255)                                    gaaCCATAAAGACAAGTTGCTTACCAGTATTCAGTAGAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGc                                                                                                        <  2:774806/148‑1 (MQ=255)                                       ccATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTGGTTAGAATGATTGGCTGc                                                                                                     >  1:96041/1‑149 (MQ=255)                                        cATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCa                                                                                                    >  1:594919/1‑149 (MQ=255)                                          tAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATg                                                                                                  >  1:433284/1‑149 (MQ=255)                                                cAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTATCATAGAAGATTATGGTGACCACAAAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCt                                                                                            <  2:794251/149‑1 (MQ=255)                                                  aTTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCaaataaa                                                                                                                                    <  1:795159/107‑1 (MQ=255)                                                  aTTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCaaataaa                                                                                                                                    >  2:795159/1‑107 (MQ=255)                                                  aTTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTga                                                                                          <  2:441145/149‑1 (MQ=255)                                                   tttGCTTACCAGTGTTCATTATAATTGCAAAGATGAATACGGTAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTTTATTAGTAAATTTTTTTGACCAAATATATAAATCACAAGTAGTTAGAATTATTGGCTGCATGTatgatggt                                                                                         >  2:409763/1‑147 (MQ=255)                                                           ccAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAACTTTCt                                                                                 >  2:573864/1‑149 (MQ=255)                                                            cAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAAAAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCtt                                                                                >  1:120429/1‑149 (MQ=255)                                                               tATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCCAATAAATAAAGCAAAAGTAGTGAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGAt                                                                             <  1:738087/149‑1 (MQ=255)                                                                 ttCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGAtgt                                                                           >  1:9203/1‑149 (MQ=255)                                                                    aGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTc                                                                        >  2:598866/1‑149 (MQ=255)                                                                       ataATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACTGTTTTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTc                                                                     >  1:445266/1‑149 (MQ=255)                                                                              cAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCa                                                              <  1:216651/149‑1 (MQ=255)                                                                                 aGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTa                                                                                                                   <  1:696216/93‑1 (MQ=255)                                                                                 aGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTa                                                                                                                   >  2:696216/1‑93 (MQ=255)                                                                                  gATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCAt                                                          <  2:433284/149‑1 (MQ=255)                                                                                        tACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTc                                                                     >  2:354323/1‑132 (MQ=255)                                                                                        tACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTc                                                                     <  1:354323/132‑1 (MQ=255)                                                                                         aCAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAACCAAAAGAAGTTAAAATAATGAAATAAATTTATTCGGAGAACTTTATTACAGATTCTTCAACATCAGCATCACatat                                                   >  2:344117/1‑149 (MQ=255)                                                                                              aaaaaCAAGTAAAGTAAATGGTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATCGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTGTCaagt                                              <  2:722407/149‑3 (MQ=255)                                                                                                aaaaaaATAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTc                                            <  2:120429/149‑1 (MQ=255)                                                                                                 aaCAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCa                                           >  2:642187/1‑149 (MQ=255)                                                                                                   cAAGTAAAGTAAATGTCTTACTTACCTTTTTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCata                                         <  1:429532/149‑1 (MQ=255)                                                                                                      gtaaagaaaaTGTTTAACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAg                                      <  2:397775/149‑1 (MQ=255)                                                                                                       tagagtaaaAGTTTTACTTACCCCTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAATAGTTAGAATGTTTATCTGCATGTATGCAGAGAAATTACTTTATGTGACTTCTACATGATCTAATGTTATCAAATTCAtataac                                     >  2:788592/4‑147 (MQ=255)                                                                                                          agtaaaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAAGAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCtaa                                  <  1:518615/149‑1 (MQ=255)                                                                                                                     tacttacCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTc                       <  2:519811/149‑1 (MQ=255)                                                                                                                        ttacCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTCAGTAAGCTAGTCCCg                    >  2:52259/1‑149 (MQ=255)                                                                                                                                gTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTa            <  1:628776/149‑1 (MQ=255)                                                                                                                                aaaaaaGTAAATTTTTTTGAACAATTAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTa            <  2:93798/144‑1 (MQ=255)                                                                                                                                  aTTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGt          >  1:369509/1‑149 (MQ=255)                                                                                                                                   ttAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGtt         <  2:124567/149‑1 (MQ=255)                                                                                                                                      gTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTAAACTAAGAAATTACTTGATGTGTCTTCAACATCAACATCTCTTATCAAATTCATATACCTAAGTAAGCTAGTCCCGAACTGGTAGTTACt      >  2:163505/1‑149 (MQ=255)                                                                                                                                         aaTTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTaa   <  2:258389/149‑1 (MQ=255)                                                                                                                                         aaTTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTaa   >  1:314460/1‑149 (MQ=255)                                                                                                                                          aTTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAAc  <  1:353985/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                              TCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTAGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAAC  >  NC_000913_3_hsa_tpiA/2106597‑2106882

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑