Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913_3_hsa_tpiA	2,106,741	A→T	100%	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%	175.3 / NA	54	S162S (TCT→TCA)	wbbH	O‑antigen polymerase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (29/25);  total (29/25)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTAGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTC  >  NC_000913_3_hsa_tpiA/2106596‑2106885                                                                                                                                                  |                                                                                                                                                 ctctATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGAc                                                                                                                                               <  1:315000/149‑1 (MQ=255)    tATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCaa                                                                                                                                            >  1:255261/1‑149 (MQ=255)    tATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGATTACCAGAACAAAGTTTGATAGCCAGAATGAATACAATAAAAACAAGTACAGTAAAATTTTTAGTTACCTCTGTATTAGTAAATTTTTTTGACCaa                                                                                                                                            >  1:829544/1‑149 (MQ=255)         acacCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCaaataaa                                                                                                                                       <  2:255261/149‑1 (MQ=255)             cTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAAtaaataaa                                                                                                                                   <  1:30606/149‑1 (MQ=255)              tACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCCCTAAAACGTCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAg                                                                                                                                  <  1:826246/149‑1 (MQ=255)              tACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAg                                                                                                                                  <  1:9108/149‑1 (MQ=255)                   tGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCaaaa                                                                                                                             <  1:122202/149‑1 (MQ=255)                         cATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGtt                                                                                                                       >  1:153308/1‑149 (MQ=255)                          aTAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTa                                                                                                                      >  1:830423/1‑149 (MQ=255)                           tAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTAGAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACGAAATAAATAAAGCAAAAGTAGTTAg                                                                                                                     >  1:758457/1‑149 (MQ=255)                             agagATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGaa                                                                                                                   >  1:512945/1‑149 (MQ=255)                              gagaTGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAAt                                                                                                                  >  1:420808/1‑149 (MQ=255)                                 atgatAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGAt                                                                                                               <  2:435890/149‑1 (MQ=255)                                  tgatAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTCTAATTCTCAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGAtt                                                                                                              <  2:771487/149‑1 (MQ=255)                                      aaCCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCt                                                                                                          >  1:493555/1‑149 (MQ=255)                                        ccATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGc                                                                                                        >  1:708612/1‑149 (MQ=255)                                         cATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCa                                                                                                       >  1:78368/1‑149 (MQ=255)                                              aGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTAt                                                                                                  >  1:830704/1‑149 (MQ=255)                                               gACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATg                                                                                                 >  1:487865/1‑149 (MQ=255)                                                  aaTTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTg                                                                                              >  2:589560/1‑149 (MQ=255)                                                  aaTTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATAATTGGCTGCATCTATGCTg                                                                                              >  2:376424/1‑149 (MQ=255)                                                   aTTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTga                                                                                             <  2:78368/149‑1 (MQ=255)                                                          tACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGATCAAATAAATACAACAAAAGTAGTTAGAATGATTGGCTGCATGTATCCTGAGAAAttt                                                                                      >  2:479965/1‑149 (MQ=255)                                                             cAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCtt                                                                                   >  1:457741/1‑149 (MQ=255)                                                                  ttCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGAtgt                                                                              <  1:476777/149‑1 (MQ=255)                                                                     aGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGATTAGAAATTTCTTGATATGTc                                                                           >  2:658387/1‑149 (MQ=255)                                                                     aGTATAACTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAACTTCTTGATGTGTc                                                                           <  1:687811/149‑1 (MQ=255)                                                                       tataATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAACCAAAAGTAGTTAGAATAATTGGCCGCATGTATGCTGAGAAATTTCTTGATGTGTCtt                                                                         <  1:36025/149‑1 (MQ=255)                                                                       tataATGGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCtt                                                                         <  1:376810/149‑1 (MQ=255)                                                                         taATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCCAAAGTAGTTACAATGATTGGCTGCATGTATGCTGAGAAATTTCTTATTGTGTCTTct                                                                       >  2:503383/1‑148 (MQ=255)                                                                         taATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCa                                                                       <  1:180155/149‑1 (MQ=255)                                                                               cAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTGTTACTTACCTTTGTATCAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCa                                                                 <  1:723538/149‑1 (MQ=255)                                                                                  aGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTATTCAACATCAGCa                                                              >  1:452672/1‑149 (MQ=255)                                                                                      gAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATTTGTCTTCAACATCAGCAtctc                                                          >  1:61518/1‑149 (MQ=255)                                                                                         tACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGAAAAAATAGTTAGAATGATTCGCTGCATTAATGCTAAGAAATTTCTTGATGTGTCTTCAAAAGCATCATCTCttt                                                       >  2:143808/1‑148 (MQ=255)                                                                                          gtaaataaaaCAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTGAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTAt                                                      <  2:625495/143‑1 (MQ=255)                                                                                                           agtaaaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTTTCCATTTAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCCACTTCATATAGcata                                     >  2:409639/1‑146 (MQ=255)                                                                                                           agtaaaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGACAAAGCATACAGCTGGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCCTCTAGCtaa                                     >  1:158839/1‑149 (MQ=255)                                                                                                                   ttttACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTa                             >  1:613994/1‑149 (MQ=255)                                                                                                                     ttaatcacCTTTGTATTAGTAAATTTTTTTGACCAAAGAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGt                           <  2:446201/143‑1 (MQ=255)                                                                                                                           acCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGaa                     <  2:542269/149‑1 (MQ=255)                                                                                                                             cTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGACTGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACt                   >  1:201139/1‑149 (MQ=255)                                                                                                                             cTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACt                   >  2:780368/1‑149 (MQ=255)                                                                                                                               ttGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTgg                 <  1:523822/149‑1 (MQ=255)                                                                                                                               ttGGATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTATAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATTAGCATCTCGTAGCATATTCATATAGCGAGGTAAGCTAGTCCCGAACTgg                 <  2:141664/149‑1 (MQ=255)                                                                                                                                  tATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAg              <  1:542609/149‑1 (MQ=255)                                                                                                                                     tAGTAAATTTATTTGAACAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTa           <  1:806386/149‑1 (MQ=255)                                                                                                                                       gTAAATTTTTTTGATTAAATAAATAAAGAAAAAGTAGTTAGAATGATTTGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGTTACTCCCGAACTGGtgctttct         >  2:223647/1‑142 (MQ=255)                                                                                                                                        tAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATTTAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACtt        >  1:159721/1‑149 (MQ=255)                                                                                                                                        tAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACtt        >  1:383210/1‑149 (MQ=255)                                                                                                                                           aTTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAAc     <  2:819360/149‑1 (MQ=255)                                                                                                                                            ttttttTGACCTAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCACCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACtccacc    <  2:565537/149‑6 (MQ=255)                                                                                                                                              ttttTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTc  <  2:613994/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                 CTCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTAGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTC  >  NC_000913_3_hsa_tpiA/2106596‑2106885

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