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% 118.5 / ‑6.7 42P725L (CCG→CTG) acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (20/21);  total (21/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.

TGACATAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCG  >  NC_000913_3_hsa_tpiA/482084‑482370
                                                                                                                                                  |                                                                                                                                            
tgacatAAACTGTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGTCTGCGCTTTTTCCTGGTCTATATCAATCTTAAACTGCagg                                                                                                                                            <  2:605608/149‑1 (MQ=255)
          tttCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCGGGGCGCAGGTTGGTGTTAATGTCGAGGATAGAAACACCCCGCGCCTGCGCTTTTTCCTGGGCGATATCAATCGTAAACTGCAGGGTATCttgct                                                                                                                                  >  2:519200/1‑146 (MQ=255)
          tttCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCa                                                                                                                                  >  1:133264/1‑149 (MQ=255)
                  cacGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGttt                                                                                                                          <  1:595106/149‑1 (MQ=255)
                  cacGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGttt                                                                                                                          <  1:686307/149‑1 (MQ=255)
                   aggACGGCGTTCGATAAAGTCGTCCACATAGCTGCCTCCCCATGCAGCGCCCAGAGGGGTGTTAAGGTCGTTGATAGAAGCACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAGCTGCAGGGTATCTTTCAGACCGTTTg                                                                                                                         <  1:358606/147‑1 (MQ=255)
                    cGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGATATCTTCCAGACCGTTTgg                                                                                                                        >  1:9910/1‑149 (MQ=255)
                     gACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGa                                                                                                                       >  2:377497/1‑149 (MQ=255)
                         gcgGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTAGTCGCTATCAATCTTAAACTGCAGGGTATCTTCGAGACCGTTTGGACGTa                                                                                                                   >  2:349310/1‑149 (MQ=255)
                            gTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGc                                                                                                                <  1:58349/149‑1 (MQ=255)
                             tCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTTGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCt                                                                                                               >  1:340245/1‑149 (MQ=255)
                                   aaGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCaa                                                                                                         >  1:554702/1‑149 (MQ=255)
                                      tCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACat                                                                                                      <  1:3587/149‑1 (MQ=255)
                                        gttagcttAGCTGCCGCCGCATGTAGCGCCCAGAGTGTTGTTAATGTCGTTGATAGAACCACCCAGCGCCTGCGCTTTTTCCTTGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACatat                                                                                                    <  2:30798/142‑1 (MQ=255)
                                            acaTAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCGAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCCGTTGGACGTACGCTGGTCAACATATCAgg                                                                                                >  2:153094/1‑149 (MQ=255)
                                            acaTAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAGACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAgg                                                                                                <  2:260917/149‑1 (MQ=255)
                                                        ccccATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGCTTCTTcgccgc                                                                                    >  2:358793/1‑149 (MQ=255)
                                                            aTGCAGCGTCCGGAGTGGTGTTAATGGCGTTGATAGAAACACCCAGCCCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCt                                                                                <  1:79527/149‑1 (MQ=255)
                                                                aGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTgcaa                                                                            <  1:109452/149‑1 (MQ=255)
                                                                 gcgcCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTgcaag                                                                           <  1:403538/149‑1 (MQ=255)
                                                                    cccAGAGTGGTGTTAATGTCGTTGATAGAACCACCCAGCGCCCGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTgcaagcaa                                                                        <  1:470592/149‑1 (MQ=255)
                                                                    cccAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTgcaagcaa                                                                        <  1:477883/149‑1 (MQ=255)
                                                                      cAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGc                                                                                            >  1:286068/1‑127 (MQ=255)
                                                                      cAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGc                                                                                            <  2:286068/127‑1 (MQ=255)
                                                                        gagTGGTGTTAATGTCGTTGATAGAAACCCCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTgg                                                                    <  1:259314/149‑1 (MQ=255)
                                                                              tgtTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATATTCCAGACCGTTTGGACGTAGGCTGGTCAACATATCAGGGTGCTTAGCGGCTTCTGCAAGTAACTGGTTAcgc                                                              >  2:557340/1‑149 (MQ=255)
                                                                                   aTGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATAGCAATCTTAACCTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGGGCCTCGCTGCGTCTGCAATCAACTGGTAACGCGACTg                                                         >  2:358851/1‑149 (MQ=255)
                                                                                     gTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGACTGAg                                                       >  2:104234/1‑149 (MQ=255)
                                                                                            aTAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCCATCTTAAACTGCAGGGTATCTTCCAGACCGTTTTGACGTACGCTGGTCAACATATCAGGGTGCTTCACTGCTTCTGCAAGCAACAGGTTACGCGCCTGAGTCAGttt                                                >  1:305023/1‑149 (MQ=255)
                                                                                                    acCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGAc                                        <  1:143841/149‑1 (MQ=255)
                                                                                                       cAGCGCCTGCGCTTTTTCCTAACCCAGAGCAATCTAAAACTCCAGGGTATCTTCCAGACCGTTTGGACGTACCCTGGCCAACATATCAGCGTGCCCCGCTCCTTCTGCAAGCAACTGGTTACGCCCCTGAGTCAGTTTTTCGTGCCCaa                                     >  2:523180/1‑149 (MQ=255)
                                                                                                        aGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAg                                    >  1:434750/1‑149 (MQ=255)
                                                                                                                       tCCTGGGCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCTACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                     <  1:610936/149‑1 (MQ=255)
                                                                                                                           ggTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGGCAACATATCAGGGTGCTCCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGc                 <  1:107356/149‑1 (MQ=255)
                                                                                                                            gTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCGATCAGCt                >  2:645882/1‑149 (MQ=255)
                                                                                                                             tCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGGCGTACGCTGGGCAACATATCAGGGGGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGACAGTTTTTCGTGAGCAAGGCCAGCCTGGTCAATCAGCtc               >  2:461899/1‑149 (MQ=255)
                                                                                                                                atatCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaa            <  1:214170/149‑1 (MQ=255)
                                                                                                                                atatCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaa            <  2:133264/149‑1 (MQ=255)
                                                                                                                                atatCAATCTCAAACTGCAGGGTATCTTCCACACCGTTGGGACGTACGCTGGTCAACATCTCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaa            <  1:186339/149‑1 (MQ=255)
                                                                                                                                       tCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtaaaagtcaaag     >  1:520885/1‑149 (MQ=255)
                                                                                                                                         ttAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCGTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCACCCTGGTCAATCAGCtcaaagtaaaagac   >  1:399281/1‑147 (MQ=255)
                                                                                                                                         ttAAACTGCAGAGTAACTACCAGACCGTTAAAACGGACGCTGGTCAACATATCAGTGTGCTTCGTTTCTTCTGCAAGCAACTGGTTACACGCCTGAGTCAGTTTTTGCTGACCAAGACCAAGCTCGACGATCTGCTCaaagccaaagcc   >  2:62504/1‑149 (MQ=255)
                                                                                                                                          tAAACTGCAGGGTATCTTCCAGCCCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCg  <  1:448101/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
TGACATAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCG  >  NC_000913_3_hsa_tpiA/482084‑482370

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