Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_hsa_tpiA 482,230 G→A 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,2300GA93.0% 130.8 / 1.1 45P725L (CCG→CTG) acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (23/17);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 5.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

TAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCA  >  NC_000913_3_hsa_tpiA/482089‑482376
                                                                                                                                             |                                                                                                                                                  
tAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCt                                                                                                                                            >  1:292968/1‑150 (MQ=255)
  aaCTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTc                                                                                                                                          >  1:159274/1‑150 (MQ=255)
     tttCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAg                                                                                                                                       >  2:212682/1‑150 (MQ=255)
     tttCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTCTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAg                                                                                                                                       >  1:656/1‑150 (MQ=255)
      ttcttcACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGa                                                                                                                                      <  2:404777/150‑1 (MQ=255)
      ttcttcACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGa                                                                                                                                      <  1:314062/150‑1 (MQ=255)
        cttcACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGAcc                                                                                                                                    <  3:29519/150‑1 (MQ=255)
              acGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTgg                                                                                                                              >  1:46674/1‑150 (MQ=255)
               cGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGCCACCCCCAGCGCCTGCTGTTTTTCCTGATCGAGATCAATCTTAATCTTCTGGTTATCTTCCAGACCGTTTGGa                                                                                                                             <  2:59398/150‑1 (MQ=255)
                        tCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTg                                                                                                                    >  2:197702/1‑150 (MQ=255)
                          gATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGt                                                                                                                  >  2:24723/1‑150 (MQ=255)
                                                ccgccCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCgct                                                                                            <  3:89131/150‑1 (MQ=255)
                                                           aGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTgcaag                                                                                 >  1:425146/1‑150 (MQ=255)
                                                              gcCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGTTTCGCTGCTTCTgcaagcaa                                                                              >  1:208590/1‑150 (MQ=255)
                                                              gcCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTgcaagcaa                                                                              >  4:226909/1‑150 (MQ=255)
                                                                 cAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGTGTGCTTCGCTGCTTCTGCAAGCAACTg                                                                           <  1:304423/150‑1 (MQ=255)
                                                                    agTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGtt                                                                        <  1:120226/150‑1 (MQ=255)
                                                                         tgtTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAcgcg                                                                   >  4:204080/1‑150 (MQ=255)
                                                                          gtTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACgcgc                                                                  <  1:212682/150‑1 (MQ=255)
                                                                              aTGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGa                                                              <  3:254814/150‑1 (MQ=255)
                                                                                  cGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCa                                                          <  3:193982/150‑1 (MQ=255)
                                                                                    ttGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGt                                                        >  2:275645/1‑150 (MQ=255)
                                                                                      gATATAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGttt                                                      <  3:56523/150‑1 (MQ=255)
                                                                                        tAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCa                                                          >  1:406389/1‑144 (MQ=255)
                                                                                        tAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCa                                                          <  2:406389/144‑1 (MQ=255)
                                                                                        tAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCCGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGCGTCAGttttt                                                    >  2:408452/1‑150 (MQ=255)
                                                                                         aGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGATGTACGCTGGGCTACAGAGACCGGTGCTTTGTTGACTGTGCAAGAAACTCGTTACGTACATGAGTCAGtttcac                                                   <  2:208348/150‑4 (MQ=255)
                                                                                         aGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTATGCGCCTAAGGCAGTTTTTc                                                   <  2:46674/150‑1 (MQ=255)
                                                                                             acacCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGa                                               <  4:21209/150‑1 (MQ=255)
                                                                                              cacCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGAc                                              >  3:60524/1‑150 (MQ=255)
                                                                                                   aGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAgg                                         >  1:99064/1‑150 (MQ=255)
                                                                                                        cTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGc                                    <  3:244343/150‑1 (MQ=255)
                                                                                                         tGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGt                                                 >  2:304790/1‑136 (MQ=255)
                                                                                                         tGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACGGGTTACGCGCCTGAGTCAGTTTTTCGt                                                 <  1:304790/136‑1 (MQ=255)
                                                                                                                tttCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCaa                            >  1:375889/1‑150 (MQ=255)
                                                                                                                    cTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAg                        <  2:159274/150‑1 (MQ=255)
                                                                                                                      ggTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGGCAATCAGCt                      >  3:91801/1‑150 (MQ=255)
                                                                                                                        tCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATc                          <  2:384615/144‑1 (MQ=255)
                                                                                                                        tCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATc                          >  1:384615/1‑144 (MQ=255)
                                                                                                                              tCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtca              >  2:112779/1‑150 (MQ=255)
                                                                                                                                    ttAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCg        >  3:151234/1‑150 (MQ=255)
                                                                                                                                     tAAACTGCAGGGTAGCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCGGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGGGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCgg       >  4:133530/1‑150 (MQ=255)
                                                                                                                                       aaCTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGtt     <  1:236354/150‑1 (MQ=255)
                                                                                                                                          tGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCa  >  3:81497/1‑150 (MQ=255)
                                                                                                                                          tGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCa  >  3:134204/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
TAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCA  >  NC_000913_3_hsa_tpiA/482089‑482376

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