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

ACTCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTAGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTT  >  NC_000913_3_hsa_tpiA/2106595‑2106879                                                                                                                                                   |                                                                                                                                           aCTCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAATATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGa                                                                                                                                          <  1:1076495/149‑1 (MQ=255)   tccaTTAACACCTACGTTTAATGCATAAAAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATAAAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGAcc                                                                                                                                        <  1:7548/146‑1 (MQ=255)       ttAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCaaat                                                                                                                                    >  2:931709/1‑149 (MQ=255)          acacCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCaaataaa                                                                                                                                 <  1:991815/149‑1 (MQ=255)                   aTGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCaaa                                                                                                                        <  1:389520/149‑1 (MQ=255)                   aTGAATGCAAAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCaaa                                                                                                                        <  1:257872/149‑1 (MQ=255)                        tGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAagtag                                                                                                                   >  1:372278/1‑149 (MQ=255)                           aTAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTa                                                                                                                >  1:625105/1‑149 (MQ=255)                            tAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCaaa                                                                                                                                     <  1:109299/127‑1 (MQ=255)                            tAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCaaa                                                                                                                                     >  2:109299/1‑127 (MQ=255)                              agagATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGaa                                                                                                             <  2:372278/149‑1 (MQ=255)                               gagaTGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAAt                                                                                                            >  1:599075/1‑149 (MQ=255)                                 gatgatAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGa                                                                                                          <  2:763278/149‑1 (MQ=255)                                  aagatAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAGTAAAGCAAAAGTAGTTAGAATGAt                                                                                                         <  1:1014742/147‑1 (MQ=255)                                        aCCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTg                                                                                                   >  1:995111/1‑149 (MQ=255)                                         ccATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCGAAAGTAGTTAGAATGATTGGCTGc                                                                                                  <  1:1014746/149‑1 (MQ=255)                                         ccATAAAGACAATTTGCTTACCAGTATTAAGTATAATTGCAAAGATGAATTCAATAAAAACAAGTAAAGTGAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTTTCTAAATAACTATAGCAAAAGTTGTTAGTATGTTTggttcc                                                                                                  >  2:345263/1‑145 (MQ=255)                                            tAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATg                                                                                               >  1:635842/1‑149 (MQ=255)                                            tAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATg                                                                                               >  1:830655/1‑149 (MQ=255)                                                gACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAACTAAACGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCCAAAATAATTACACTTGTTTGCTGNATGtcta                                                                                           >  2:1147511/1‑146 (MQ=255)                                                  cAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCt                                                                                         <  2:625105/149‑1 (MQ=255)                                                   aaTGTGCATACCAGGATTCAGTATAAGGGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTg                                                                                        <  1:241841/149‑1 (MQ=255)                                                              cAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCtt                                                                             >  1:1064942/1‑149 (MQ=255)                                                              cAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCtt                                                                             >  1:598775/1‑149 (MQ=255)                                                               aGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGGTTGGCTGCATGTATGCTGAGAAATTTCTTg                                                                            >  2:72789/1‑149 (MQ=255)                                                                    tCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGAATGTATGCTGAGAAATTTCTTGAtttg                                                                       >  2:452326/1‑149 (MQ=255)                                                                        tataATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTg                                                                                                   <  2:820284/117‑1 (MQ=255)                                                                        tataATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTg                                                                                                   >  1:820284/1‑117 (MQ=255)                                                                         ataATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTc                                                                  <  1:76598/149‑1 (MQ=255)                                                                               gCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATc                                                            >  1:137975/1‑149 (MQ=255)                                                                                  aaGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGc                                                         <  1:825963/149‑1 (MQ=255)                                                                                  aaGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACATTTGTATTAGTAAATTTTTTTGAACAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGATGCATGTATGCTGAGAAATTTCTTGATTTGTGTTAAACATCAga                                                         >  2:839666/1‑148 (MQ=255)                                                                                         aTACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCtt                                                  >  1:666935/1‑149 (MQ=255)                                                                                           aCAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTAt                                                <  1:931709/149‑1 (MQ=255)                                                                                             aaTAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCaaaa                                                                                                                       >  1:238365/1‑76 (MQ=255)                                                                                             aaTAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCaaaa                                                                                                                       <  2:238365/76‑1 (MQ=255)                                                                                             aaTAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCa                                              >  1:1164785/1‑149 (MQ=255)                                                                                               tAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGAt                                                                          >  2:617309/1‑119 (MQ=255)                                                                                               tAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGAt                                                                          <  1:617309/119‑1 (MQ=255)                                                                                                 aaaaCAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAAtt                                          <  2:598775/149‑1 (MQ=255)                                                                                                  aaaCAAGTAAAGTAAATGTTTTACTTACCATTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTc                                         >  1:697653/1‑149 (MQ=255)                                                                                                     cAAATAAAATAAATTTTTTACTTACCTTTGTATTAGTAAATTTTTTTAACCAAATAAATAAAGCAAAAGTAGTTATAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCata                                      <  1:148264/149‑1 (MQ=255)                                                                                                        gtaaagtaaaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAg                                   <  1:997615/149‑1 (MQ=255)                                                                                                        gaaaagtaaaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAg                                   <  2:49922/147‑1 (MQ=255)                                                                                                             gtaaaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCAtata                                    >  1:428716/1‑143 (MQ=255)                                                                                                             gtaaaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCAtata                                    <  2:428716/143‑1 (MQ=255)                                                                                                             gtaaaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCtaag                              >  1:362217/1‑149 (MQ=255)                                                                                                             gtaaaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCtaag                              >  1:822616/1‑149 (MQ=255)                                                                                                              taaaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACAAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCtaagt                             >  1:230787/1‑149 (MQ=255)                                                                                                                aaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCtaagtaa                           >  1:624563/1‑149 (MQ=255)                                                                                                                 aTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGTAACAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATGAGCATCTCTTATCAAATTCATATAGCtaagtaag                          <  2:884557/149‑1 (MQ=255)                                                                                                                   gTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTc                                                                  >  1:1082569/1‑107 (MQ=255)                                                                                                                   gTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTc                                                                  <  2:1082569/107‑1 (MQ=255)                                                                                                                    ttttACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTa                       >  2:921316/1‑149 (MQ=255)                                                                                                                              cTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACt             >  1:1076797/1‑149 (MQ=255)                                                                                                                                     ttAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTGGTACCGAACTGgttgtt      >  2:350168/1‑149 (MQ=255)                                                                                                                                      tAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTa     <  1:646472/149‑1 (MQ=255)                                                                                                                                       aGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCCGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTAc    <  1:921316/149‑1 (MQ=255)                                                                                                                                         tAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACtt  >  1:896728/1‑149 (MQ=255)                                                                                                                                         tAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAATTACtt  >  1:974676/1‑149 (MQ=255)                                                                                                                                         tAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCCCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACtt  >  1:518000/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                           ACTCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTAGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTT  >  NC_000913_3_hsa_tpiA/2106595‑2106879

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