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,7410AT100.0% 133.0 / ‑6.7 44S162S (TCT→TCAwbbHO‑antigen polymerase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (22/21);  total (23/21)
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.

TCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTAGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAA  >  NC_000913_3_hsa_tpiA/2106597‑2106881
                                                                                                                                                 |                                                                                                                                             
tctcTTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGGAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTGTTTTAGag                                                                                                                                            >  2:311527/1‑148 (MQ=255)
   aTTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCaaa                                                                                                                                        <  2:391792/149‑1 (MQ=255)
    ttAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCACTATTCAGTATAATTGCAAATATGAATACAATAAAAACAAGTAAATTAAATGTTTTACTTACTTTTGTATTAGTAAATTTTTTTGACCaaat                                                                                                                                       >  2:455148/1‑149 (MQ=255)
           aatctaTGAAAACATAAGAGAAGACAACCAAAAAGAATACTTTTTTACCACTATTCAATATAATTAAAAAGATGAATACAATAAAATCAAGTAAAGTAAAGGTTTTACTTACCTTTGTATTAGTAAATTTTTTTTACCAAAtaaataaa                                                                                                                                <  2:197727/144‑1 (MQ=255)
                  gAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAag                                                                                                                         <  1:527552/149‑1 (MQ=255)
                    aTGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAagta                                                                                                                       <  2:491087/149‑1 (MQ=255)
                         tAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAg                                                                                                                  >  2:593619/1‑149 (MQ=255)
                               atgatAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGAt                                                                                                            <  1:668184/149‑1 (MQ=255)
                                  atAACCATAAGTACAATTTGCTTACCATTATTCAGTATAAGTGCAAAGATGAATACAATAAAAACAAGTAAAGGAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTgg                                                                                                         <  2:581768/149‑1 (MQ=255)
                                  atAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAATGATGATTACAATAAAAACAATCAAAGTAAATGTTTTACTTATCTTTGTATTAGTAAATTTTTTTGACCAAATCAATATAGCAAAAGTAGTAAGACTGATTgg                                                                                                         >  2:198096/1‑149 (MQ=255)
                                      ccATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAGAAACAAGTAAAGTAAATGTTTTACTTTCCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAATCAGAATTAGTTAGAATGTTTGTCTga                                                                                                     >  2:618951/1‑148 (MQ=255)
                                       cATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCa                                                                                                    >  1:339052/1‑149 (MQ=255)
                                         tAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAATTAAATCAACCAAAACTAGTTAGAATGATTGCCTGCATg                                                                                                  >  2:366885/1‑149 (MQ=255)
                                            aGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTAt                                                                                               >  1:457758/1‑149 (MQ=255)
                                             gACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATg                                                                                              >  2:559786/1‑149 (MQ=255)
                                              aCAAATTTCTTAACAGTATTAAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTATACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGc                                                                                             <  2:419947/149‑1 (MQ=255)
                                                 atctGCCTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTga                                                                                          <  1:216930/146‑1 (MQ=255)
                                                  tttGCTTACCATGATTCTATATAATTGCAAAGATTAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTgag                                                                                         <  1:40344/149‑1 (MQ=255)
                                                          ccAGTATTCAGTATAATTGCAAAGATGACTACTTTAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCt                                                                                 >  2:147375/1‑149 (MQ=255)
                                                           cAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCtt                                                                                >  1:554767/1‑149 (MQ=255)
                                                           cAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCtt                                                                                <  1:289053/149‑1 (MQ=255)
                                                              tatgtcatataATAGAAAAGATGAATACAATAAAAACAAGTAAAGTAAATGATTTACTTACCTTTGTATTAGTAAATTTTTTTCACCAAATCAATAAAGCAAAAGTATTGAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGAt                                                                             <  2:664793/142‑1 (MQ=255)
                                                                   aGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTc                                                                        >  1:106820/1‑149 (MQ=255)
                                                                        aaacacAAAAAAGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCaa                                                                  <  1:237058/144‑1 (MQ=255)
                                                                           tGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACAt                                                               <  2:373883/149‑1 (MQ=255)
                                                                           tGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGGGTCTTCAACAt                                                               >  1:139975/1‑149 (MQ=255)
                                                                                aGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTCAAATAAGTAAATTTTATTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCa                                                          >  1:237073/1‑149 (MQ=255)
                                                                                           aTAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCATCATCTCTTATCaa                                               >  2:649347/1‑149 (MQ=255)
                                                                                            tAAAAACAAGTAAAGTAAAGGTTTTACTTAGCTTTATATTAGTAACTTTGTTTGACCGAATAAAGAAAGCAAAAGTAGTTAGTATGATTGGCTGCATGTATGGTGAGAGATTTCTTGAGGTGTCTTCAACATCAGCATCTCTTATCaaa                                              <  2:304196/149‑1 (MQ=255)
                                                                                             aaaacCAATCAAAGTAAATGTTCAACTCACGATTGTATTAGTAAATTTTTTTTACCAAATAAATAAATCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTATAACACTTTTTACATCAGCATCTCTTATCAAAt                                             <  2:667958/149‑1 (MQ=255)
                                                                                             aaaaaCAGCGAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTACAATGATTGGCTCCACGTATGCTGAGAAATCTCCTGATCTGCCTCCAACATCACCATCTCTTATCAAAt                                             <  2:330460/149‑1 (MQ=255)
                                                                                              aaaaCAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGAAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAAtt                                            >  1:281985/1‑149 (MQ=255)
                                                                                              aaaaCAAGTAAAGTAAATGTTTTACATACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAAtt                                            <  1:426494/149‑1 (MQ=255)
                                                                                                 aCAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCat                                         <  2:690716/149‑1 (MQ=255)
                                                                                                      taaagtaaaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGc                                    <  1:102137/149‑1 (MQ=255)
                                                                                                       aaagtaaaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAATCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCt                                   <  1:655020/149‑1 (MQ=255)
                                                                                                         agtaaaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTTAGAAATTTTTAGTTGTAACGACAAGATCAGCATCTCTTATCAAATTCATATAGCtaa                                 >  1:645367/1‑149 (MQ=255)
                                                                                                             aaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCtaagtaa                             >  1:258791/1‑149 (MQ=255)
                                                                                                                       ttacCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTACCg                   >  1:17427/1‑149 (MQ=255)
                                                                                                                            tttGTATTAGTAAATTTTTTTTTTTTATTTTTTAAAGCAAAAATAGTTAGAATGATTGGCTGCATGTATGCTGGGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCTTTTATCTATTTAATCTAGTCCAGAACtt              >  2:266597/1‑148 (MQ=255)
                                                                                                                               gTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTgcga           >  2:108129/1‑146 (MQ=255)
                                                                                                                                tATTAGTAAATTTTTTTGAACAATTAAGTAAAGCAAAAGGAGATAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATATCTAATTAACCTAGTCCCGAACTgttcg          >  2:81812/1‑145 (MQ=255)
                                                                                                                                  ttAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTAGATGTGTCTTCATCATCAGCATCTCTTATCAAATTCATATACCTAAGAAAGCTAGTCCCGAAACGGTAGtt        >  2:95909/1‑149 (MQ=255)
                                                                                                                                        aaTTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTAAGCTTAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTaa  <  2:470739/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
TCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTAGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAA  >  NC_000913_3_hsa_tpiA/2106597‑2106881

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: