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,2300GA98.0% 170.8 / NA 50P725L (CCG→CTG) acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (21/28);  minor base C (1/0);  total (22/28)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGC  >  NC_000913_3_hsa_tpiA/482088‑482375
                                                                                                                                              |                                                                                                                                                 
atAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCAGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTAt                                                                                                                                             >  2:153459/1‑149 (MQ=255)
     ctctcttcACACGACCGCGGTCGATACAGCCGTTCCCATAGCTCCCGCCCCCTGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCCGCGCCTGCGTTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTcc                                                                                                                                        <  1:295736/146‑1 (MQ=255)
     cTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCCGGGTATCTTcc                                                                                                                                        >  2:364653/1‑149 (MQ=255)
      cttcttcTCACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCa                                                                                                                                       <  2:370520/148‑1 (MQ=255)
              cccGACCTCGGTCGATAAAGGCGTTCACATTGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGttt                                                                                                                               <  2:121993/147‑1 (MQ=255)
              cacGACCGCGGTCGATACAGTCATACAGATAGCTGCCGCCCCATGAAGCGCCCAGACTGGTGTTAATGTCGTCGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCATGGTATCTTCCTGACCGttt                                                                                                                               <  2:698855/149‑1 (MQ=255)
              cacGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGttt                                                                                                                               >  1:684930/1‑149 (MQ=255)
               ccGACCTCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGCAAAACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTg                                                                                                                              <  1:368469/148‑1 (MQ=255)
               acGCCCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGGGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTg                                                                                                                              <  2:564962/149‑1 (MQ=255)
                  aCCGCGGTCGATAAAGTCGTGCACATAGCTCCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGAc                                                                                                                           <  2:357638/149‑1 (MQ=255)
                     gcgGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTa                                                                                                                        <  2:156629/149‑1 (MQ=255)
                        gTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGc                                                                                                                     <  1:597779/149‑1 (MQ=255)
                             tAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTc                                                                                                                <  1:512439/149‑1 (MQ=255)
                                  tCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACat                                                                                                           <  2:817346/149‑1 (MQ=255)
                                        acaTAGCTGCCGCCCCATGCCGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCGGGTCGATATCAATCTTAAACTGCAGGGTAACTTCCAAACCGTTTGGACGTACGCTGTGCAACATATCAgg                                                                                                     >  2:568031/1‑149 (MQ=255)
                                        acaTAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAgg                                                                                                     <  1:640798/149‑1 (MQ=255)
                                              cTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCtt                                                                                               >  1:517048/1‑149 (MQ=255)
                                               tcccgccCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTc                                                                                              <  1:135238/147‑1 (MQ=255)
                                                 ccgccCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCgc                                                                                            >  2:745959/1‑149 (MQ=255)
                                                       cATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCagg                                                                                                                                                 <  2:742920/90‑1 (MQ=255)
                                                            aGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTgcaa                                                                                 <  1:745959/149‑1 (MQ=255)
                                                                cccAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGtt                                                                                                                                >  2:624525/1‑98 (MQ=255)
                                                                cccAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGtt                                                                                                                                <  1:624525/98‑1 (MQ=255)
                                                                    gagTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTgg                                                                         <  1:694137/149‑1 (MQ=255)
                                                                            ttAATGTCGTTGATAGGAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACgcgc                                                                 >  1:105450/1‑149 (MQ=255)
                                                                              aaTGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGa                                                                                                                            <  2:137944/88‑1 (MQ=255)
                                                                              aaTGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGa                                                                                                                            >  1:137944/1‑88 (MQ=255)
                                                                                        aTAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGttt                                                     >  2:257463/1‑149 (MQ=255)
                                                                                          aGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTACTGCAAGCGACTAGTTACGCGCCTGAGTCAGgtttt                                                   >  1:228988/1‑149 (MQ=255)
                                                                                            aaaCACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCg                                                 >  1:39147/1‑149 (MQ=255)
                                                                                                acCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGAc                                             >  1:796051/1‑149 (MQ=255)
                                                                                                          tGCGCTTTTTCCTGGTCGATATCAATCTTAAACGGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGc                                   <  1:339177/149‑1 (MQ=255)
                                                                                                            cgcTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCt                                 <  2:517048/149‑1 (MQ=255)
                                                                                                            cgcTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCt                                 <  2:684930/149‑1 (MQ=255)
                                                                                                             gcTGTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTg                                <  2:493353/149‑1 (MQ=255)
                                                                                                               tttttCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGt                              >  1:110025/1‑149 (MQ=255)
                                                                                                               tttttCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGt                              <  1:482324/149‑1 (MQ=255)
                                                                                                                      tGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGCAGAATGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATcca                       >  2:137899/1‑147 (MQ=255)
                                                                                                                        gTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCt                     >  1:762264/1‑149 (MQ=255)
                                                                                                                                 aaTCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGc                      <  2:254737/139‑1 (MQ=255)
                                                                                                                                 aaTCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGc                      >  1:254737/1‑139 (MQ=255)
                                                                                                                                 aaTCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagccaa            >  1:761713/1‑149 (MQ=255)
                                                                                                                                  attttAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtcaaa           <  1:570399/146‑1 (MQ=255)
                                                                                                                                  aTCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtcaaa           <  1:570468/149‑1 (MQ=255)
                                                                                                                                  aTCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCACTTTTTCGTGACCTAGGCCAGCCTGGTCAATCAGCtcaaagtcaaa           <  1:656914/149‑1 (MQ=255)
                                                                                                                                   tCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATGTCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtcaaag          >  1:780525/1‑149 (MQ=255)
                                                                                                                                   tCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGt                                                <  2:255843/111‑1 (MQ=255)
                                                                                                                                   tCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGt                                                <  2:258452/111‑1 (MQ=255)
                                                                                                                                   tCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGt                                                >  1:258452/1‑111 (MQ=255)
                                                                                                                                   tCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGt                                                >  1:255843/1‑111 (MQ=255)
                                                                                                                                           tGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGc  >  1:95408/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 < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑