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	98.3%	184.3 / NA	58	P725L (CCG→CTG)	acrB	multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (28/29);  minor base C (0/1);  total (28/30)
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

TTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAC‑GCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGC  >  NC_000913_3_hsa_tpiA/482094‑482375                                                                                                                                         |                                                                                                                                                   tttCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCTTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCa                                                                                                                                        >  2:390138/1‑149 (MQ=255)   tcttcAGCCGACCGCGTGCGATAAGTTCGTTCACATAGTGGCCGCCCCATGCTGCTCCCAGAGTGGTGTTAATGTCGTGGATCGAAACACCCAGCGCCTGCTCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGa                                                                                                                                      <  2:652784/149‑1 (MQ=255)     ttcACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGAcc                                                                                                                                    <  1:522314/149‑1 (MQ=255)                  ggTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACg                                                                                                                       >  1:392356/1‑149 (MQ=255)                    tCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCt                                                                                                                     <  2:75857/149‑1 (MQ=255)                     cGATAAATTCGTCCACAAAGCTTCCTCCCCCCGGAGGCCCTATATTTGTCTTAATGTCGTTGATAGAAACCCCCAGCGCCTGCGCTTTTTCCTGGTCGGTATTAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTg                                                                                                                    <  2:44317/149‑1 (MQ=255)                     cGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTg                                                                                                                    <  2:543667/149‑1 (MQ=255)                     cGATAAAGTCGTTCACATAGCTGCCGCCCCACGCAGCGCCCAGCTTTTTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTg                                                                                                                    <  2:599704/149‑1 (MQ=255)                                ttCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATc                                                                                                         >  2:564378/1‑149 (MQ=255)                                   acaTAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAgg                                                                                                      >  2:703053/1‑149 (MQ=255)                                       aGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCCGCGCCTGCGCTTTTCCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGCCCGTTTGGACGTACGCTGGTCAACATATCAGGGTGc                                                                                                  <  1:487771/149‑1 (MQ=255)                                        gCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACGTATCAGGGTGCt                                                                                                 >  1:322949/1‑149 (MQ=255)                                          tGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTc                                                                                               >  2:534644/1‑149 (MQ=255)                                           gccgccCCATGCGGCGCCAAGAGTGCAGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCg                                                                                              <  2:635442/149‑1 (MQ=255)                                                  cATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATc                                                                                                         <  1:980/131‑1 (MQ=255)                                                  cATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATc                                                                                                         >  2:980/1‑131 (MQ=255)                                                   aggCAGCGCCCAGAGTGGTGTGAATGTCGTTGATAGAAGCACCCAGCGCCTGCGCTTTGTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCt                                                                                      <  1:802232/147‑1 (MQ=255)                                                       aGCGCCCAGAGTTGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAATCTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTgcaa                                                                                  <  1:353173/149‑1 (MQ=255)                                                            ccAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAAc                                                                             <  2:307366/149‑1 (MQ=255)                                                                  tggtgTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGt                                                                                                                  >  2:287647/1‑106 (MQ=255)                                                                  tggtgTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGt                                                                                                                  <  1:287647/106‑1 (MQ=255)                                                                  tggtgTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTa                                                                       >  1:342074/1‑149 (MQ=255)                                                                  tggtgTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTa                                                                       >  1:757359/1‑149 (MQ=255)                                                                       ttAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAC‑gcgc                                                                 <  2:394821/149‑1 (MQ=255)                                                                         aaTGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTAGCCTGTTCAACATGTCAGGGAGCTGCGCTGCTTCTGCAAGCGACTGGTTAC‑GCGCCt                                                               >  1:617866/1‑149 (MQ=255)                                                                             tCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAC‑GCGCCTGAGt                                                           >  1:73185/1‑149 (MQ=255)                                                                                      gAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAC‑GCGCCTGAGTCAGTTTTTc                                                  >  2:371628/1‑149 (MQ=255)                                                                                      gAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGGGCTTCGCTGCTTCTGCAAGCAACTGGTTAC‑GCGCCTGAGTCAGTTTTTc                                                  >  1:500920/1‑149 (MQ=255)                                                                                       aaaCACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAC‑GCGCCTGAGTCAGTTTTTCg                                                 >  1:351240/1‑149 (MQ=255)                                                                                          cacCCATCGCCTGCGCTCTTTCCTGGTAGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCTTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAC‑GCGCCTGAGTCAGTTTTTCGTGa                                              <  1:390138/149‑1 (MQ=255)                                                                                           acCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAC‑GCGCCTGCGGCAGTTTTTGGTGAc                                             >  2:555851/1‑149 (MQ=255)                                                                                             ccAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAc                                                                      >  1:734698/1‑123 (MQ=255)                                                                                             ccAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAc                                                                      <  2:734698/123‑1 (MQ=255)                                                                                             ccAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAC‑GCGCCTGAGTCAGTTTTTCGTGACCa                                           >  1:364197/1‑149 (MQ=255)                                                                                                        gcTTTTCCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAC‑GCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTg                                <  1:92729/149‑1 (MQ=255)                                                                                                          tttttCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAC‑GCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGt                              <  2:392356/149‑1 (MQ=255)                                                                                                              tCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAC‑GCGGCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                          <  2:757359/149‑1 (MQ=255)                                                                                                              tCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAC‑GCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                          <  2:342074/149‑1 (MQ=255)                                                                                                              tCCTGGGCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAC‑GCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                          <  2:107056/149‑1 (MQ=255)                                                                                                                 tGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAA‑CAACCCCTGGCGGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGGCAATCAg                       >  1:583295/1‑149 (MQ=255)                                                                                                                  gctCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAC‑GCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGc                      <  1:59379/147‑1 (MQ=255)                                                                                                                   gTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCtt                                                                                        >  1:422868/1‑83 (MQ=255)                                                                                                                   gTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCtt                                                                                        <  2:422868/83‑1 (MQ=255)                                                                                                                   gTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAC‑GCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCt                     >  2:780668/1‑149 (MQ=255)                                                                                                                    tCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAC‑GCGCCTGa                                                             >  2:130843/1‑108 (MQ=255)                                                                                                                    tCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAC‑GCGCCTGa                                                             <  1:130843/108‑1 (MQ=255)                                                                                                                      gATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGt                                                                                                                  >  1:137466/1‑54 (MQ=255)                                                                                                                      gATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGt                                                                                                                  <  2:137466/54‑1 (MQ=255)                                                                                                                       atatCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAC‑GCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaa                 <  1:534644/149‑1 (MQ=255)                                                                                                                              tCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAC‑GCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtcaaag          >  1:653003/1‑149 (MQ=255)                                                                                                                                 cacaCTGCCGGGTATCTTCCAGTCCGTTTGGACGTACGGTGTTCAACATATCAGGGTGTTTCGCTGCTTCTGCAAGTAACTGGTTAC‑GCCCCTGAGTCAGTTTTTCTTCACCAAGGCCAGCCTGCTCAATCAGCTCAAAGTCAAAGCCg       <  2:241587/146‑1 (MQ=255)                                                                                                                                   aaCTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAC‑GCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTg                                <  1:87347/122‑1 (MQ=255)                                                                                                                                   aaCTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAC‑GCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTg                                >  2:87347/1‑122 (MQ=255)                                                                                                                                    ccTGCAGGGTATCGTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAC‑GCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCACTCAGCTCAAAGTCAAAGCCGGtt    <  2:84712/148‑1 (MQ=255)                                                                                                                                     cTGCAGGGTATCTTCCAGACCGTTTGGACGGACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAC‑GCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTg   <  1:770067/149‑1 (MQ=255)                                                                                                                                      tGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAC‑GCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGc  >  1:604692/1‑149 (MQ=255)                                                                                                                                      tGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAC‑GCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGc  >  1:698337/1‑149 (MQ=255)                                                                                                                                      tGCAGGGTATCATCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAC‑GCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAACGTCAACTCCGGTTGc  >  1:532345/1‑149 (MQ=255)                                                                                                                                         |                                                                                                                                                   TTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAC‑GCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGC  >  NC_000913_3_hsa_tpiA/482094‑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: ‑