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	92.0%	68.4 / ‑4.6	25	S162S (TCT→TCA)	wbbH	O‑antigen polymerase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (13/10);  minor base G (1/1);  total (14/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTAGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTA  >  NC_000913_3_hsa_tpiA/2106595‑2106858                                                                                                                                                   |                                                                                                                      aCTCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGa                                                                                                                     >  1:151932/1‑149 (MQ=255)     tATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCaa                                                                                                                 >  1:33915/1‑149 (MQ=255)      aTTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCaaa                                                                                                                >  2:360856/1‑149 (MQ=255)            acCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATAAGTAAATTTTTTTGACCAaataacta                                                                                                          >  2:224508/1‑149 (MQ=255)                  gATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCaa                                                                                                    >  1:314561/1‑149 (MQ=255)                  gATGAATGCATAAGAGAATATAAACATACAGACAATTTGCTTACCAGTATTCAATATAATTGCAAAGATGAATACAATAAAAACAAGTAAATTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCaa                                                                                                    <  1:362743/149‑1 (MQ=255)                    agAATGAAAAAGAGATGATAAACATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCaaaa                                                                                                  <  2:33915/148‑1 (MQ=255)                           aTAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAATATAATTGCAAAGATGAATACCATAAACACAAGTAAAGTAAATGTTTTACTTACCTTTGTAGTAGTAAATTGTTTTGACCAGATAAATAAAGCAAAAGTAGTTa                                                                                           >  2:37319/1‑149 (MQ=255)                               gagaTGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAAt                                                                                       >  1:226085/1‑149 (MQ=255)                                               acaCAACTGGCTTATCAGCATGCATGATAAATACAAAGACGAATAAATTAAAAACAAGTAAAGTAAATGTTTTACTTACGTTTGTATTAGTAAATTTTTTTGACCCAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTAt                                                                       <  1:321130/147‑1 (MQ=255)                                                          ttACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAACATAGGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAAtt                                                            >  1:289561/1‑149 (MQ=255)                                                            aCCAGTATTCAGTATAATTGCAAAGATGAAAACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTc                                                          <  2:135127/149‑1 (MQ=255)                                                              cAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCtt                                                        >  1:147147/1‑149 (MQ=255)                                                                 tATTTAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGAt                                                     <  2:42522/149‑1 (MQ=255)                                                                                   aGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCa                                   <  1:282107/149‑1 (MQ=255)                                                                                         aTACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTGCGTTTGATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCtt                             <  1:97697/149‑1 (MQ=255)                                                                                          tACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAATTAGTTAGAATGATTGACTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTa                            >  2:340208/1‑149 (MQ=255)                                                                                          tACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTa                            >  1:211834/1‑149 (MQ=255)                                                                                              aTAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGGAATTTGTTGATGTGTCTTCAACATCAGCTTCTCTTATCaa                        >  2:234005/1‑149 (MQ=255)                                                                                               tAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATGTTGTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCaaa                       <  1:78284/149‑1 (MQ=255)                                                                                                       agtaaagtaaATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCAtata               <  1:252319/149‑1 (MQ=255)                                                                                                        gtaaagtaaaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAg              >  2:329780/1‑149 (MQ=255)                                                                                                            agtaatTGTCTTACTTACCTTTGTATTAGTCAATTTTTTTGACCAAAGAAATAAAGCAAAAGTAGTTAGAATGATTGGATGCATGTATGCTGAGATATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATTAAATTCATATAGCtaa          <  1:344511/149‑1 (MQ=255)                                                                                                                 aTGTTTTACTTATCTTTGTATTAGTAAATTTTTTTGACCATATAATTAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGATATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCtaagtaag     <  2:203791/149‑1 (MQ=255)                                                                                                                    ttttACTTACCTTTGTATTAGTAAAATTATTGCGCCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTa  >  1:101726/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                      ACTCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTAGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTA  >  NC_000913_3_hsa_tpiA/2106595‑2106858

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