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

CTCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTAGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTT  >  NC_000913_3_hsa_tpiA/2106596‑2106879                                                                                                                                                  |                                                                                                                                           ctctATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGAc                                                                                                                                         <  1:219451/149‑1 (MQ=255)          cacCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATACAAACACCGCATTTTGACGGGTTACTTGCCTTTGTATTAGTAAGTTTTTTTTACCaattaaat                                                                                                                                >  1:19893/1‑149 (MQ=255)             cTACGATGAATGCATAAGAAATGATAACCATACAGACAATTTTCTTACCAGTATTCAGTATAATTGCAGGGATAAATAGCATAAATACAAGTAAAGTAAATGTAGTGCTAACCTTTGTATTAGAAAATTTTTTTTACAAAAtaattaaa                                                                                                                             >  2:538228/1‑149 (MQ=255)              tACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAACGTAAATGTTTTACTTACATTTGTCCCAGTAAATTTAGATGAGCAAAtaaataaac                                                                                                                            >  1:630626/1‑148 (MQ=255)                          ctAAGAGATGATAACCATAAATCCAATTTCCTTACCAGTATTCAGTATAATTGGAAATATGAATACAATAAAAACAAGTAAATGAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTa                                                                                                                <  2:128938/148‑1 (MQ=255)                          aTAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGCCTAAATAACTAAACCAAAAGTAATTa                                                                                                                >  2:170656/1‑149 (MQ=255)                           tAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAg                                                                                                               >  1:68400/1‑149 (MQ=255)                                        ccATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGc                                                                                                  >  1:615506/1‑149 (MQ=255)                                               gACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTTTTTTTTTATATAAAGCAAAAGTAGTAAGAATGATTGCCTGCTTGTATg                                                                                           >  2:411851/1‑149 (MQ=255)                                                aCAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGc                                                                                          <  1:487199/149‑1 (MQ=255)                                                   atgggCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAAAGTTTTACTTACCTTTGTATTAATAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTga                                                                                       <  1:191588/145‑1 (MQ=255)                                                         ttACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAAtt                                                                                 >  2:372438/1‑149 (MQ=255)                                                            ccAGGATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATAATTGGCTGCATGTATGCTGAGAAATTTCt                                                                              <  1:574666/149‑1 (MQ=255)                                                                  ttCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGAtgt                                                                        <  1:449112/149‑1 (MQ=255)                                                                   tCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGAtgtg                                                                       >  1:607702/1‑149 (MQ=255)                                                                           aTTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATATATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAAc                                                               <  1:500952/149‑1 (MQ=255)                                                                               cAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCa                                                           <  1:92357/149‑1 (MQ=255)                                                                                          aCAATAAAAACAAGTAAAATAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAATCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTAt                                                <  2:473152/149‑1 (MQ=255)                                                                                                           agtaaaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCtaa                               <  2:68400/149‑1 (MQ=255)                                                                                                              aaactttttACTTACCCTTGTATTAGTCAATTTTTTTGACAAAACAAATAACGCAGAAGTAGTTAGAATGATTGGCTGCATGGATGCTGAGAAATTTCTTGATGTGTCGTCAACATCAGCATCTCTTATCAAATTCATATAGCtgagta                            <  2:348086/144‑1 (MQ=255)                                                                                                                 tGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGc                         <  1:372438/149‑1 (MQ=255)                                                                                                                   ttttACTTACCTTTGTAATAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATTGTCAAGTAAGCTa                       >  2:233253/1‑149 (MQ=255)                                                                                                                             cTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAGGTAGTTAGAATGATTGGCTGAATGTATGCTGACAAAGTTCTTGATGTGTCTTCAAAATCAGCATCTCTTATCAAATTCATATAGCTTATTAAGCTAGTCCTGAACt             >  2:450473/1‑149 (MQ=255)                                                                                                                             cTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACt             >  1:584300/1‑149 (MQ=255)                                                                                                                                        tAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGAtgt                                                                        >  1:482634/1‑79 (MQ=255)                                                                                                                                        tAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGAtgt                                                                        <  2:482634/79‑1 (MQ=255)                                                                                                                                        tAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACtt  >  2:143343/1‑149 (MQ=255)                                                                                                                                        tAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGGCTAATTATTAAATTCTTGATTTGTCTTCAATCTCTTAACCGCGGATCAAATGCATATAGCTAAGTAAGCTTATTCCGAACTGGTAGTTACtt  >  2:651882/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                           CTCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTAGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTT  >  NC_000913_3_hsa_tpiA/2106596‑2106879

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