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_tpiA	482,230	0	G	A	100.0%	162.0 / NA	50	P725L (CCG→CTG)	acrB	multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/24);  total (26/24)

TGACATAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGC  >  NC_000913_3_hsa_tpiA/482084‑482375                                                                                                                                                   |                                                                                                                                                  tgacatAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCagg                                                                                                                                                 >  1:526349/1‑149 (MQ=255)         aCTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTc                                                                                                                                         <  1:857313/149‑1 (MQ=255)         aCTTTCTGCGCACGACCGCGGGCGCATAAGTCGTTGACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATGGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTc                                                                                                                                         <  2:140320/149‑1 (MQ=255)           tttCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGCAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCa                                                                                                                                       <  1:775959/149‑1 (MQ=255)                  acacGACCGCGGTCGATAAAGTCGTTCACATAGCTGCAGCCCCATGCAGCGCCCAGAGTGGCGCCAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGtt                                                                                                                                >  1:811306/1‑149 (MQ=255)                       aCCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGAc                                                                                                                           <  2:669378/149‑1 (MQ=255)                           cgGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTAc                                                                                                                       >  1:546720/1‑149 (MQ=255)                              tCGATAACGTCGTTCACATAGCTGCCGCCCCATGCAGCCCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCACTTTTTCCTGGTCGAGATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTAGCCATACGCt                                                                                                                    >  2:202665/1‑149 (MQ=255)                                    aaGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCaa                                                                                                              <  1:858617/149‑1 (MQ=255)                                     cgTCATTCACATGGCAACCGCCCCATGCAGCTCCCAGAGTGTTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCGGGTCGACATCAATCTCAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAAc                                                                                                             <  1:465973/148‑1 (MQ=255)                                          ttCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATc                                                                                                        <  2:526349/149‑1 (MQ=255)                                                 aGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGc                                                                                                 >  1:352966/1‑149 (MQ=255)                                                  gCTGCCGCCCCAGGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCt                                                                                                <  2:577287/149‑1 (MQ=255)                                                      ccgccCCAGGCAGCGGCCGGGGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCgc                                                                                            <  2:352966/149‑1 (MQ=255)                                                       cgccCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCgct                                                                                           <  1:288821/149‑1 (MQ=255)                                                        gccCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTTCTGGTTTTTATCAATCTTAAACTGTAGGGTATCTTCCAAACCGTTTGGACGTTCGCTTGTCACCATGTCAGGGTGCTTCgctg                                                                                          >  2:803248/1‑149 (MQ=255)                                                          cccATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGa                                                                                                                            >  2:605510/1‑113 (MQ=255)                                                          cccATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGa                                                                                                                            <  1:605510/113‑1 (MQ=255)                                                            cATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCGGACCGGTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTc                                                                                      >  2:761405/1‑149 (MQ=255)                                                              tGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTg                                                                                    <  1:555837/149‑1 (MQ=255)                                                                 aGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTCTCTTCCACACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCCCTGCTTCTtcaa                                                                                 >  2:674003/1‑149 (MQ=255)                                                                 aGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTACTGGTCGATATCCATATTCAACTGCAGGGTATCTTGCAAACCACTTGCTCGTACGCTGGCAAGCTTCTCTGGTTGCTTCGCTGCTTATgcga                                                                                 >  2:389673/1‑147 (MQ=255)                                                                              gtgtTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCACCTGGTTAcg                                                                    >  1:259632/1‑149 (MQ=255)                                                                               tgtTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGtt                                                                       >  1:782189/1‑145 (MQ=255)                                                                               tgtTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGtt                                                                       <  2:782189/145‑1 (MQ=255)                                                                                   aaTGTAGTGGATAGAAACACGCAGCGCCTGCTCTGTTTCCTCGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCt                                                               <  2:299637/149‑1 (MQ=255)                                                                                         gTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCa                                                         <  1:254748/149‑1 (MQ=255)                                                                                                gAAACACCCAGCGCCTGCTGTTTTTCCTGGTCGATATCATCCCTAATCGGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTc                                                  <  2:309993/149‑1 (MQ=255)                                                                                                gAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCCATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTc                                                  <  1:775243/149‑1 (MQ=255)                                                                                                 aaaCACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCGAGACCGTTTGGAAGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAGACAACTGTTTACGCGCCTGGGTCAGTCTTTCg                                                 >  1:210532/1‑149 (MQ=255)                                                                                                 aaaCACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGAAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTg                                                                                                  >  2:204977/1‑100 (MQ=255)                                                                                                 aaaCACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACCGAAGGGAATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTg                                                                                                  <  1:204977/100‑1 (MQ=255)                                                                                                  aaCACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGt                                                <  2:811306/149‑1 (MQ=255)                                                                                                      cccAGCGCATGCGCTTTTTACTGGTCGATATCAATCTTAAACTGCAGAGTATCTTCCAGACCGTTTGGACGTACGCTAATCAACATATAAGGGTACTTCGCTGCTTCTCCAATCAACTGATTACTCACCTGATACAGTTTTTCATGaac                                            >  2:84823/1‑147 (MQ=255)                                                                                                            gcCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCaa                                          <  2:59335/145‑1 (MQ=255)                                                                                                            gcCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCaa                                          >  1:59335/1‑145 (MQ=255)                                                                                                            gcCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGcc                                      <  2:546720/149‑1 (MQ=255)                                                                                                               tGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCC‑GAGGCAGTTTTTCGTGACCAATGCCAGcc                                  >  2:276615/1‑149 (MQ=255)                                                                                                                    tttttCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAGCTAGTTACGCGCCTGCGTCAGTTTTTCTTTACCACGGACAGccggat                              >  2:76140/1‑145 (MQ=255)                                                                                                                           tGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGATTACGCGCCTGAGTCAGTTTTTCGTGACCAAGCCCAGCCTGGTCAAGCAg                       >  2:227695/1‑149 (MQ=255)                                                                                                                             gTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCt                     >  1:576630/1‑149 (MQ=255)                                                                                                                              tCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtc                    <  1:520190/149‑1 (MQ=255)                                                                                                                                gATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCCGAGTCAGTTTTACGTGACCAAGGCCAGCCTGGCCAATCAGCccaa                  >  2:263062/1‑149 (MQ=255)                                                                                                                                   atCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCACGTATCCCGTGTCGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagt               >  1:201551/1‑149 (MQ=255)                                                                                                                                      aaTCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTgg                                                                         <  2:256490/88‑1 (MQ=255)                                                                                                                                      aaTCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTgg                                                                         >  1:256490/1‑88 (MQ=255)                                                                                                                                        tCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaacggcaaag          >  1:213294/1‑149 (MQ=255)                                                                                                                                         cTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCAAGCCTGGTCAATCAGCTCAAAGTCAAAGc         >  1:218889/1‑149 (MQ=255)                                                                                                                                               cTGCAGCGTATCCTCCAGCCCGTTTGGTCGGACGCGGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTg   <  2:757381/149‑1 (MQ=255)                                                                                                                                                tGCAGGGTATCTTCCAGACGGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTAGATACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTGAAGCGGCTCAAAGTCAAAGCCGGTTga  >  2:731272/1‑148 (MQ=255)                                                                                                                                                tGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGGCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCGGCTCAAAGTCAAAGCCGCTTGc  >  2:303604/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                  TGACATAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGC  >  NC_000913_3_hsa_tpiA/482084‑482375

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