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% 136.5 / NA 48P725L (CCG→CTG) acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/25);  total (23/25)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ATAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGC  >  NC_000913_3_hsa_tpiA/482088‑482375
                                                                                                                                              |                                                                                                                                                 
atAATCTTTCTTCACACGGCCGCGGTCGATAAAGTCGTTCAAATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCAAAATTTCACAACTTGATCTTGTCGATATCAATCTTAAACTGCAGGGTAt                                                                                                                                             >  2:487481/1‑149 (MQ=255)
atAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTAt                                                                                                                                             >  1:594432/1‑149 (MQ=255)
      tttCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGTTCGTTTTCTATCTTAACCTCCAGGGTATATTcag                                                                                                                                       >  2:534985/1‑147 (MQ=255)
      tttCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCa                                                                                                                                       <  1:87909/149‑1 (MQ=255)
           tcACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGAGAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAGACTGCAGGGTATCTTCCAGACCg                                                                                                                                  <  2:49431/149‑1 (MQ=255)
             acacGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGCTAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGGATCTTCCAAACCGtt                                                                                                                                >  2:330334/1‑149 (MQ=255)
                 gACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGa                                                                                                                            <  2:52264/149‑1 (MQ=255)
                        gTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGAGGATAGTGGCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGc                                                                                                                     <  1:537751/149‑1 (MQ=255)
                           gATAACGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGGTAATGTCGTTGATAGAAACACCCAGCGCCTGCACTTTTTCCTGGTCGACATCAATCTTAAACTGCATGCTATCTTTCACACCGTTTGGACGTaacatgg                                                                                                                  >  2:608508/1‑143 (MQ=255)
                               gaGTCGTCCACCGAGCTGCCGCCCCACGCCGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCTCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCaa                                                                                                              <  2:325501/148‑1 (MQ=255)
                               aaGTTGTTCACATAGCTGTCGCCCCATGCAGCGCCCAGAGGGGTGTTAATGTCGTTGATAGAAACACCCAGCGCGTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCaa                                                                                                              <  2:75318/149‑1 (MQ=255)
                               aaGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGCCGTCCGCTGTTCaa                                                                                                              >  1:552116/1‑149 (MQ=255)
                                       cacaTAGCTGCCGCCCCATGCAGCGCCCAGAGTACGTTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTTTTTTTTTTTATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGCCGTTCGCTGGTCAGCATCTCAg                                                                                                      >  2:493005/1‑149 (MQ=255)
                                        acaTAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTCCGCTGGTCAACATCTCAgg                                                                                                     >  1:193132/1‑149 (MQ=255)
                                        acaTAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAgg                                                                                                     >  1:587721/1‑149 (MQ=255)
                                        acaTAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACGGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAgg                                                                                                     <  2:557099/149‑1 (MQ=255)
                                            aGCTGCCGCCCCATGCAGCGCCAAGAGTGGTGTTAATGTCGTTGATAGCAACACCCAGTGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGc                                                                                                 <  2:73907/149‑1 (MQ=255)
                                               tGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGGCGTTGATAGAAACACCCAGCGACTGCGCATTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTGCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTc                                                                                              >  1:168126/1‑149 (MQ=255)
                                               tGCCGCCCCATGCAGCGCCCAGAGGGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTc                                                                                              >  1:596466/1‑149 (MQ=255)
                                                 ccgccCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCgc                                                                                            >  2:6893/1‑149 (MQ=255)
                                                 ccgccCCATGCAGCGCCCAGAGTGGTGTTAATGACGTTGATATAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCgc                                                                                            <  1:126540/149‑1 (MQ=255)
                                                  cgccCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGCCCGTTTGGACGTACGCCAGTCAACATATCAGGCTGCTTCcct                                                                                           >  2:450417/1‑149 (MQ=255)
                                                  cgccCCATGCAGCGCCCAGAGTGGGTTTAATGTCGGTGATAGAAACACCCAGCGCCTGCGCTTTTTCCGGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCgct                                                                                           <  2:516450/149‑1 (MQ=255)
                                                      ccATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCgctgcct                                                                                       >  1:152299/1‑147 (MQ=255)
                                                      ccATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCtt                                                                                       >  1:233732/1‑149 (MQ=255)
                                                                 ccAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGGCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAAc                                                                            <  2:643244/149‑1 (MQ=255)
                                                                  agggggtgtggTAATGTGGGTGATGGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATAGCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTACTTCGCTGCTTCTGCAAGCAACt                                                                           <  2:375139/142‑1 (MQ=255)
                                                                          tgtTTATGTTGTTGATAGAAACACCCAGTGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGAGGTACGCTGGTCAACATATTAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAcgc                                                                   <  1:371744/149‑1 (MQ=255)
                                                                          tgtTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAccc                                                                   >  1:254170/1‑147 (MQ=255)
                                                                            ttAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACgcgc                                                                 <  1:93763/149‑1 (MQ=255)
                                                                                    gTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGAGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCa                                                         >  2:603369/1‑149 (MQ=255)
                                                                                       gATAGAAACACCCAGTGCTTGGCCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCACtt                                                      <  2:120376/149‑1 (MQ=255)
                                                                                            aaaCACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTCCTTCTGCAAGCAACTGGTTACGCGCCTGAGCCAGTTTTTct                                                 >  2:461211/1‑148 (MQ=255)
                                                                                              acacTCAGCGCCTGCGCTTTTTCTGGGGCGATATTAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTg                                               <  1:370807/149‑1 (MQ=255)
                                                                                               cacCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTCCTGCGATGAGCTGGTTACGCGCCTGAGTCAGTTTTTCGTGa                                              <  1:510205/149‑1 (MQ=255)
                                                                                                   cAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAACCAACTGTTTACACTCCTGAGTCAGTTTTTCTTGACCaa                                          >  1:5108/1‑149 (MQ=255)
                                                                                                                   gccTGGTCGATATCTATCTTAGACTGCAGGGTATCTTCCAGTGCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCGTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                          <  2:460886/148‑1 (MQ=255)
                                                                                                                     cTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCa                        <  1:378374/149‑1 (MQ=255)
                                                                                                                     cTGGTCGATAGTAATCTGAAACTGCAGGGTATCCTCCAGATCGGTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCa                        <  2:115007/149‑1 (MQ=255)
                                                                                                                      gggTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAg                       >  2:465238/2‑149 (MQ=255)
                                                                                                                       ggTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGc                      <  1:599480/149‑1 (MQ=255)
                                                                                                                            atatCAATCTTAAACTGCAGTGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaa                 <  2:445178/149‑1 (MQ=255)
                                                                                                                              atCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTTGTCAACATCTCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCGGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagt               <  2:91693/149‑1 (MQ=255)
                                                                                                                              atCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagt               >  1:369931/1‑149 (MQ=255)
                                                                                                                                caagcgtAAACTGGAGGGGATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtca             <  1:318180/143‑1 (MQ=255)
                                                                                                                                       aaaCTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCgg      >  2:252528/1‑149 (MQ=255)
                                                                                                                                          cTGCAGGGTATCTGCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCCAAGCCGGTTg   <  1:672158/149‑1 (MQ=255)
                                                                                                                                           tGCAGGGTATCTTCCAGACCGATTGGACGTACTCTGGTCAACATATCAGGGTGCTTCGCTCCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGGCTGGTCAATCAGCGTAAAATCAACGCCGCTTGc  >  2:250973/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
ATAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGC  >  NC_000913_3_hsa_tpiA/482088‑482375

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