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

TGACATAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCC‑TGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAG  >  NC_000913_3_hsa_tpiA/482084‑482367                                                                                                                                                   |                                                                                                                                           tgacatAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCagg                                                                                                                                          >  2:355700/1‑149 (MQ=255)     atAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTAt                                                                                                                                      >  2:296303/1‑149 (MQ=255)           tttCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCa                                                                                                                                <  2:360146/149‑1 (MQ=255)                tcACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCg                                                                                                                           >  2:481724/1‑149 (MQ=255)                     cGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAAATGCAGGGTATCTTCCAGACCGTTTgg                                                                                                                      >  1:10961/1‑149 (MQ=255)                              tCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCt                                                                                                             <  1:355700/149‑1 (MQ=255)                                gATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGGTAATGTGGTTGATAGAAACACCGAGCGCCGGCGCTTTTTCCTGAGGGGCGCCGATGTGAGACTGTAGGGTTTCTTCGGGACCGTTTGGACGTCCGCAgg                                                                                                           >  1:143519/1‑149 (MQ=255)                                                tAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGCGTg                                                                                           >  2:382088/1‑149 (MQ=255)                                                                     cccAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTgcaagcaa                                                                      <  1:296303/149‑1 (MQ=255)                                                                      ccAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAAc                                                                     >  2:503424/1‑149 (MQ=255)                                                                      ccAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAAc                                                                     >  1:154388/1‑149 (MQ=255)                                                                                 ttAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACgcgc                                                          <  1:481724/149‑1 (MQ=255)                                                                                 ttAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACgcgc                                                          >  1:267725/1‑149 (MQ=255)                                                                                   aaTGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCC‑t                                                       <  1:503424/149‑1 (MQ=255)                                                                                          ttGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCC‑TGAGTCAg                                                >  1:282594/1‑149 (MQ=255)                                                                                           tGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCC‑TGAGTCAGt                                               >  1:287869/1‑149 (MQ=255)                                                                                               aGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCC‑TGAGTCAGttttt                                           <  1:387038/149‑1 (MQ=255)                                                                                                gAAACACCCAGCGCCTTCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCC‑TGAGTCAGTTTTTc                                          <  2:10961/149‑1 (MQ=255)                                                                                                    cagccAGCGCCCGCGCTTTTGCCTGGTCGATAACAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGGACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCC‑TGAGTCAGTTTTTCGTGa                                      <  2:497634/146‑1 (MQ=255)                                                                                                     acCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTGCACATATCAGGGTGCTTCCCTGCTTCTGCAACCAACTCGTTCACCGCC‑TGACTCATTTTTTCGTGAc                                     >  2:12288/1‑149 (MQ=255)                                                                                                        cAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACAGGTTACGCGCC‑TGAGTCAGTTTTTCGTGACCaa                                  >  2:396736/1‑149 (MQ=255)                                                                                                         aGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCC‑TGAGTCAGTTTTTCGTGACCAAg                                 <  1:380107/149‑1 (MQ=255)                                                                                                                    tttttCCTGGTCGATATCAATCTTAAACTGCAGGGTCTCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCC‑TGAGTCAGTTTTTCGTGGCCAAGGCCAGCCTGGt                      >  1:500926/1‑149 (MQ=255)                                                                                                                           tGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAA‑ACGCCACGAACCGTTACTCCGCTTTACGTTACCAAGTCCAGCCTGGTCAATCAg               >  1:245671/1‑149 (MQ=255)                                                                                                                              tCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCC‑TGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtc            >  1:318263/1‑149 (MQ=255)                                                                                                                                  tatCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCC‑TGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaag        >  2:437725/1‑149 (MQ=255)                                                                                                                                    tCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCgc                                                                                     >  2:10373/1‑71 (MQ=255)                                                                                                                                    tCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCgc                                                                                     <  1:10373/71‑1 (MQ=255)                                                                                                                                    tCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCC‑TGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGCCAATCAGCtcaaattc      >  2:227609/1‑149 (MQ=255)                                                                                                                                        tCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACATACGCTGGTCAACAAATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCC‑TGAGTCAGTTTTTCGTTACCAAGGCCAGCCTGGTCACTGAGCtcaaagcccacg  >  2:493021/1‑143 (MQ=255)                                                                                                                                                   |                                                                                                                                           TGACATAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCC‑TGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAG  >  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: ‑