Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_hsa_tpiA 482,230 G→A 100% P725L (CCG→CTG)  acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA482,2300GA100.0% 117.0 / ‑6.7 43P725L (CCG→CTG) acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (22/19);  total (23/20)
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.

TGACATAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCC  >  NC_000913_3_hsa_tpiA/482084‑482369
                                                                                                                                                  |                                                                                                                                           
tgacatAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCagg                                                                                                                                           <  1:227880/149‑1 (MQ=255)
    atAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGGAACACCCAGCGCCGGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTAt                                                                                                                                       >  1:49400/1‑149 (MQ=255)
         cTTTCTTCACACGACCGCGGTCGATAAAGTGGTTCACATAGCTGACGCCCCATGCAGCGCCCAGAGTGGTGTTAATATCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCGGGTCGATATCAATATTAAACTGCAGGGTGTCTTcc                                                                                                                                  >  2:483040/1‑149 (MQ=255)
         cTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAACCACCCAGTTCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTACCTTcc                                                                                                                                  >  1:269451/1‑149 (MQ=255)
                  cacgcacgcgGTCGATAAACTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGttt                                                                                                                         <  1:13832/143‑1 (MQ=255)
                     caCCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCACAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGa                                                                                                                      <  1:157006/148‑1 (MQ=255)
                      tcccctgTCGATAAAGTCTGTCACATAGCTGCCGCCCCATGGAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTAGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGAc                                                                                                                     <  2:385911/143‑1 (MQ=255)
                        cgcgGGCGATAAAGGCGTGCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGC‑TTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTa                                                                                                                  <  2:319269/149‑1 (MQ=255)
                         gcgGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTa                                                                                                                  <  1:460682/149‑1 (MQ=255)
                             tCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCAAGAACGTTTGGACGTACGCt                                                                                                              >  1:460907/1‑149 (MQ=255)
                             tCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAACGCGATGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTActct                                                                                                              >  1:299497/1‑149 (MQ=255)
                                  aaaGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCa                                                                                                         >  1:516746/1‑149 (MQ=255)
                                   aaGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCaa                                                                                                        <  2:320393/149‑1 (MQ=255)
                                    aGTCGTTCACATAGTTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAAc                                                                                                       <  2:606584/149‑1 (MQ=255)
                                              aTAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGt                                                                                             >  1:233745/1‑149 (MQ=255)
                                                      ccccccATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCgct                                                                                     <  1:332667/147‑1 (MQ=255)
                                                         cccATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCgctgct                                                                                  >  1:480184/1‑149 (MQ=255)
                                                          ccATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCtt                                                                                 >  1:134495/1‑149 (MQ=255)
                                                           cATGCAGCGCCCAGAGTGGTGTTGATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTc                                                                                <  1:433348/149‑1 (MQ=255)
                                                              gCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTgc                                                                             <  1:553694/149‑1 (MQ=255)
                                                                aGCGCCCAGAGGGGTGTTAATGTCGTGGATAGAAACCCCCCGCGCCTGCGCTTTTTCCTTGTCGATATCAATCTTAAACTGCAAGGTCTTTTCCATACCGTTTGGCCGTGCGCTTGTCAACATATCAGGGTGCTTCGCTGCTTCTgcaa                                                                           <  1:241357/149‑1 (MQ=255)
                                                                    cccAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTgcaagcaa                                                                       >  1:189050/1‑149 (MQ=255)
                                                                           tggtgTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTCGCCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTa                                                                <  1:423129/149‑1 (MQ=255)
                                                                                ttAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACTTACACTGGGCAACATATCAGCGTGCTTCGCTGCTTCTGCAACCAACTGGTTACgcgc                                                           >  2:219613/1‑149 (MQ=255)
                                                                                  aaTGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCt                                                         <  1:68685/149‑1 (MQ=255)
                                                                                        gTTAATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACACGCCTGAGTCa                                                   >  1:500660/1‑149 (MQ=255)
                                                                                            agagagACACCGAGCGCCTGCACTTGTTCCAGGTCGATATCAATCTTAAACTGCAGGGTATCTTGCAGACCGTTTGGAAGTACGCTGATCAAACTATCAGGGTGATTCGCTGCTTCTGCAAGGAACAGGTTACGAGCGCGAGACAGttt                                               >  2:8509/3‑149 (MQ=255)
                                                                                            aTAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGttt                                               >  1:571705/1‑149 (MQ=255)
                                                                                                   cacCCCGGGCCTGCCCTTTGTCCTGGGTCATAGCAATCGTAAACTGCAGTGTATGTTCCCGACCGGTTGGACGTATGCTGGTTAACATATCAGGGTGCTTCGCTTTTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGa                                        <  2:156009/149‑1 (MQ=255)
                                                                                                   cacCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGa                                        <  2:439577/149‑1 (MQ=255)
                                                                                                        aGCGCCTGCTCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGTGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAg                                   <  1:38531/149‑1 (MQ=255)
                                                                                                        aGCGCCCGCGGTTTGTCCTGGTCGATATCAATCTTAAGCTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAg                                   <  2:80323/149‑1 (MQ=255)
                                                                                                           gcCTGCGCTTTTTCATGGTCGATATCAATCTTAAAAAGCAGGGCCTATTCCAGACAGTGAGGACGTACGCTGGTCAACATATCATGCTGCTTCACTGCTTCTGCAAGAACCCGGTTACGCGCCTGAGTCACTTTTTCCTTACCAACGcc                                >  2:225422/1‑149 (MQ=255)
                                                                                                             cTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGAAGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCACGCAACTGGTTACGCCCCTGAGTCAGTTTTTCGTGACCAAGGCCAg                              >  2:537274/1‑149 (MQ=255)
                                                                                                                   tttttCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGTCTGGt                        >  1:213835/1‑149 (MQ=255)
                                                                                                                          tGGTCGATATCAATCTTAAACTGCAGGGTAGCTTCCAGACCGTTTGGACGTACGATGGTCAACATATCAGGGTGCTTTGCTGCTTCTTGAAAGAACTGGTTACGCGCCTGAGCCAGCTTTTCGTGACCACGGCCAGCCTGGTAAATCaa                 >  2:99386/1‑148 (MQ=255)
                                                                                                                            gTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTCGTTACGCGCCTGAGTCAGTTATTCGTGGGCGAGGCGAGCCGCGTCaagtaggg               >  2:427927/1‑143 (MQ=255)
                                                                                                                                ataaCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGTCTTAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaa           <  2:565258/149‑1 (MQ=255)
                                                                                                                                   tCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtc        >  2:83249/1‑149 (MQ=255)
                                                                                                                                      aTCTTAAACTGCTGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGCGGCTGTACCCGGGCTAGTCAATCAGCtcaaagtcaaa     <  2:132177/149‑1 (MQ=255)
                                                                                                                                       tCTTAAACTGCAGGGTATCTTCCAGCCCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtcaaag    <  1:195875/149‑1 (MQ=255)
                                                                                                                                       tCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtcaaag    >  1:239795/1‑149 (MQ=255)
                                                                                                                                       tCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCCGGTCAATCAGAccaaagtcaaag    >  2:267490/1‑149 (MQ=255)
                                                                                                                                         ttAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTCGTTTTTATACAAGGTTTCTTTGTTGCGCCTGCACGCAACTGGTTAAGCGCCTGAGTAAGTTTTTCGTAACCCAGACCATCCTAGTCAAACAGCtcaaaaacaaagac  >  2:283233/1‑147 (MQ=255)
                                                                                                                                                  |                                                                                                                                           
TGACATAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCC  >  NC_000913_3_hsa_tpiA/482084‑482369

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