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	97.8%	150.5 / NA	46	S162S (TCT→TCA)	wbbH	O‑antigen polymerase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (24/21);  minor base C (0/1);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTAGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAA  >  NC_000913_3_hsa_tpiA/2106600‑2106881                                                                                                                                              |                                                                                                                                             aTTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAATTTTTTTTGACCaaa                                                                                                                                       >  1:130821/1‑149 (MQ=255) aTTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAATTTTTTTTGACCaaa                                                                                                                                       >  1:130899/1‑149 (MQ=255)     acacCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGAc                                                                                                                                           <  1:950083/141‑1 (MQ=255)     acacCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGAc                                                                                                                                           >  2:950083/1‑141 (MQ=255)       gccTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAACGATGAATACAATAGAAACAAGTAAAGTAAATGTTTTACTGACCTTTGTATTAGTAAATTTATTTGACCAaataaata                                                                                                                                 >  2:968323/2‑149 (MQ=255)         cTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAAtaaataaa                                                                                                                               >  2:636269/1‑149 (MQ=255)          tACGATGAATGCATAAGAGATGATAACCATGAAGACAATTTGCTTACCAGTATTCAGTATAATTGCACAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAGTTTTTTTGAACAAATAGATAAAg                                                                                                                              >  1:420487/1‑149 (MQ=255)                      aTAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAATTTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTa                                                                                                                  <  2:130899/149‑1 (MQ=255)                      aTAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAATTTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTa                                                                                                                  <  2:130821/149‑1 (MQ=255)                         aaagaTGAGAACCATAAAAACAATTTGCTTACCAGTATTAAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGaa                                                                                                               <  2:243014/147‑1 (MQ=255)                             atgatAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGAt                                                                                                           <  1:46575/149‑1 (MQ=255)                                 tAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTAACTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGc                                                                                                       <  1:987688/149‑1 (MQ=255)                                     cATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTTCa                                                                                                   >  1:210568/1‑149 (MQ=255)                                          aGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCACAAGTAGTTAGAATGATTGGCTGCGTGTAt                                                                                              >  1:862366/1‑149 (MQ=255)                                                tttGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGAAAGAGTAGTTAAAATGAATGGCTGCATGTTTGCTgag                                                                                        >  1:939773/1‑149 (MQ=255)                                                         cAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCtt                                                                               <  1:636269/149‑1 (MQ=255)                                                          aGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTg                                                                              >  2:769195/1‑149 (MQ=255)                                                           gTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGa                                                                             >  1:592197/1‑149 (MQ=255)                                                              ttCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGAtgt                                                                          <  1:223057/149‑1 (MQ=255)                                                              ttCAGTATAATTGCAAAGATGAATAAAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGAtgt                                                                          <  2:508417/149‑1 (MQ=255)                                                                 aGTATAATTGCAAAGATGAATACAATAAAAACAAGTATAGTAAATGTTTTACTTAACTTTGTATTAGTAAATGTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTc                                                                       <  2:247834/149‑1 (MQ=255)                                                                    ataATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGACTGCATGTATGCTGAGTAATTTCTTGATGTGTCTTc                                                                    >  1:280325/1‑149 (MQ=255)                                                                    agaaTTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTc                                                                    <  2:787708/147‑1 (MQ=255)                                                                      aaTTGTAAAGAAGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTTGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCaa                                                                  <  2:465420/149‑1 (MQ=255)                                                                         tGCAAATATGAATACAATAAAAACATGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATGTTTTTGACCAAATAAACCAAGCAAAAGTAGTTAGAATCATTGGCTGCATGTATGCTGACAAATTTCTTGATGTGTCTTCAACAt                                                               <  1:948577/149‑1 (MQ=255)                                                                           cAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCa                                                             >  1:639388/1‑149 (MQ=255)                                                                                      aCAATGAAAACGAGTGAAGGAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAACTAAAGCAAAAGTAGTTAGAATGATTGGCTTCAGGTAGGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTAt                                                  <  1:101518/149‑1 (MQ=255)                                                                                      aCAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAAttt                                                                                  <  2:226918/117‑1 (MQ=255)                                                                                      aCAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAAttt                                                                                  >  1:226918/1‑117 (MQ=255)                                                                                        aaTAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCa                                                >  1:528116/1‑149 (MQ=255)                                                                                         aTAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGa                                                                             >  2:482454/1‑119 (MQ=255)                                                                                         aTAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGa                                                                             <  1:482454/119‑1 (MQ=255)                                                                                         aTAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCaa                                               >  1:723930/1‑149 (MQ=255)                                                                                          tAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCaaa                                              >  1:856034/1‑149 (MQ=255)                                                                                          tAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCaaa                                              >  1:11691/1‑149 (MQ=255)                                                                                             aaaCACTTAAAGTAAACGTTTTACTTACCTTTGTATTAGTAAATTTGTTTGACCAAATAAATAAAGCAAAAGGAGTTAGCATGACTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTc                                           <  1:323701/149‑1 (MQ=255)                                                                                               aCAAGTAAAGTAAATGTTTTACTTAGCTTTGTATTAGTCACTTTTTTCGACCAAATAAATAAAGCACAAGAAATGAGAATGCTTTCTGGCAGTTGTGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCat                                         <  1:1031184/149‑1 (MQ=255)                                                                                                    taaagtaaaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGc                                    >  2:181758/1‑149 (MQ=255)                                                                                                     aaagtaaaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTTTCTAATTAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATTTTTCTTCCACATTAGCATCTCTTTTTAAATTTATATAGCt                                   >  2:344159/1‑149 (MQ=255)                                                                                                       tgtgagtGTTTTACTCACGTTGGTATTAGTAAATTTCTCTGACCAAATAAATAAATCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATGAAATTCATATAGCtaa                                 <  2:644866/143‑1 (MQ=255)                                                                                                        gtaaaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCACCATCAGCATCCCTTATCAATTTCATTTAGAtaat                                >  2:934197/1‑148 (MQ=255)                                                                                                                     ttacCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTTTCAAATTCATATAGCTAAGTAAGCTAGTCCCg                   >  1:361004/1‑149 (MQ=255)                                                                                                                         cTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACt               <  2:528116/149‑1 (MQ=255)                                                                                                                         cTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAAAGTTTGGCGGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGCACt               >  1:300722/1‑149 (MQ=255)                                                                                                                          tttGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTg              <  2:319106/149‑1 (MQ=255)                                                                                                                                      aaTTTTTTTGACCAAATAAATAAAGCAAAAGGAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTaa  <  2:280325/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                             ATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTAGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAA  >  NC_000913_3_hsa_tpiA/2106600‑2106881

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