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

TGACATAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTG  >  NC_000913_3_hsa_tpiA/482084‑482374
                                                                                                                                                  |                                                                                                                                                
tgacatAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATCGCTGCCGCCCCATGTAGTGCCCAGAGTGGTGTTAATGCCGTTGATAGAAAGACCAAGCGCCTGCGCTTTTTCCTGGTCGATATTAGTCGTAAACTGCagg                                                                                                                                                <  2:649881/149‑1 (MQ=255)
tgacatAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCagg                                                                                                                                                <  2:136201/149‑1 (MQ=255)
  acatAAACTTTCTTCACACGACCGCGGTCGATACAACAGATCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCCATCTTAAACTGCAGGGt                                                                                                                                              >  2:516853/1‑149 (MQ=255)
    atAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCAATGCAACGACCAAAGTGTTGTTAATCTCGCTGATAACAACACCAAGCGCCGGCTCGCTCTCCTGGTCGATATCAATCTTAAACTGCaggcgat                                                                                                                                            >  2:334513/1‑145 (MQ=255)
    atAAACTTTCTCCACACGACCGCGGACGATCACGTCGACCACATATCTGCCGCCCCATGCAGCGCCCAGAGTGGTCCTACTGCCGGTGATTGAAACACCCAGCGCCTGCACTTTTTCCTAGTCGATATCAATCTTAAACTGCCGGGCAt                                                                                                                                            >  2:126236/1‑149 (MQ=255)
        acgttcttcACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTc                                                                                                                                        <  2:500162/146‑1 (MQ=255)
         cTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTcc                                                                                                                                       <  2:455987/149‑1 (MQ=255)
         cTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTcc                                                                                                                                       >  2:102919/1‑149 (MQ=255)
          tttCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCa                                                                                                                                      <  1:101130/149‑1 (MQ=255)
              ttcACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGAcc                                                                                                                                  <  1:543225/149‑1 (MQ=255)
                 acacGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGtt                                                                                                                               >  2:458361/1‑149 (MQ=255)
                  cacGACCGCGGTCGATAAAACCGTTCACATACCCGCAGACCCCCCCAGCGCCCAGAGTGGCGTCAATGTCGTCGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGCCGATATCAATCTTAAACCCCAGGGCCTCTTACCGACCGttt                                                                                                                              >  2:305205/1‑149 (MQ=255)
                           ggTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATGTTCCAGACCGTTTGGACGTACg                                                                                                                     >  2:44123/1‑149 (MQ=255)
                           ggTCGATAAAGTCGTTCACATAGCTGCCGCCCCACGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTAGATATCAATCTTAAACTGCAGGGGATCTTCCAGACCGTTTGGACGTACg                                                                                                                     >  1:645297/1‑149 (MQ=255)
                             tCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTCCAGGGTATCTTCCAGCCCGTCTGGCCGTACGCt                                                                                                                   >  2:203627/1‑149 (MQ=255)
                             tCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTCGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCt                                                                                                                   >  2:687761/1‑149 (MQ=255)
                                    aGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAAc                                                                                                            >  1:310730/1‑149 (MQ=255)
                                     gTCGTTCACATAGCTGCCGCCCCATGCAGCGCCGAGCGTGGTTTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACa                                                                                                           <  1:494732/149‑1 (MQ=255)
                                         ttCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATc                                                                                                       >  1:456776/1‑149 (MQ=255)
                                              aTAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAAGACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAGCATATCAggga                                                                                                  >  1:474612/1‑148 (MQ=255)
                                              aTAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGt                                                                                                  >  2:276365/1‑149 (MQ=255)
                                                       gccCCATGCAGCGCCCAGAGTTGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCgctg                                                                                         >  2:423540/1‑149 (MQ=255)
                                                       gccCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTgg                                                                                                                            >  1:407402/1‑114 (MQ=255)
                                                       gccCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTgg                                                                                                                            <  2:407402/114‑1 (MQ=255)
                                                        ccccATGCAGCGCCCAGAGTGGTGTTAATTTCGTTGATAGAACCACAGAGCGACTGCGCTTTGTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCgctgc                                                                                        <  1:497770/149‑1 (MQ=255)
                                                                     ccAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAAc                                                                           >  1:231772/1‑149 (MQ=255)
                                                                     ccAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAAc                                                                           >  1:556804/1‑149 (MQ=255)
                                                                       agagTGGTGTTAATGTCGTGGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTCGTCAACATATCAGGGGTCTTCTCTACTTCTGCACGCAACTg                                                                         >  2:609476/1‑149 (MQ=255)
                                                                              tgtTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGa                                                                                                                           >  1:345045/1‑92 (MQ=255)
                                                                              tgtTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGa                                                                                                                           <  2:345045/92‑1 (MQ=255)
                                                                                            aTAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAcgc                                                                  >  2:562243/1‑135 (MQ=255)
                                                                                            aTAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAcgc                                                                  <  1:562243/135‑1 (MQ=255)
                                                                                            aTAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGttt                                                    >  1:490876/1‑149 (MQ=255)
                                                                                             tAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGtttt                                                   <  1:687761/149‑1 (MQ=255)
                                                                                                 aaCACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACatat                                                                                                        <  1:653502/92‑1 (MQ=255)
                                                                                                 aaCACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACatat                                                                                                        >  2:653502/1‑92 (MQ=255)
                                                                                                  acacCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAAc                                                                           <  2:294621/120‑1 (MQ=255)
                                                                                                  acacCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAAc                                                                           >  1:294621/1‑120 (MQ=255)
                                                                                                  acacCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTg                                              <  2:490876/149‑1 (MQ=255)
                                                                                                   cacCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGa                                             <  2:500751/149‑1 (MQ=255)
                                                                                                      ccAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCa                                          >  1:611123/1‑149 (MQ=255)
                                                                                                       aggcgcCTGCGCTTTTCCCTGGTCGAGATCAATCTTAAACTGCAGGGTATCTTGCAGACCTTTTGGACGTAGGCTGGGCAACATATAAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCATGGGTCAGTTTTTCGTGCCCac                                         <  2:551640/147‑2 (MQ=255)
                                                                                                            ccTGCGCTTTTTTTTGGTTGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCa                                    >  2:394293/1‑149 (MQ=255)
                                                                                                                 gcTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGCATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTg                               >  2:188889/1‑149 (MQ=255)
                                                                                                                   tttttCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGt                             <  1:306914/149‑1 (MQ=255)
                                                                                                                   tttttCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGt                             >  2:73553/1‑149 (MQ=255)
                                                                                                                          tGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAg                      >  1:450097/1‑149 (MQ=255)
                                                                                                                          tGGTCGATATCAATCGTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAg                      <  2:295679/149‑1 (MQ=255)
                                                                                                                                     caTCTCAAACCACAGGGTACCTTCGATACCGTTTGGACGTACGCTGGCCAACGTATCCGAGTGCTTTGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTTACTTTTTCGTGACGAAGCCCAGCATGGTCAATCAGCtcatagtcaa           <  2:660172/148‑1 (MQ=255)
                                                                                                                                       tCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtcaaag         >  1:676313/1‑149 (MQ=255)
                                                                                                                                           aaaCTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCgg     <  2:102491/149‑1 (MQ=255)
                                                                                                                                             aCTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGtt   <  2:310730/149‑1 (MQ=255)
                                                                                                                                             aCTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGtt   <  1:677375/149‑1 (MQ=255)
                                                                                                                                              cTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTg  <  1:276365/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TGACATAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTG  >  NC_000913_3_hsa_tpiA/482084‑482374

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