Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_hsa_tpiA 2,106,741 A→T 100% 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,7410AT98.4% 94.0 / ‑6.4 33S162S (TCT→TCAwbbHO‑antigen polymerase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base T (16/15);  total (16/17)
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
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTAGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTT  >  NC_000913_3_hsa_tpiA/2106597‑2106879
                                                                                                                                                |                                                                                                                                          
tctATTAACACCTACGATGAATGCATAAGAGATGATAACCATCAAGACAATTTGCTTACCAATATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGCATTAGTAAATTTTTTTGAcc                                                                                                                                        >  1:161490/1‑149 (MQ=255)
tctATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGAcc                                                                                                                                        >  1:591204/1‑149 (MQ=255)
tctATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACGTACCTTTGATTTAGTAAATTTTTTTGAcc                                                                                                                                        >  1:325723/1‑149 (MQ=255)
            tACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATAAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAACTGAAg                                                                                                                            >  2:453531/1‑149 (MQ=255)
                 tGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCaaaa                                                                                                                       <  1:558990/149‑1 (MQ=255)
                  gAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAag                                                                                                                      >  1:143476/1‑149 (MQ=255)
                    agcaaTAAGAAATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAATTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAagta                                                                                                                    <  2:423479/145‑1 (MQ=255)
                      tcagaaGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAagtagt                                                                                                                  <  2:322483/145‑1 (MQ=255)
                       cATAAGAGATGATAACCATAAAGACAATTTGCTAACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAACACAAGTAAAGTAAATGTTTTACTTACCATTGTAATATTAAATTTTTTTCAACCAAAAAATAAAGTAAAagtagtg                                                                                                                 >  2:70680/1‑148 (MQ=255)
                                  atAACCATAAAGACATTTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAGTAATAACGAGTCTAGTGTCTGTTTTAGGTGGCTTTGTATTAGTAAATTTTTTTGATCAATTAATTTAAGTAAAAGTCGTTAGATTGATTgg                                                                                                      >  2:237098/1‑149 (MQ=255)
                                             ggaaaattGCTTACCCGTATCCAGAATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTCGAATGATAGGATGCATGTATg                                                                                           >  2:208453/7‑149 (MQ=255)
                                                    tgtttACCAGTAATAAGTATAATTACAAAGATGAATACAAGAAAAACAAGTAAAGTAAATTTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGaa                                                                                    <  1:386086/146‑1 (MQ=255)
                                                     gCTTACCATTTCCCAGTATAACTGTAAAGATAAATAAAATGACTACAAGTAAAGTAAATGTTTTTCTTACCTTTGTATTAGTAAATTTTTTGACCAAATAACTAACCCAAAAGTAGTAAGAATGTTTCGATGCATGTATGCTgagttat                                                                                  <  2:382005/149‑5 (MQ=255)
                                                           cAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCtt                                                                             >  1:12385/1‑149 (MQ=255)
                                                             gTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGa                                                                           <  2:377930/149‑1 (MQ=255)
                                                                     tataATTGCAAAGATGAAGAAAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCtt                                                                   <  1:436462/149‑1 (MQ=255)
                                                                             cAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCa                                                           <  2:289654/149‑1 (MQ=255)
                                                                                 gATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCAt                                                       <  1:302085/149‑1 (MQ=255)
                                                                                   tGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCGGCAtct                                                     >  1:61514/1‑149 (MQ=255)
                                                                                   tGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCAtct                                                     <  1:367499/149‑1 (MQ=255)
                                                                                        aaaaTAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTAt                                                <  1:607841/147‑1 (MQ=255)
                                                                                         ctaTAAAAACAACCAACGTAAATTTTTAACTGACTTTTGTATTAGGAAATTTTTTTGACTAAATAAACAAAGCAAAAGTATGTAGAATGATTGTCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATc                                               <  2:358227/147‑1 (MQ=255)
                                                                                           aTAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCaa                                             >  1:477438/1‑149 (MQ=255)
                                                                                               aaaCAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTc                                         >  1:362159/1‑149 (MQ=255)
                                                                                                      taaagtaaaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTATTCACCATCAGCGTCTCTTATCAAATTCAtatatc                                  >  2:377911/1‑147 (MQ=255)
                                                                                                        aagtaaaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTTTTCAACATCAGCATCACTTATCAAATTCATATAGCtt                                >  2:182230/1‑148 (MQ=255)
                                                                                                              aTGTTTTACTTACATTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCtaagtaag                          <  1:581545/149‑1 (MQ=255)
                                                                                                                 ttttACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTa                       <  2:534413/149‑1 (MQ=255)
                                                                                                                           ccttGTATTAGTAAATTTTTTTCCATACATGAATAAAGAACCAGCTTTTCCAAGTATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACt             <  1:312699/147‑1 (MQ=255)
                                                                                                                           cTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACt             <  1:170627/149‑1 (MQ=255)
                                                                                                                              tGTATTAGTAAATTTTTTAGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGCGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGt          <  1:641754/149‑1 (MQ=255)
                                                                                                                                tATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAg        >  2:258959/1‑149 (MQ=255)
                                                                                                                                      tAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACtt  >  1:412707/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                          
TCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTAGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTT  >  NC_000913_3_hsa_tpiA/2106597‑2106879

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