Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913_3_hsa_tpiA	3,905,355	G→A	7.4%	D81D (GAC→GAT)	bglF ←	fused beta‑glucoside‑specific PTS enzymes: IIA component/IIB component/IIC component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_hsa_tpiA	3,905,355	0	G	A	7.4%	200.4 / 4.5	67	D81D (GAC→GAT)	bglF	fused beta‑glucoside‑specific PTS enzymes: IIA component/IIB component/IIC component
Reads supporting (aligned to +/- strand):  ref base G (33/29);  new base A (2/3);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

AGTGAGCGCCAGAGCCAGCATACCTTTCAAGATCCCGGTTGCCGCCATCAAACCGATCAGAGGCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAA  >  NC_000913_3_hsa_tpiA/3905208‑3905504                                                                                                                                                    |                                                                                                                                                      aGTGAGCGCCAGAGCCAGCATACCGTTCAAGATCCCGGTTGCCGCCATCAAACCGATCAGAGGCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTc                                                                                                                                                     <  2:477606/150‑1 (MQ=255)     aGCGCCAGAGCCAGCATACCTTTCAAGATCCCGGTTGCCGCCATCAAACCGATCAGAGGCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAgg                                                                                                                                                 >  2:149226/1‑150 (MQ=255)               ccAGCATACCTTTCAAGATCCCGGTTGCCGCCATCAAACCGATCAGAGGCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCacac                                                                                                                                       <  4:282994/150‑1 (MQ=255)               ccAGCATACCTTTCAAGACCCCGGTTGCCGCCATCAAACCGATCAGAGGCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCacac                                                                                                                                       >  4:34892/1‑150 (MQ=255)                  gCATACCTTTCAAGATCCCGGTTGCCGCCATCAAACCGATCAGAGGCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGt                                                                                                                                    <  3:357120/150‑1 (MQ=255)                  gCATACCTTTCAAGATCCCGGTTGCCGCCATCAAACCGATCAGAGGCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGt                                                                                                                                    <  4:382446/150‑1 (MQ=255)                      aCCTTTCAAGATCCCGGTTGCCGCCATCAAACCGATCAGAGGCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAAc                                                                                                                                >  4:12130/1‑150 (MQ=255)                        cTTTCAAGATCCCGGTTGCCGCCATCAAACCGATCAGAGGCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCg                                                                                                                              >  1:291227/1‑150 (MQ=255)                        cTTTCAAGATCCCGGTTGCCGCCATCAAACCGATCAGAGGCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCg                                                                                                                              >  3:160791/1‑150 (MQ=255)                        cTTTCAAGATCCCGGTTGCCGCCATCAAACCGATCAGAGGCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCg                                                                                                                              >  1:205771/1‑150 (MQ=255)                        cTTTCAAGATCCCGGTTGCCGCCATCAAACCGATCAGAGGCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCg                                                                                                                              >  2:400236/1‑150 (MQ=255)                          ttCAAGATCCCGGTTGCCGCCATCAAACCGATCAGAGGCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGcc                                                                                                                            <  4:250327/150‑1 (MQ=255)                           tCAAGATCCCGGTTGCCGCCATCAAACCGATCAGAGGCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCa                                                                                                                           <  1:543995/150‑1 (MQ=255)                            cAAGATCCCGGTTGCCGCCATCAAACCGATCAGAGGCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAg                                                                                                                          >  3:250898/1‑150 (MQ=255)                                aTCCCGGTTGCCGCCATCAAACCGATCAGAGGCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAg                                                                                                                      >  1:39796/1‑150 (MQ=255)                                 tCCCGGTTGCCGCCATCAAACCGATCAGAGGCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGa                                                                                                                     <  1:400501/150‑1 (MQ=255)                                    cGGTTGCCGCCATCAAACCGATCAGAGGCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACAt                                                                                                                  >  4:243621/1‑150 (MQ=255)                                      gTTGCCGCCATCAAACCGATCAGAGGCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCg                                                                                                                <  2:455896/150‑1 (MQ=255)                                      gTTGCCGCCATCAAACCGATCAGAGGCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCg                                                                                                                >  3:156094/1‑150 (MQ=255)                                           cgccATCAAACCGATCAGAGGCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCcac                                                                                                           >  3:146996/1‑150 (MQ=255)                                            gccATCAAACCGATCAGAGGCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCcaca                                                                                                          <  1:30449/150‑1 (MQ=255)                                             ccATCAAACCGATCAGAGGCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACAt                                                                                                         <  4:346640/150‑1 (MQ=255)                                             ccATCAAACCGATCAGAGGCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACAt                                                                                                         <  2:490664/150‑1 (MQ=255)                                              cATCAAACCGATCAGAGGCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATg                                                                                                        >  3:365105/1‑150 (MQ=255)                                                   aaCCGATCAGAGGCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCATCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTAc                                                                                                   >  3:341314/1‑150 (MQ=255)                                                   aaCCGATCAGAGGCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCATCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTAc                                                                                                   >  3:183621/1‑150 (MQ=255)                                                   aaCCGATCAGAGGCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTCTAAAGCATGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTAc                                                                                                   <  3:379254/150‑1 (MQ=255)                                                    aCCGATCAGAGGCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTAcc                                                                                                  <  3:259349/150‑1 (MQ=255)                                                        aTCAGAGGCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATg                                                                                              >  1:385971/1‑150 (MQ=255)                                                         tCAGAGGCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGa                                                                                             <  4:199443/150‑1 (MQ=255)                                                           agagGCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGacc                                                                                           >  1:11802/1‑150 (MQ=255)                                                             agGCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGaccac                                                                                         <  3:232511/150‑1 (MQ=255)                                                               gCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCt                                                                                       <  4:336471/150‑1 (MQ=255)                                                                   aaaaaTACCTGAAATAACATAAACAAAGCTGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGaaa                                                                                   <  3:308918/150‑1 (MQ=255)                                                                    aaaaTACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAAc                                                                                  >  2:431564/1‑150 (MQ=255)                                                                    aaaaTACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAAc                                                                                  >  4:146317/1‑150 (MQ=255)                                                                        tACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGc                                                                              <  1:276326/150‑1 (MQ=255)                                                                          ccTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCa                                                                            >  1:562057/1‑150 (MQ=255)                                                                                      tAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTccacc                                                                >  1:187633/1‑150 (MQ=255)                                                                                        aaCAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTCGCGCTTTTTCGTCAAGGCCTGCCACACTGTCAACCGCCAGGAAGACATCGGCCACATGGTTACCTCTGACCACCTGAACCTGGCCACCGCTTTccaccac                                                              <  4:159875/150‑1 (MQ=255)                                                                                          cacaggggTTTAGCAGATTACCTTTATCATCAGTGTCCGGTGCCTGTTGCGCTTTTTTGTCAAGGCCTGCCACACTTTTAACCGCCAGTATGTCATCGGCCTCCTGGTTAACTCTGCCCATCTGAAACTGTCCACAGCGTTccaccacca                                                            <  2:572491/144‑1 (MQ=255)                                                                                          aaaggCGTTTTAGCAGAATACCTATAGCATCATTTATGGGAGTCTGTTGCGCATTTTCGTCAAGGCCTGCCACACTGTTAAGCGCCAGGAAGGCATCGGCCACATGGTTACATATGACCACCTGAAACTGGCCACCGCTTTccaccacca                                                            <  2:295525/146‑1 (MQ=255)                                                                                            aaGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCata                                                          <  2:371947/150‑1 (MQ=255)                                                                                              gCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCataat                                                        <  2:584534/150‑1 (MQ=255)                                                                                              gCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCataat                                                        >  2:461880/1‑150 (MQ=255)                                                                                                  tttAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCataataata                                                    >  4:215386/1‑150 (MQ=255)                                                                                                    tAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGTCCTTCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATAcc                                                  >  4:333246/1‑150 (MQ=255)                                                                                                       cAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCggg                                               <  3:243621/150‑1 (MQ=255)                                                                                                        aGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCgggg                                              >  2:503371/1‑150 (MQ=255)                                                                                                        aGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCgggg                                              <  4:365105/150‑1 (MQ=255)                                                                                                           ttACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTc                                           >  3:298522/1‑150 (MQ=255)                                                                                                               cTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAATTCGCCCACACAGTTAACCGCCAGGAAGACATCGGCCACATGATTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCtttt                                       >  2:26023/1‑150 (MQ=255)                                                                                                                    tcatcaTTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGt                                  <  3:107836/150‑1 (MQ=255)                                                                                                                      atcatTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCAGCATAATAATACCGGGGGTCTTTTTCAGTAc                                >  1:376640/1‑150 (MQ=255)                                                                                                                      atcatTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCttttttattat                                >  4:294339/1‑144 (MQ=255)                                                                                                                      atcatTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTCAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCAGGGGTCTTTTTCAGTAc                                >  2:375223/1‑150 (MQ=255)                                                                                                                       tcatTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTAcc                               >  2:135495/1‑150 (MQ=255)                                                                                                                          ttttCCGGTGCCTGTTGCGCTTTTTCATCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACctct                            <  4:341314/150‑1 (MQ=255)                                                                                                                          ttttCCGGTGCCTGTTGCGCTTTTTCATCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACctct                            <  4:183621/150‑1 (MQ=255)                                                                                                                                ggTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACCTCTGCTTgc                      >  2:461471/1‑150 (MQ=255)                                                                                                                                     cTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACCTCTGCTTGCgcttt                 <  2:504588/150‑1 (MQ=255)                                                                                                                                         tGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATATTACCGGGGGTCTTTTTCAGTACCTCTGCTTGCgctttgctt             <  2:194216/150‑1 (MQ=255)                                                                                                                                           cgcTTTTTCGTTAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACCTCTGCTTGCGCTTTGCTTTc           >  2:80834/1‑150 (MQ=255)                                                                                                                                             cTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCAt         <  1:579569/150‑1 (MQ=255)                                                                                                                                               ttttCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATcc       <  3:296776/150‑1 (MQ=255)                                                                                                                                               ttttCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATcc       >  4:136883/1‑150 (MQ=255)                                                                                                                                               ttttCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATcc       >  1:199341/1‑150 (MQ=255)                                                                                                                                                    gTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCtttaa  >  2:43359/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                      AGTGAGCGCCAGAGCCAGCATACCTTTCAAGATCCCGGTTGCCGCCATCAAACCGATCAGAGGCGTAAAAATACCTGAAATAACATAAACAAAGCGGTTTAGCAGATTACCTTTATCATCATTTTCCGGTGCCTGTTGCGCTTTTTCGTCAAGGCCTGCCACACTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAA  >  NC_000913_3_hsa_tpiA/3905208‑3905504

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: ‑