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

CTCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTAGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTT  >  NC_000913_3_hsa_tpiA/2106596‑2106879                                                                                                                                                  |                                                                                                                                           ctctATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGAc                                                                                                                                         <  1:554713/149‑1 (MQ=255)      ttAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGATCaaat                                                                                                                                    >  2:490907/1‑149 (MQ=255)              tACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGAAAATTTTTTTGACCAAAtaaataaaa                                                                                                                            >  2:115797/1‑148 (MQ=255)                      aTGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAagta                                                                                                                    <  1:550219/149‑1 (MQ=255)                                 atgatAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGAt                                                                                                         >  1:189567/1‑149 (MQ=255)                                    atAACCAGAAAAACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTgg                                                                                                      <  2:398057/149‑1 (MQ=255)                                          aTAAAGACAATTTGCTTAGCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCCAAAGTAGTTAGAATGATTGGCTGCAt                                                                                                <  2:379664/149‑1 (MQ=255)                                          aTAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCAt                                                                                                >  1:61195/1‑149 (MQ=255)                                           tAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATg                                                                                               >  1:523542/1‑149 (MQ=255)                                                  aaTTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTg                                                                                        >  1:304052/1‑149 (MQ=255)                                                           aCCAGTATTCAGTATAATTGCAAAGATGAATAAAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTc                                                                               <  1:178171/149‑1 (MQ=255)                                                            ccAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTATGATGATTGGCTGCATGTATGCTGAGAAATTTCt                                                                              >  1:388258/1‑149 (MQ=255)                                                             gaGTATTCGGTATAATTTGTAAGATGAATACAATACAAACCAGTAAAGTAACTGTTTTACTTACCTTTGTAGTAGTAAATTTTTTTGTCCTTTCATATAAAGCAAAAGTAGTTAGAATGATTGGCTGTATGTATGCTGAGAAATGTCtt                                                                             <  2:280476/148‑1 (MQ=255)                                                             cAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCtt                                                                             >  1:374963/1‑149 (MQ=255)                                                                                aaaGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAg                                                          >  1:404274/1‑149 (MQ=255)                                                                                      taaTACAATATAAAGCAGTAAAGTAATTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCAtctc                                                    <  2:337404/148‑1 (MQ=255)                                                                                      gAATAAAATAAAAACAAGTAAATTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCAtctc                                                    <  1:332600/149‑1 (MQ=255)                                                                                          actaTAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGGTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTAt                                                <  1:112064/146‑1 (MQ=255)                                                                                             aTAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTATTCAACATCAGCATCTCTTATCaa                                             >  1:356412/1‑149 (MQ=255)                                                                                                aaaaCTAGTAAAGTAAAGTTTCTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATCTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAAtt                                          <  2:40205/149‑1 (MQ=255)                                                                                                        taaaggaaaTGTTTTACTTAACTTTGGAGTAGGAGATTGTTTTGACGAAATAAATAAAGGAAAAGTAGTTAGAATGATTGGCTGCATGTATTCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCGCTTATCAAATTCATATAGc                                  <  2:160280/149‑1 (MQ=255)                                                                                                               aaTGTTTTACTTACCTTTGTATTAGTAAAGTTTTTTGACCAAATAAATAAAGCAAAAGTAGGTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGGCTTCAACATCAGCATCTCTTATCAAATTCATATAGCtaagtaa                           <  2:264022/149‑1 (MQ=255)                                                                                                                        cttaaCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTccc                  <  1:23236/149‑1 (MQ=255)                                                                                                                         ttacCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCg                 >  1:458752/1‑149 (MQ=255)                                                                                                                             ttttGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACt             <  2:473224/148‑1 (MQ=255)                                                                                                                              tttGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAATCTAGTCCCGAACTg            >  1:267485/1‑149 (MQ=255)                                                                                                                              tttGTATTAGTAAATTGTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTg            <  1:207841/149‑1 (MQ=255)                                                                                                                                 gTATTAGTAAATTTTTTTGACCAATTAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTa         <  1:43816/149‑1 (MQ=255)                                                                                                                                 gTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTa         <  2:278681/149‑1 (MQ=255)                                                                                                                                 gTATGAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTa         <  2:361326/149‑1 (MQ=255)                                                                                                                                       gTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACt   <  1:532400/149‑1 (MQ=255)                                                                                                                                        tAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATATCTAAGTAAGCTAGTCCCGAACTGGTAGTTACtt  >  2:373575/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: ‑