Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_hsa_tpiA	4,588,333	C→T	R172H (CGC→CAC)	yjiA ←	metal‑binding GTPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_hsa_tpiA	4,588,333	0	C	T	88.9%	141.7 / 5.3	67	R172H (CGC→CAC)	yjiA	metal‑binding GTPase
Reads supporting (aligned to +/- strand):  ref base C (2/5);  new base T (26/30);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

TTTATCTGCGATAAAGTGGAAACGCGGTTTGGTGCTGACGACGTTTTCTTCCAGCATAAAACCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGCGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCAGAATGCGGTCGGCGTAGCCAACCTGTGACTGGGCGATGGTGAACTGGTTCATCTGCTCATCGGCATGTACCGCATCCACCAGCGCAATCACGCCGT  >  NC_000913_3_hsa_tpiA/4588188‑4588477                                                                                                                                                  |                                                                                                                                                 tttATCTGCGATAAAGTGGAAACGCGGTTTGGTGCTGACGACGTTTTCTTCCAGCATAAAACCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGtt                                                                                                                                               >  1:781542/1‑149 (MQ=255) tttATCTGCGATAAAGTGGAAACGCGGTTTGGTGCTGACGACGTTTTCTTCCAGCATAAAACCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGtt                                                                                                                                               >  2:534275/1‑149 (MQ=255)  gtATCTGAGATAAGGGGCAAACGCGGTTTGGGGCTGACGACGTTGTCTTCCAGCATAAAACCGTTGGTGTTGAACAGCAGGCCCAGTTCGATGGCGCCGTGGGGGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTc                                                                                                                                              <  1:464341/148‑1 (MQ=255)    aTCTGCGATAAAGTGGAAACGCGGTTTGGTGCTGACGACGTTTTCTTCCAGCATAAAACCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGCGTTCAt                                                                                                                                            <  1:389225/149‑1 (MQ=255)        gCGATAAAGTGGAAACGCGGTTTGGTGCTGACGACGTTTTCTTCCAGCATAAAACCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGCGTTCATGCAg                                                                                                                                        >  2:386241/1‑149 (MQ=255)            tAAAGTGGAAACGCGGTTTGGTGCTGACGACGTTTTCTTCCAGCATAAAACCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGGGCCAGGTGTTAATGCAGtttt                                                                                                                                    >  2:845461/1‑149 (MQ=255)            tAAAGTGGAAACGCGGTTTGGTGCTGACGACGTTTTCTTCCAGCATAAAACCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGtttt                                                                                                                                    <  2:651987/149‑1 (MQ=255)            tAAAGGGGAAACGCGGTGTGGTGCTGACGACGTTTTCTTCCAGCATAAAACCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGtttt                                                                                                                                    <  2:657212/149‑1 (MQ=255)             aaaGTGGAAACGCGGTTTGGTGCTGACGACGTTCTCTTCCAGCATAAAACCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGttttt                                                                                                                                   <  2:1060903/149‑1 (MQ=255)                gTGGAAACGCGGTTTGGTGCTGACGACGTTTTCTTCCAGCATAAAACCGTTGGTGTTGAAAAGCAGGCCCAGATCGATGTCGCCGTGGTTGACGGGGTAGACCGGTGCGCGGGGGTTGAGGCGGGCGAGGAGTTCATGCAGTTTTTCTg                                                                                                                                >  1:289029/1‑149 (MQ=255)                gTGGAAACGCGGTTTGGTGCTGACGACGTTGTCTTCCAGCATAAAACCGTTGGTGTTGAACAGCAGGCCCGGATCGATGTCTCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGAGAGCCAGGTGTTCCGGGAGTGTTTCTg                                                                                                                                >  1:592158/1‑149 (MQ=255)                 tGGAAACGCGGTTTGGTGCTGACGACGTTTTCTTCCAGCATAAAACCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGTTTTTCTGc                                                                                                                               >  2:384293/1‑149 (MQ=255)                   gAAACGCGGTTTGGTGCTGACGACGTTTTCTTCCAGCATAAAACCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGTTTTTCTGCtt                                                                                                                             <  1:676609/149‑1 (MQ=255)                      aCGCGGTTTGGTGCTGACGACGTTTTCTTCCAGCATAAAACCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGTTTTTCTGCTTCGc                                                                                                                          <  2:1061403/149‑1 (MQ=255)                      aCGCGGTTTGGTGCTGACGACGTTTTCTTCCAGCATAAAACCGTTGGTGTTGAAAAGCAGGCCCAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGTTTTTCGGCTTCGc                                                                                                                          >  1:864830/1‑149 (MQ=255)                          ggTTTGGTGCTGACGACGTTTTCTTCCAGCATAAAACCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGc                                                                                                                      >  1:677551/1‑149 (MQ=255)                           gTTTGGTGCTGACGACGTTTTCTTCCAGCATAAAACCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGACGGTGTAGACCTCGGCGCTGGCGGTATTGCGCGACAGGAAGTCATGCAGTTTTTCTTCTTCGCCTGCg                                                                                                                     >  2:289834/1‑149 (MQ=255)                               ggTGCTGACGACGTTTTCTTCCAGCATAAAACCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGt                                                                                                                 <  1:22391/149‑1 (MQ=255)                                  gCTGACGATGTTTTCTTCCAGCATAAAACCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCgg                                                                                                              <  1:1101619/149‑1 (MQ=255)                                   cTGACGACGTTTTCTTCCAGCATAAAACCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTCCGCCGCCGGt                                                                                                             >  2:13756/1‑149 (MQ=255)                                     gacgacGTTTTCTTCCAGCATAAAACCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGttt                                                                                                           >  2:903530/1‑149 (MQ=255)                                      acgacgTTTTCTTCCAGCATAAAACCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGtttt                                                                                                          >  2:257788/1‑149 (MQ=255)                                      acgacgTTTTCTTCCAGCATAAAACCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGtctt                                                                                                          >  2:1089690/1‑149 (MQ=255)                                                      gCATAAAACCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCAGAATGCgg                                                                                          >  1:774802/1‑149 (MQ=255)                                                        aTAAAACCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCAGAATGCGGTc                                                                                        <  1:801932/149‑1 (MQ=255)                                                            aaCCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTGATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCAGAATGCGGTCGGCg                                                                                    >  1:403071/1‑149 (MQ=255)                                                            aaCCGTTGGTGTTGAACAGCAGGCCCAGATCGATGGCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCAGAATGCGGTCGGCg                                                                                    >  2:1088139/1‑149 (MQ=255)                                                             aCCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCCGAATGCGGTCGGCGt                                                                                   >  2:281658/1‑149 (MQ=255)                                                              ccGTTGGTGTTTAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCAGAATGCGGTCGGCGTa                                                                                  <  1:1199687/149‑1 (MQ=255)                                                              ccGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCAGAATGCGGTCGGCGTa                                                                                  <  1:1199095/149‑1 (MQ=255)                                                              ccGTTGGTGGTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCAGAATGCGGTCGGCGTa                                                                                  <  1:62341/149‑1 (MQ=255)                                                                gTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGGGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCAGAATGCGGTCGGCGTAGc                                                                                >  1:1008039/1‑149 (MQ=255)                                                                      gtTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCAGAATGCGGTCGGGGTAGCCAACCt                                                                          >  2:572030/1‑149 (MQ=255)                                                                      gtTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCAGAATGCGGTCGGCGTAGCCAACCt                                                                          >  1:646199/1‑149 (MQ=255)                                                                          aaCGGCAGGCCCAGATCGTTGGCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGCGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGTGGAGCTTAATGCGGTCGGCGTAGCCAACCTGTGa                                                                      <  2:1199027/149‑1 (MQ=255)                                                                                  gCCCAGATCGATGTCGCCGTGGGTGACGGTGTAGATCTCGGCGGGGGCGGTCATGCGCACCAGGTGTTCATGCAGTTTTTCTGCTTCACCTGCGGCGTCGGTTTTGGTCAGCAGAATTCGGTCGGCGGCGCCAACCCGTCACTGGTCgt                                                              >  2:889258/1‑148 (MQ=255)                                                                                    ccAGATCGATGTCGCCGTGTGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGGAGTTTTTCTGCTTCGCCTGCGACGTCGGTGTTGGTGATCAGGATGAGGACGGCGTGGACACCCGGTGGATGGGCGATg                                                            >  2:823710/1‑149 (MQ=255)                                                                                    ccAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCAGAATGCGGTCGGCGTAGCCAACCTGTGACTGGGCGATg                                                            >  1:1231458/1‑149 (MQ=255)                                                                                    ccAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGACATCGGTTTTGGTCAGCAGAATGCGGTCGGCGTAGCCAACCTGTGACTGGGCGATg                                                            >  1:988182/1‑149 (MQ=255)                                                                                    ccAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTACGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCAGAATGCGGTCGGCGTAGCCAACCTGTGACTGGACGATg                                                            >  1:330187/1‑149 (MQ=255)                                                                                     cAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGCGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCAGAATGCGGTCGGCGTAGCCAACCTGTGACTGGGCGATgg                                                           <  1:386241/149‑1 (MQ=255)                                                                                      aGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCAGAATGCGGTCGGCGTAGCCAACCTGTGACTGGGCGATGGt                                                          <  1:534275/149‑1 (MQ=255)                                                                                         ttgATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCAGAATGCGGTCGGCGTAGCCAACCTGTGACTGGGCGATGGTGaa                                                       <  2:952434/147‑1 (MQ=255)                                                                                         tCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCAGAATGCGGTCGGCGTAGCCAACCTGTGACTGGGCGATGGTGaa                                                       <  1:903530/149‑1 (MQ=255)                                                                                         gcGATGGCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGGCGTCGGTTTTGGTCAGCAGAATGCGGTCGGCGTAGCCAATCTGTGACTGGGCGATGGTGaa                                                       <  1:623615/148‑1 (MQ=255)                                                                                          cGATGTCGCCGTGGTGGCCGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGCGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGGCAGCAGAATGCGGTCGGCGTAGCCAACCTGTGACTGGGCGATGGTGAAc                                                      <  2:1187325/149‑1 (MQ=255)                                                                                          cGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCAGAATGCGGTCGGCGTAGCCAACCTGTGACTGGGCGATGGTGAAc                                                      <  2:1005139/149‑1 (MQ=255)                                                                                             tGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCACGCCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCAGAATGCGGTCGGCGTAGCCAACCTGTGACTGGGCGATGGTGAACTgg                                                   >  1:985181/1‑149 (MQ=255)                                                                                             ggTCGCCGTGGGTTACGGTTTAGTCCTGTGCGCGGTCGTTTATGCGCGCCTGTTGGTCAGGCAGTTTTTCTGCGTCGCCTGCGACGTCGGTTTTGGTCAGCAGAATGCGGTCGGCGTAGCCAACCTGTGACTGGGCGATGGTGAACTgg                                                   <  1:570767/148‑1 (MQ=255)                                                                                               tCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCAGAAGGCGGTCGGCGTAGCCAACCTGTGACTGGGCGATGGTGAACTGGtt                                                 <  1:220254/149‑1 (MQ=255)                                                                                                  ccGTGGGGGAGGGTGGAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTTATGCAGGTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCAGAATGCGGTCGGCGTAGCCAACCTGTGACTGGGCGATGGTGAACTGGTTCAt                                              <  1:13756/149‑1 (MQ=255)                                                                                                   cGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCAGAATGCGGTCGGCGTAGCCAACCTGTGACTGGGCGATGGTGAACTGGTTCATc                                             <  1:637734/149‑1 (MQ=255)                                                                                                   cGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGAGCGCCAGGTGTTCATGCTGTTTTTCCGCGTCGCCTGCGACGTCGGCTTCGCCCATCAGATTGCGGTCGGCGTAGCCAACCTGTGACTGTGCGATGGTGAACTGGTTCAtg                                             <  2:942383/149‑2 (MQ=255)                                                                                                      gggTGACCGGGTAGACCGGTGCGCGCGCGTTTATGCGGGCTAGGCCTTCATGCACTTTTTCTGCTTCGCCTGCGACATCTGTTTTGGTCAGCAGAATGCGGTCGGCGTATCCAACCTGTGACTGGGCGATGGTGAACTGGTTCATCTGc                                          <  2:289029/149‑1 (MQ=255)                                                                                                          gACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCAGAATGCGGTCGGCGTAGCCAACCTGTGACTGGGCGATGGTGAACTGGTTCATCTGCTCAt                                      >  2:312037/1‑149 (MQ=255)                                                                                                              gtgtAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCAGAATGCGGTCGGCGTAGCCAACCTGTGACTGGGCGATGGTGAACTGGTTCATCTGCTCATCGGc                                  <  2:1231458/149‑1 (MQ=255)                                                                                                                  aGACCGGTGCGCGGGCGTTGATGCGCGCCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCAGAATGCGGTCGGCGTAGCCAACCTGTGACTGGGCGATGGTGAACTGGTTCATCTGCTCATCGGCATGt                              >  1:400326/1‑149 (MQ=255)                                                                                                                      cGGTGCGCGGGCGTTGATGCGCGCCAGGCGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCAGAATGCGGTCGGCGTAGCCAACCTGTGACTGGGCGATGGTGAACTGGTTCATCTGCTCATCGGCATGTACCg                          >  1:87366/1‑149 (MQ=255)                                                                                                                          gcgtGGGCGTTGATGCGCGCCAGTTGTTTATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCAGAATGCGGTCGGCGTAGCCAACCTGTGACTGGGCGATGGTGAACTGGTTCATCTGCTCATCGGCATGTACCGCATc                      <  1:1185886/149‑1 (MQ=255)                                                                                                                                     gATGCGCGCCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCAGAATGCGGTCGGCGTAGCCAACCTGTGACTGGGCGATGGTGAACTGGTTCATCTGCTCATCGGCATGTACCGCATCCACCAGCGCaa           <  1:257788/149‑1 (MQ=255)                                                                                                                                       tGCGCGCCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCAGAATGCGGTCGGCGTAGCCAACCTGTGACTGGGCGATGGTGAACTGGTTCATCTGCTCATCGGCATGTACCGCATCCACCAGCGCAAtg         >  1:260892/1‑148 (MQ=255)                                                                                                                                        gcgcgcCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCAGAATGCGGTCGGCGTAGCCAACCTGTGACTGGGCGATGGTGAACTGGTTCATCTGCTCATCGGCATGTACCGCATCCACCAGCGCAATCa        <  1:743477/149‑1 (MQ=255)                                                                                                                                         cgcgcCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGACATCGGTTTTGGTCAGCAGAATGCGGTCGGCGTAGCCAACCTGTGACTGGGCGATGGTGAACTGGTTCATCTGCTCATCGGCATGTACCGCATCCACCAGCGCAATCAc       <  2:988182/149‑1 (MQ=255)                                                                                                                                           cgcAAGGGGTTCATGCAGTTTGTCTGCTTCGCCTGCGACGTCGGTTTCGGTCAGCAGAATGCGGTCGGCGTAGCCAACCTGTGACTGGGCGATGGTGAACTGGTTCATCTGCTCATCGGCATGTACCGCATCCACCAGCGCAATCACGc     <  2:1241616/149‑1 (MQ=255)                                                                                                                                            gcCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCAGAATGCGGTCGGCGTAGCCAACCTGTGACTGGGCGATGGTGAACTGGTTCATCTGCTCATCGGCATGTACCGCATCCACCAGCGCAATCACGcc    <  1:690226/149‑1 (MQ=255)                                                                                                                                            gcCAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCAGAATGCGGTCGGCGTAGCCAACCTGTGACTGGGCGATGGTGAACTGGTTCATCTGCTCATCGGCATGTACCGCATCCACCAGCGCAATCACGcc    <  1:1088139/149‑1 (MQ=255)                                                                                                                                              cAGGTGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCAGAATGCGGTCGGCGTAGCCAACCTGTGACTGGGCGATGGTGAACTGGTTCATCTGCTCATCGGCATGTACCGCATCCACCAGCGCAATCACGCCGt  <  1:691180/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                 TTTATCTGCGATAAAGTGGAAACGCGGTTTGGTGCTGACGACGTTTTCTTCCAGCATAAAACCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGACGGTGTAGACCGGTGCGCGGGCGTTGATGCGCGCCAGGCGTTCATGCAGTTTTTCTGCTTCGCCTGCGACGTCGGTTTTGGTCAGCAGAATGCGGTCGGCGTAGCCAACCTGTGACTGGGCGATGGTGAACTGGTTCATCTGCTCATCGGCATGTACCGCATCCACCAGCGCAATCACGCCGT  >  NC_000913_3_hsa_tpiA/4588188‑4588477

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