Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_hsa_tpiA	1,755,785	G→T	V30F (GTT→TTT)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_hsa_tpiA	1,755,785	0	G	T	94.3%	224.1 / 6.1	70	V30F (GTT→TTT)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base T (30/36);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.69e-01

CAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGG  >  NC_000913_3_hsa_tpiA/1755640‑1755929                                                                                                                                                  |                                                                                                                                                 cAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTAt                                                                                                                                              <  3:264994/150‑1 (MQ=255)  aGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATg                                                                                                                                             <  1:564529/150‑1 (MQ=255)       gcCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAGTCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCt                                                                                                                                        <  4:464103/150‑1 (MQ=255)       gcCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTCAAATGCTGGACGCTGGCATGAACTTTATGCGTCt                                                                                                                                        <  4:388926/150‑1 (MQ=255)       gcCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCt                                                                                                                                        <  1:643472/150‑1 (MQ=255)          cAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGaa                                                                                                                                     <  2:527250/150‑1 (MQ=255)                gCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGACATGAACTTTATGCGTCTGAACTtctc                                                                                                                               <  1:577951/150‑1 (MQ=255)                 cAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTtctct                                                                                                                              <  3:130351/150‑1 (MQ=255)                  aaGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACtt                                                                                                                                  >  4:71619/1‑145 (MQ=255)                  aaGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACtt                                                                                                                                  <  3:71619/145‑1 (MQ=255)                       tctcCCATCCTTCTCACCTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATCCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGt                                                                                                                        <  2:750282/150‑1 (MQ=255)                         tcCCATCCTTCTCCACTGAAAGACCACCACTGTCATGACAAAGACTAAAACTGTTTCCACCCTCGGCGCGAAAACCGAATCTCAAGAGATGTTAGCTAAAATGCTGGACGCTCGCATGAACTTTATCCGTCTGAACTTCTCTCATGGtca                                                                                                                      <  2:200905/150‑3 (MQ=255)                           ccATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACt                                                                                                                    >  1:259903/1‑150 (MQ=255)                             aTCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGc                                                                                                                                            <  2:488551/124‑1 (MQ=255)                             aTCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGc                                                                                                                                            >  1:488551/1‑124 (MQ=255)                              tCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTtc                                                                                                                                 <  1:216904/134‑1 (MQ=255)                              tCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTtc                                                                                                                                 >  2:216909/1‑134 (MQ=255)                              tCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTtc                                                                                                                                 >  2:216904/1‑134 (MQ=255)                              tCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTtc                                                                                                                                 <  1:216909/134‑1 (MQ=255)                                   ctcAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGaa                                                                                                            <  3:450086/150‑1 (MQ=255)                                    tcAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAac                                                                                                           >  1:365106/1‑150 (MQ=255)                                     cAACTTAAAGACTAAGACTGTCATGAAAAAGACCCAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAaca                                                                                                          >  2:533331/1‑150 (MQ=255)                                     cAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAGCTTCTCTCATGGTGACTATGCAGAaca                                                                                                          <  1:577339/150‑1 (MQ=255)                                      aaCTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAacac                                                                                                         >  1:536732/1‑150 (MQ=255)                                         ttAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGt                                                                                                      >  2:715091/1‑150 (MQ=255)                                               actaagactGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAgcgc                                                                                                >  4:465973/1‑150 (MQ=255)                                                 taagactGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCAt                                                                                              >  2:645712/1‑150 (MQ=255)                                                 taagactGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCAt                                                                                              <  3:98842/150‑1 (MQ=255)                                                   agactGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTc                                                                                            <  4:182411/150‑1 (MQ=255)                                                        gTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAAt                                                                                       >  3:136231/1‑150 (MQ=255)                                                        gTCAAGAAAAAGACCAAAATAGTTTGCACCACCGGACCGAAAACCGAATATGAAGAGAAGTTAGCTAAAATGTTGGACGCTGACATCAACTCTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGATCAACGCATTCAGatt                                                                                       >  2:387326/1‑148 (MQ=255)                                                           aTGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCAAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGaa                                                                                                            >  1:712448/1‑126 (MQ=255)                                                           aTGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCAAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGaa                                                                                                            <  2:712448/126‑1 (MQ=255)                                                                                   aCCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGt                                                            <  1:715091/150‑1 (MQ=255)                                                                                       gctgACGGAAAACCGGCTCTGAAGAGATGTTAGCTAAGAGGCTGGAAGCTGGCATGAACTTTATGCGTCTGAACTGCTCGCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTagaa                                                        <  4:284099/147‑1 (MQ=255)                                                                                         ggACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAAcc                                                      >  1:653894/1‑150 (MQ=255)                                                                                           aCCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAAccgc                                                    <  2:653894/150‑1 (MQ=255)                                                                                             cGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAAccgccg                                                  >  1:416479/1‑150 (MQ=255)                                                                                             cGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAAccgccg                                                  <  1:146221/150‑1 (MQ=255)                                                                                              gAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGTGCAACGTGATGAGCAAAACTGGTAAAACcgccgc                                                 >  2:682448/1‑150 (MQ=255)                                                                                               aaaaCCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCt                                                >  4:307756/1‑150 (MQ=255)                                                                                                  aCCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATc                                             <  4:136231/150‑1 (MQ=255)                                                                                                   ccGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATcc                                            >  4:179998/1‑150 (MQ=255)                                                                                                   ccGAATCTGAAGAGATGTTAGCTAAAAGGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATcc                                            >  1:577467/1‑150 (MQ=255)                                                                                                      aaTCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGc                                                                                                                <  1:390748/79‑1 (MQ=255)                                                                                                      aaTCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGc                                                                                                                >  2:390748/1‑79 (MQ=255)                                                                                                       aTCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCt                                        <  3:179998/150‑1 (MQ=255)                                                                                                        tCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGCACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCtt                                       <  3:360816/150‑1 (MQ=255)                                                                                                         cTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTg                                      >  2:91558/1‑150 (MQ=255)                                                                                                          tGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGa                                     <  1:643796/150‑1 (MQ=255)                                                                                                           gAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGAt                                    >  2:731811/1‑150 (MQ=255)                                                                                                           gAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGAt                                    <  4:35949/150‑1 (MQ=255)                                                                                                            aaGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATa                                   >  2:270446/1‑150 (MQ=255)                                                                                                                 aTGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAAccgccg                                                  <  1:334357/130‑1 (MQ=255)                                                                                                                 aTGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAAccgccg                                                  >  2:334357/1‑130 (MQ=255)                                                                                                                     tAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGtgatga                                                                       <  2:75816/105‑1 (MQ=255)                                                                                                                     tAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGggatga                                                                       >  1:75816/1‑105 (MQ=255)                                                                                                                       gCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCg                        <  2:577467/150‑1 (MQ=255)                                                                                                                       gCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTGCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCg                        <  2:482401/150‑1 (MQ=255)                                                                                                                            aaTGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAAt                   >  2:22046/1‑150 (MQ=255)                                                                                                                            aaTGCTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAAt                   <  1:645712/150‑1 (MQ=255)                                                                                                                                cTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGc               <  3:465973/150‑1 (MQ=255)                                                                                                                                cTGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGc               <  1:91558/150‑1 (MQ=255)                                                                                                                                 tGGATGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCa              >  2:537941/1‑150 (MQ=255)                                                                                                                                 tGGACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCa              >  4:217827/1‑150 (MQ=255)                                                                                                                                   gACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGc                                                                                                                <  4:21503/50‑1 (MQ=255)                                                                                                                                   gACGCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGc                                                                                                                >  3:21503/1‑50 (MQ=255)                                                                                                                                      gCTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATg         >  2:63791/1‑150 (MQ=255)                                                                                                                                       cTGGCATGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGa        <  4:346072/150‑1 (MQ=255)                                                                                                                                             tGAACTTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTgg  >  4:147013/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                 CAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGG  >  NC_000913_3_hsa_tpiA/1755640‑1755929

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