Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_hsa_tpiA	1,757,058	G→C	G454A (GGT→GCT)	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,757,058	0	G	C	92.2%	208.6 / 10.0	64	G454A (GGT→GCT)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base G (4/1);  new base C (28/31);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

GCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAA  >  NC_000913_3_hsa_tpiA/1756913‑1757200                                                                                                                                                  |                                                                                                                                               gCTCATCAGTTGGTATTGAGCAAAGGCCTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCa                                                                                                                                            <  2:617880/150‑1 (MQ=255)  cTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCAc                                                                                                                                           >  1:377364/1‑150 (MQ=255)   tcatcaGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACt                                                                                                                                          >  4:419084/1‑150 (MQ=255)           tGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGa                                                                                                                                  <  4:226159/150‑1 (MQ=255)              tACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCg                                                                                                                               >  2:197577/1‑150 (MQ=255)                 tGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCa                                                                                                                            >  2:564352/1‑150 (MQ=255)                   aGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCg                                                                                                                               <  2:37186/145‑1 (MQ=255)                   aGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCg                                                                                                                               >  1:37186/1‑145 (MQ=255)                    gCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCacta                                                                                                                         >  2:222325/1‑150 (MQ=255)                     cAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCactac                                                                                                                        >  3:406217/1‑150 (MQ=255)                     cAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCactac                                                                                                                        <  4:306218/150‑1 (MQ=255)                      aaaGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCactact                                                                                                                       <  1:663245/150‑1 (MQ=255)                      aaaGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTACAGAGCGATCTTACACACAAAGGTGACGTTGTAGTTATAGTTTCTGCTGCACTGGGACCGAGGTGCactact                                                                                                                       <  2:424143/150‑1 (MQ=255)                            gTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACAcc                                                                                                                 <  1:101894/150‑1 (MQ=255)                            gTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACAcc                                                                                                                 <  1:758744/150‑1 (MQ=255)                                tgCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCAt                                                                                                             <  1:197577/150‑1 (MQ=255)                                          tGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGt                                                                                                   >  1:547779/1‑150 (MQ=255)                                           gTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTAACGTc                                                                                                  >  4:473140/1‑150 (MQ=255)                                                agagATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACAcc                                                                                                                 <  1:816673/130‑1 (MQ=255)                                                agagATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACAcc                                                                                                                 >  2:816673/1‑130 (MQ=255)                                                agagATCACTTCTACTCATGATTTCTACCGTCTGGGTAAAGACCTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTGGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTa                                                                                             <  3:2655/150‑1 (MQ=255)                                                  agaTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAat                                                                                           >  2:658848/1‑150 (MQ=255)                                                       aCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGc                                                                                      <  4:358729/150‑1 (MQ=255)                                                             aCTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTtgt                                                                                <  3:299020/150‑1 (MQ=255)                                                               tgatgaTTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGa                                                                              >  2:206650/1‑150 (MQ=255)                                                               tgatgaTTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGa                                                                              <  1:222325/150‑1 (MQ=255)                                                                 atgatTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGaat                                                                            >  2:151090/1‑150 (MQ=255)                                                                 atgatTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGaat                                                                            >  2:821330/1‑150 (MQ=255)                                                                       tCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCt                                                                                                <  4:234972/124‑1 (MQ=255)                                                                       tCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCt                                                                                                >  3:234972/1‑124 (MQ=255)                                                                         tACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGt                                                                    <  4:406217/150‑1 (MQ=255)                                                                         tACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGt                                                                    <  3:63377/150‑1 (MQ=255)                                                                             gTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTatat                                                                >  4:205782/1‑150 (MQ=255)                                                                                 gggTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGaa                                                            <  2:602372/150‑1 (MQ=255)                                                                                  ggTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATAtt                                                                                        <  1:174152/121‑1 (MQ=255)                                                                                  ggTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATAtt                                                                                        >  2:174152/1‑121 (MQ=255)                                                                                      aaGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGcc                                                       >  4:188925/1‑150 (MQ=255)                                                                                       aGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGccc                                                      <  1:302977/150‑1 (MQ=255)                                                                                           cTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGAt                                                  <  2:340399/150‑1 (MQ=255)                                                                                             ggCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATgg                                                <  1:658848/150‑1 (MQ=255)                                                                                                 ctGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTCATTTGTATATCGAAGCGcccatcaggg                                               <  1:618049/147‑8 (MQ=255)                                                                                                 ctGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGcccatcaggg                                               >  2:618049/1‑140 (MQ=255)                                                                                                      gagCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTa                                                                                             >  2:93759/1‑96 (MQ=255)                                                                                                      gagCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTa                                                                                             <  1:93759/96‑1 (MQ=255)                                                                                                      gagCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCttttt                                       <  2:676264/150‑1 (MQ=255)                                                                                                       agCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTtgt                                                                                >  1:179398/1‑108 (MQ=255)                                                                                                       agCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTtgt                                                                                <  2:179398/108‑1 (MQ=255)                                                                                                           gTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTAACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTtttatt                                  <  3:473140/150‑1 (MQ=255)                                                                                                            tCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTtttattt                                 >  2:615752/1‑150 (MQ=255)                                                                                                             cTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGt                                                                    <  1:715433/114‑1 (MQ=255)                                                                                                             cTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGt                                                                    >  2:715433/1‑114 (MQ=255)                                                                                                              tGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTaa                               >  3:308706/1‑150 (MQ=255)                                                                                                              tGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTaa                               >  3:308534/1‑150 (MQ=255)                                                                                                                gCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATc                             >  3:73874/1‑150 (MQ=255)                                                                                                                   cacaAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTTGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGAt                          <  2:377364/150‑1 (MQ=255)                                                                                                                    acaAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATa                         >  1:287789/1‑150 (MQ=255)                                                                                                                    acaAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATa                         >  2:156033/1‑150 (MQ=255)                                                                                                                            gACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAg                 <  1:206650/150‑1 (MQ=255)                                                                                                                                         tGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCttttt                                       >  3:495328/1‑115 (MQ=255)                                                                                                                                         tGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCttttt                                       <  4:495328/115‑1 (MQ=255)                                                                                                                                           gTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGCAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAaaatcaaa  >  2:88294/1‑150 (MQ=255)                                                                                                                                           gTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGgcg                                             <  2:194461/107‑1 (MQ=255)                                                                                                                                           gTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGgcg                                             >  1:194461/1‑107 (MQ=255)                                                                                                                                           gTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAaaatcaaa  >  3:117762/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                               GCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAA  >  NC_000913_3_hsa_tpiA/1756913‑1757200

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