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

TCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTAGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACC  >  NC_000913_3_hsa_tpiA/2106597‑2106883                                                                                                                                                 |                                                                                                                                               tctATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGAcc                                                                                                                                            >  2:72272/1‑149 (MQ=255)                       gCATAAGAGAACATAACAATAAAGAAGATGTCCCTACCAGTATTTATGATAATTGCAAAGATGAATACAATAAAAACCAGTATAGTAAATTTTTTACTTACCGTTGTATTAGTAAATTTCTGTGATCAAATCAATAAAGCAAAagtagt                                                                                                                      <  2:315192/149‑1 (MQ=255)                              atatgatAACCATAAAGACAATTTGCGTACAAGAATTCAGTATAATTGCAAAGATAAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATg                                                                                                               <  2:68130/147‑1 (MQ=255)                                  gatAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTg                                                                                                           <  2:460194/149‑1 (MQ=255)                                    tAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATATATAAAGCAAAAGTAGTTAGAATGATTGGc                                                                                                         >  1:576847/1‑149 (MQ=255)                                     aaCCATAAATACAATTTGCTTACCATGATTGAGTATAATTGCAAAGATGAATAAAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCt                                                                                                        <  2:186387/149‑1 (MQ=255)                                            aaaaCAATTTGCTTACCAGTATTTAGTATAATTGAAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTa                                                                                                 <  2:168189/146‑1 (MQ=255)                                            aaGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAATAAAAATTAATTAGAATGATTGGCTGCATTTa                                                                                                 >  2:179662/1‑149 (MQ=255)                                             aGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAGAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTAt                                                                                                >  2:638294/1‑149 (MQ=255)                                                    ttGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTgaga                                                                                         <  1:489116/149‑1 (MQ=255)                                                      gCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGaaa                                                                                       >  1:636125/1‑149 (MQ=255)                                                       cTTACAAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAAt                                                                                      <  2:42231/149‑1 (MQ=255)                                                          aCCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTc                                                                                   <  2:636125/149‑1 (MQ=255)                                                               tATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTTCATGTATGCTGAGAAATTTCTTGAt                                                                              >  1:413719/1‑149 (MQ=255)                                                                       ataATTCCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAATTAAATAAAGCAAAAGTAGTTGGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTc                                                                      <  1:606808/149‑1 (MQ=255)                                                                                 aGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCa                                                            >  1:408961/1‑149 (MQ=255)                                                                                      aaTACAATAAAAACAAGTAGAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATctct                                                       <  1:638294/149‑1 (MQ=255)                                                                                        tACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTa                                                     <  2:63415/149‑1 (MQ=255)                                                                                         aCAATAAAAACAAGTAAAGTAAATGGTTTACTTACCTTTGTATTCGCCATTTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTAt                                                    <  1:153893/149‑1 (MQ=255)                                                                                            aTAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTCTTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCaa                                                 >  1:600149/1‑149 (MQ=255)                                                                                            aTAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCaa                                                 <  2:408961/149‑1 (MQ=255)                                                                                             tAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCaaa                                                <  2:8454/149‑1 (MQ=255)                                                                                             tAAAAACAAATAAAGTAAAGATTTAACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAACCTATATCGCTGTCTTCAACATCAGCATCTCTTATCaaa                                                <  2:482499/149‑1 (MQ=255)                                                                                                    aagtaaagttAATGTTTTACTTAACTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCACCATCTCTTATCAAATTCatat                                         <  2:591336/149‑1 (MQ=255)                                                                                                          agtaaaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCtaa                                   <  2:600149/149‑1 (MQ=255)                                                                                                          agtaaaTGTTTTACTTACCTTGTTATTAATAACTTTTTTTGACCAGAAAAGAAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCtaa                                   <  2:71935/149‑1 (MQ=255)                                                                                                           gtaaaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGATCTATTAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTTAGAAATTTCTTCATGTGTCTTCAACATCAGCATCTCTTATCACATTCAAATATCtact                                  >  2:325866/1‑147 (MQ=255)                                                                                                            taaaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCtaatt                                 <  1:34602/149‑3 (MQ=255)                                                                                                              aaTGTTTTACTTACCTGTGTAGTAGAAAATATTTTTGACCAAATGAATAAAGCAAAAGTAGTTAAAATGATTGGCTGCATGTGTGCTGAAAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCttagtca                               >  2:518767/1‑147 (MQ=255)                                                                                                                tGTTTTACTAAACTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGc                             <  2:15226/149‑1 (MQ=255)                                                                                                                            cTTTGTATTAGTAAATTTTTTTTTTTTATTAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAAATTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATTTATCTAATTAAGCTAGTCCCGAACt                 >  2:552797/1‑149 (MQ=255)                                                                                                                              ttGTATTAGTAAATTTTTTTGACCAAATACATAAAGAAAACGGAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTgg               <  1:461621/149‑1 (MQ=255)                                                                                                                                  agtAGTAAATTTTTTTGAACAGATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGt           <  2:656059/147‑1 (MQ=255)                                                                                                                                           ttttttTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTCTGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAccc  >  1:151842/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                               TCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTAGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACC  >  NC_000913_3_hsa_tpiA/2106597‑2106883

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