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%	187.2 / NA	53	S162S (TCT→TCA)	wbbH	O‑antigen polymerase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (31/22);  total (31/22)

TATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTAGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTC  >  NC_000913_3_hsa_tpiA/2106599‑2106885                                                                                                                                               |                                                                                                                                                 tATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCaa                                                                                                                                            >  2:503458/1‑149 (MQ=255)   ttAACACCTACGATGAATGCATAAGCGATGATAGGAATTAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAACTTTTTTGACAaaat                                                                                                                                          >  1:898268/1‑149 (MQ=255)      acacCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCGGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCaaataaa                                                                                                                                       >  1:592826/1‑149 (MQ=255)       cacCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCaaataaat                                                                                                                                      >  1:769040/1‑149 (MQ=255)       cacCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCaaataaat                                                                                                                                      >  2:664851/1‑149 (MQ=255)          cTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAAtaaataaa                                                                                                                                   <  1:552394/149‑1 (MQ=255)          cTACAATAAATACATAAGAGAGGTGAACCATAAAGACAAGTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAATAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAAtaaataaa                                                                                                                                   >  2:84172/1‑149 (MQ=255)              gATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACTAAATAAATAAAGCaa                                                                                                                               >  2:827787/1‑149 (MQ=255)                   acgCATAAGAGATGATAACCATACAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGGTGAATACAATAAAAACAAATAAAGTAAATCTTTTACTTACATTTGTATTAGTACATTTTTTTTACCAAATAAATAAAGCAAAagta                                                                                                                          >  2:626571/3‑149 (MQ=255)                       ataataGAAGATAACCATAAAGACAAGTAGCTTAACAGTAGTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTa                                                                                                                      <  2:858458/149‑1 (MQ=255)                          agagATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGaa                                                                                                                   <  2:609249/149‑1 (MQ=255)                            agaTGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATg                                                                                                                 >  2:742777/1‑149 (MQ=255)                             tatgatAACCATAAAGACAATTTTCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTCAATGTTTTGCTTACCCTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGa                                                                                                                <  1:764403/148‑1 (MQ=255)                               tgatAACCATAAAGACAATTTGCATACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGAtt                                                                                                              <  2:275968/149‑1 (MQ=255)                                   aaCCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCt                                                                                                          <  2:769040/149‑1 (MQ=255)                                       aTAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTAAGTTTGTATTAGTAAATTTTTTTGAACAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCAt                                                                                                      >  1:185771/1‑149 (MQ=255)                                           aGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTAt                                                                                                  >  1:376345/1‑149 (MQ=255)                                             aCAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGc                                                                                                >  2:337985/1‑149 (MQ=255)                                               aaTTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAAAATGATTGTCTGCTTATATGTTg                                                                                              >  2:296239/1‑149 (MQ=255)                                               aaTTTGCTTACCAGTATTCAGTATAATTGCAAAAATGAATACAATAAAAACAAGTAAAGTAAATGTTTTCTTTACCTTTGTATTAGCAAATTTTTTTGACCACAAAAATAAAGCAAAAGAAGTTAGAATGATAGGCTGCAGGTTTGcgt                                                                                              >  2:296238/1‑147 (MQ=255)                                                      ttACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAAtt                                                                                       >  2:398562/1‑149 (MQ=255)                                                       taacAGTATTGAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAAttt                                                                                      <  2:186456/146‑1 (MQ=255)                                                        aCCAGTATTCAGTATAATTGCAAATATGAATACAATAAAAACAAGTAAAGTAAATGGTTTACTTACCTTTGTATGAGTAAATTTTTTTGATCATTTTAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTc                                                                                     <  2:248371/149‑1 (MQ=255)                                                         ccAGTATTCAGTATAATTGCCTAGAGTTGTACAAGAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAGTAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCt                                                                                    <  2:185771/149‑1 (MQ=255)                                                         ccAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCt                                                                                    >  2:433578/1‑149 (MQ=255)                                                         ccAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTAACCAAATAAATAAAGCAAAAGTAGTTAAATTGATTGGCTGCATGTATGCTGAGAAATTTCt                                                                                    >  1:194335/1‑149 (MQ=255)                                                          cAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCtt                                                                                   <  1:146558/149‑1 (MQ=255)                                                                  aGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATTTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTc                                                                           <  1:682035/149‑1 (MQ=255)                                                                  aGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTc                                                                           <  2:601490/149‑1 (MQ=255)                                                                     ataATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTc                                                                        <  1:337985/149‑1 (MQ=255)                                                                          tGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACAt                                                                   <  2:509625/149‑1 (MQ=255)                                                                           gCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACATTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATc                                                                  >  2:603172/1‑149 (MQ=255)                                                                                      tACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTa                                                       <  1:421636/149‑1 (MQ=255)                                                                                           tAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCaaa                                                  >  1:799907/1‑149 (MQ=255)                                                                                           tAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCaaa                                                  <  1:503458/149‑1 (MQ=255)                                                                                             aaaaCAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATTAGCACCTCTTATCAAAtt                                                >  2:739636/1‑149 (MQ=255)                                                                                                aCAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTTTTTTTTTTAATAAAGAAAAAGTAATTAGAATGTTTAGCTGCTTGTATGCTGACAAATTTTTTGATTTTTCTTCAACATCAGCATCTCTTATTAAATTCat                                             >  2:897498/1‑149 (MQ=255)                                                                                                        agtaaaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCtaa                                     <  1:662800/149‑1 (MQ=255)                                                                                                        agtaaaTGGTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTCATCAAATTCATATAGCtaa                                     <  1:827665/149‑1 (MQ=255)                                                                                                          taaaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCtaagt                                   >  1:518320/1‑149 (MQ=255)                                                                                                            aaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCtaagtaa                                 >  1:587745/1‑149 (MQ=255)                                                                                                              tGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGATCTAATAAATAAATCAAAAGTATTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCtaagtaatc                               >  2:815140/1‑147 (MQ=255)                                                                                                              tGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGc                               <  1:626571/149‑1 (MQ=255)                                                                                                                ttttACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTa                             >  2:664102/1‑149 (MQ=255)                                                                                                                ttttACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTa                             >  1:634164/1‑149 (MQ=255)                                                                                                                      ttacCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCGTCTCTTATCAAAGTCATATAGCTAAGTAAGCTAGTCCCg                       >  1:298458/1‑149 (MQ=255)                                                                                                                      ttacCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCg                       >  1:508280/1‑149 (MQ=255)                                                                                                                      ttacCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCg                       >  2:648760/1‑149 (MQ=255)                                                                                                                      ttacCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCg                       >  1:539508/1‑149 (MQ=255)                                                                                                                             tGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGt                <  2:592826/149‑1 (MQ=255)                                                                                                                                    gTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACt         <  1:84172/149‑1 (MQ=255)                                                                                                                                     tAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTCCtt        >  2:341304/1‑149 (MQ=255)                                                                                                                                           ttttTGACCAAATAAATAAAGGAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTc  <  2:745637/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                 TATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTAGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTC  >  NC_000913_3_hsa_tpiA/2106599‑2106885

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