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%	139.6 / NA	44	P725L (CCG→CTG)	acrB	multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/19);  total (25/19)

TGACATAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAG  >  NC_000913_3_hsa_tpiA/482084‑482367                                                                                                                                                   |                                                                                                                                          tgacatAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCagg                                                                                                                                         <  1:559487/149‑1 (MQ=255) tgacatAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGATCACATAGATGCCGCCCAATGCAGCGCCCAGAGTGGTGTTAATGTAGTTGGTACGAGGTGACATCGACAGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCagg                                                                                                                                         >  1:21554/1‑149 (MQ=255)   acatAGACTTGCATCCCACGCCCGCGGTCGATAAAGTCGTTTACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGGCGTTGATAGAAACACCCAGCGCCTGCGCTTTCTCCTGGTCGATATCAATCTTAAACTGCAGGGt                                                                                                                                       <  2:47079/149‑1 (MQ=255)         aCTTTCTGCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTc                                                                                                                                 <  1:316754/149‑1 (MQ=255)          cTTTCTTCACACGACCGCGGTGGATAAGGTCGTTCACATAGCTGCCGTCCTATGCAGCGCGCAGGGTGGTGTTAATGTCGTTGATAGTAACCCCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATTTTCAACTGCAGAGTATCttgc                                                                                                                                >  2:666698/1‑147 (MQ=255)           tttCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCa                                                                                                                               <  2:63566/149‑1 (MQ=255)               ttcACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGAcc                                                                                                                           >  1:625189/1‑149 (MQ=255)                           cagTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATCACGTGCCTGTGCGGCAGGGTATCTTCCAGACCGTTTGGACGTAc                                                                                                               <  2:697516/147‑1 (MQ=255)                             gTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGc                                                                                                             >  1:15430/1‑149 (MQ=255)                                  tAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTGTCTTGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTc                                                                                                        >  1:365144/1‑149 (MQ=255)                                       tCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGGAGGCTATCTTCCAGACCGTTTGGACGTACGCTGGTCCACat                                                                                                   >  2:254430/1‑149 (MQ=255)                                         gTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACatat                                                                                                 >  2:711886/1‑149 (MQ=255)                                         gTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACatat                                                                                                 >  2:711942/1‑149 (MQ=255)                                             acaTAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGGTAATGTCGTTGAGAGAAACACCCAGCACCTGAGCTTTTTCCTGGTCGATATCAATCTTAAACTGAAGAGTATATTCAAGACCGTTTGGACGTACGCGGGTAAACatattaag                                                                                             >  2:457226/1‑145 (MQ=255)                                               aTAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTCGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCCGGTCAACCTCTCCGGGt                                                                                           >  1:27432/1‑149 (MQ=255)                                                      ccgccCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCgc                                                                                    <  1:711942/149‑1 (MQ=255)                                                      ccgccCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTGAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCgc                                                                                    <  1:711886/149‑1 (MQ=255)                                                           ccATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCtt                                                                               >  1:628349/1‑149 (MQ=255)                                                             aTGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGAGCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCt                                                                             <  2:694858/149‑1 (MQ=255)                                                             aTGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCt                                                                             <  1:751256/149‑1 (MQ=255)                                                                  gcgcCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTgcaag                                                                        >  1:421631/1‑149 (MQ=255)                                                                        agagTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTg                                                                  <  1:544564/149‑1 (MQ=255)                                                                                gtTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCCACATATCAGGGTGCATCGCTGCTTCTCCAAGCAACTGGTTAcaca                                                          >  2:557569/1‑146 (MQ=255)                                                                                 ttAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACCCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAACCAACTGGGTACgcgc                                                         >  1:656045/1‑149 (MQ=255)                                                                                   aaTGTCGTTGATAGCAATACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCt                                                       >  2:279405/1‑149 (MQ=255)                                                                                   aaTGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGGCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCt                                                       <  1:622116/149‑1 (MQ=255)                                                                                         gTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCGGCAAGCAACTCGTTAAGCGCCTGCGTCa                                                 >  2:426502/1‑149 (MQ=255)                                                                                             aTAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATc                                                                                                <  1:689933/98‑1 (MQ=255)                                                                                             aTAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTACAGACCGTTTGGACGTACGCTGGTCAACATATc                                                                                                >  2:689933/1‑98 (MQ=255)                                                                                                gAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGCCCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCCGAGTCAGTTTATc                                          >  1:23392/1‑149 (MQ=255)                                                                                                caaaCACCCAGCGACTGCGCTTTTTCCTGGTCGGTATCAAGCTTAAACTTCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGAGGCTTCGCTGCTTCTGCAAGCAGCTGGTTACCGGCCTGGTTCAGTTTTTc                                          >  2:258599/2‑149 (MQ=255)                                                                                                       ccAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATGAGGGGGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCa                                   >  1:103945/1‑149 (MQ=255)                                                                                                         aGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATTTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAg                                 <  2:330967/149‑1 (MQ=255)                                                                                                          gcgcCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAgg                                <  2:551750/149‑1 (MQ=255)                                                                                                              cTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAg                            >  2:710206/1‑149 (MQ=255)                                                                                                               tGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGc                           >  1:612135/1‑149 (MQ=255)                                                                                                                 cccTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCt                         <  1:115331/147‑1 (MQ=255)                                                                                                                  gcTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTg                        <  1:279405/149‑1 (MQ=255)                                                                                                                  gcTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGAAACTGGTTACGCACCTGAGTCAGATTTTAGTGACAACAGcctgcctg                        >  2:25151/1‑149 (MQ=255)                                                                                                                       ttCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCaa                   <  2:628349/149‑1 (MQ=255)                                                                                                                          cTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCa                <  1:187560/149‑1 (MQ=255)                                                                                                                             gTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCt             >  2:710016/1‑149 (MQ=255)                                                                                                                                   atCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagt       <  1:710206/149‑1 (MQ=255)                                                                                                                                       aTCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtcaaa   <  2:612135/149‑1 (MQ=255)                                                                                                                                        tCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtcaaag  >  2:373697/1‑149 (MQ=255)                                                                                                                                        tCTTAAACTGCAGGGTATCTTCCAGACCGTGTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCACCTGAATCACTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaaatcaaag  >  2:230192/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                          TGACATAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAG  >  NC_000913_3_hsa_tpiA/482084‑482367

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