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_tpiA2,106,7410AT96.0% 79.6 / NA 25S162S (TCT→TCAwbbHO‑antigen polymerase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (6/18);  minor base G (1/0);  total (7/18)
Fisher's exact test for biased strand distribution p-value = 2.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTAGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACC  >  NC_000913_3_hsa_tpiA/2106596‑2106883
                                                                                                                                                 |                                                                                                                                              
ctctATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTGTACTTACCTTTGTATTAGTAAATTTTTTTGAc                                                                                                                                             <  1:505007/149‑1 (MQ=255)
 tctATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGAcc                                                                                                                                            >  2:536557/1‑149 (MQ=255)
 tccaTTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGAcc                                                                                                                                            <  2:486152/146‑1 (MQ=255)
                 aTGAATGCATCAGAGCTGATAACCATAAAGACAATTTCCTTACCAGGATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAAGATTTTGTACAAATTCAATCAAGCaaa                                                                                                                            >  2:70537/1‑149 (MQ=255)
                                     aaCCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACACTAAAAACAAGTAAAGTAACTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATATAGCAGACATAGTGACAATGTATGGCt                                                                                                        >  2:190249/1‑149 (MQ=255)
                                                             aGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTg                                                                                <  1:424861/149‑1 (MQ=255)
                                                                 ttCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGAtgt                                                                            <  1:565830/149‑1 (MQ=255)
                                                                       ataATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTc                                                                      <  1:500467/149‑1 (MQ=255)
                                                                         aaTTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCaa                                                                    <  1:310125/149‑1 (MQ=255)
                                                                           ttGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACa                                                                  <  1:99939/149‑1 (MQ=255)
                                                                               aaaGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAg                                                              >  1:219679/1‑149 (MQ=255)
                                                                                     gAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCAtctc                                                        >  2:25741/1‑149 (MQ=255)
                                                                                        tACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCaaataaa                                                                                                                                     <  2:409660/69‑1 (MQ=255)
                                                                                        tACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCaaataaa                                                                                                                                     >  1:409660/1‑69 (MQ=255)
                                                                                                 aaCAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCa                                            <  2:139041/149‑1 (MQ=255)
                                                                                                 aaCAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCa                                            <  2:139325/149‑1 (MQ=255)
                                                                                                  aCAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTTTTTTTTTAAATAAAGCAAAAGTAGTTAGAATGATTGCCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCATCATAAGCATCTCTTATCAATTTCat                                           >  2:496888/1‑149 (MQ=255)
                                                                                                   cAAGTAAAGTAAATGTTTTACTTACGTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTCGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCata                                          <  2:4623/149‑1 (MQ=255)
                                                                                                    aagtaaagtaAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCatat                                         <  1:419922/149‑1 (MQ=255)
                                                                                                     agtaaagtaaATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCAtata                                        <  1:207763/149‑1 (MQ=255)
                                                                                                          agtaaaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCtaa                                   <  1:536557/149‑1 (MQ=255)
                                                                                                                tGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGc                             <  1:25741/149‑1 (MQ=255)
                                                                                                                              ttGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTgg               <  1:496888/149‑1 (MQ=255)
                                                                                                                                    tAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTa         <  2:390652/149‑1 (MQ=255)
                                                                                                                                           ttttttTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAAcc  <  2:78088/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
CTCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTAGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACC  >  NC_000913_3_hsa_tpiA/2106596‑2106883

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: