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

ACTCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAA‑TTTTTTAGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTT  >  NC_000913_3_hsa_tpiA/2106595‑2106879                                                                                                                                                    |                                                                                                                                           aCTCTATTAAACACTACGCGGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAA‑TTTTTTTGAc                                                                                                                                         <  2:319656/150‑1 (MQ=255)      aTTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAA‑TTTTTTTGACCaaat                                                                                                                                    <  2:162567/150‑1 (MQ=255)           aaacTACGATGAATGAATAAGAGAGGATAACCATAAAGACAATTTGCTTACCAGTATTCACTCTAATGTGAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAA‑TTTTTTTGACCAaataaata                                                                                                                               <  1:323495/147‑1 (MQ=255)            acCTACGATGAATGCAGAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAA‑TTTTTTTGACCAAataaataa                                                                                                                              <  2:595244/150‑1 (MQ=255)                   aTGAATGCATAATAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAA‑TTTTTTTGACCAAATAAATAAAGCaaaa                                                                                                                       <  1:556299/150‑1 (MQ=255)                     gAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTTAGTATTATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTGACTTACCTTTGTATTGGTAAA‑TTTTTTTGACCAAATAAATAAAGCAAAagt                                                                                                                     <  2:536115/150‑1 (MQ=255)                       aTGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAA‑TTTTTTTGACCAAATAAATAAAGCAAAagtag                                                                                                                   >  2:532327/1‑150 (MQ=255)                       aTGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAA‑TTTTTTTGACCAAATAAATAAAGCAAAagtag                                                                                                                   <  1:231655/150‑1 (MQ=255)                            tAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAA‑TTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGa                                                                                                              >  1:385493/1‑150 (MQ=255)                             aaGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCGTTGTATTAGTAAA‑TTTTTTTGACCAAATAAATAAAGGAGAAGTAGTTACaa                                                                                                             >  1:479803/1‑150 (MQ=255)                                  atgatAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAA‑TTTTTTTGa                                                                                                                                          <  1:553776/116‑1 (MQ=255)                                      aaaCCCTAAAGACAATTTGCCTAACAGTATTTAGTATAATTGCAAAGATGAATATAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTTTATTAGTAAATTTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGc                                                                                                     <  2:18571/149‑1 (MQ=255)                                                   aaTTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAA‑TTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTga                                                                                       >  2:279180/1‑150 (MQ=255)                                                              cAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAA‑TTTTTTTGACCaaataaa                                                                                                                                 >  2:118343/1‑97 (MQ=255)                                                              cAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAA‑TTTTTTTGACCaaataaa                                                                                                                                 <  1:118343/97‑1 (MQ=255)                                                                  aTTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAA‑TTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGAtgt                                                                        <  2:385493/150‑1 (MQ=255)                                                                  aTTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAA‑TTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGAtgt                                                                        <  1:532327/150‑1 (MQ=255)                                                                       gTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAA‑TTTTTTTGACAAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCtt                                                                   >  1:155078/1‑150 (MQ=255)                                                                          taATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAA‑TTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCaa                                                                <  2:408698/150‑1 (MQ=255)                                                                               gCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAA‑TTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCa                                                           <  2:155078/150‑1 (MQ=255)                                                                                          tACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAA‑TTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTAt                                                >  1:317858/1‑150 (MQ=255)                                                                                             aaTAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAA‑TTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCt                                                                                         >  1:126918/1‑106 (MQ=255)                                                                                             aaTAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAA‑TTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCt                                                                                         <  2:126918/106‑1 (MQ=255)                                                                                              aTAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAA‑TTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCaaa                                            <  1:279180/150‑1 (MQ=255)                                                                                                        gtaaagtaaaTGTTTTACTTACCTTTGTATTAGTAAA‑TTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGc                                  <  2:613684/150‑1 (MQ=255)                                                                                                             gtaaaTGTTTTACTTACCTTTGTATTAGTAAA‑TTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCtaagt                             >  1:132040/1‑150 (MQ=255)                                                                                                              taaaTGTTTTACTTACCTTTGTATTAGTAAA‑TTTTTTTGACCAAATAAGTAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCtaagta                            >  1:215194/1‑150 (MQ=255)                                                                                                                    ttttACTTACCTTTGTATTAGTAAA‑TTTTTTTGTCTAAATAAATAAAGCAAAAGTAGTTAGAAAGATTGACTGCATGTATGCTGAGAAATTTCTTGATGTGTCATCAACATCAGCATCTCTTATCAAATTAAAATATCTAACTAAGCTat                      >  2:229032/1‑149 (MQ=255)                                                                                                                       tacttacCTTTGTATTAGTAAA‑TTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGACTGGCTGCATTGATGCTGATAAATTGCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTcc                   <  2:406695/150‑1 (MQ=255)                                                                                                                       tacttacCGTTATATTAGTAAA‑TTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTATAAACATTGGCGGCATTTATGGTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTcc                   <  2:406756/150‑1 (MQ=255)                                                                                                                              cTTTGTATTAGTAAA‑TTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTg            >  1:179297/1‑150 (MQ=255)                                                                                                                              cTTTGTATTAGTAAA‑TTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTg            >  1:166366/1‑150 (MQ=255)                                                                                                                                        gTAAA‑TTGTTTTGACCAAGTAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACtt  <  2:596979/150‑1 (MQ=255)                                                                                                                                           aa‑TTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCg                 <  2:268217/132‑1 (MQ=255)                                                                                                                                           aa‑TTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCg                 >  1:268217/1‑132 (MQ=255)                                                                                                                                                    |                                                                                                                                           ACTCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAA‑TTTTTTAGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTT  >  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: ‑