Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913_3_hsa_tpiA	4,111,662	G→A	14.6%	L179L (CTG→TTG)	yiiQ ←	DUF1454 family putative periplasmic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_hsa_tpiA	4,111,662	0	G	A	14.6%	164.9 / 17.6	68	L179L (CTG→TTG)	yiiQ	DUF1454 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (27/31);  new base A (6/4);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 5.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

TAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCG  >  NC_000913_3_hsa_tpiA/4111516‑4111811                                                                                                                                                    |                                                                                                                                                      tAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGccc                                                                                                                                                    <  1:92560/150‑1 (MQ=255) tAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGccc                                                                                                                                                    >  4:62904/1‑150 (MQ=255)      gAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCtttt                                                                                                                                               <  1:509725/150‑1 (MQ=255)          aTGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGCTCAGCCCCTTTTCGcc                                                                                                                                           <  3:340315/150‑1 (MQ=255)           tGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCg                                                                                                                                          <  2:367515/150‑1 (MQ=255)                   ccGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCcgc                                                                                                                                  <  1:561916/150‑1 (MQ=255)                   ccGCAGGGCAGCAGGTCTGGGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACC‑TTCAAGCGAGTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCcgc                                                                                                                                  >  4:103124/1‑150 (MQ=255)                      ccgggCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGac                                                                                                                               >  2:13057/3‑150 (MQ=255)                         ggCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGacaac                                                                                                                            >  2:570131/1‑150 (MQ=255)                         ggCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGacaac                                                                                                                            >  4:313292/1‑150 (MQ=255)                         ggCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGacaac                                                                                                                            <  1:477486/150‑1 (MQ=255)                             gcagGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATaa                                                                                                                        <  3:190115/150‑1 (MQ=255)                              cagGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAAc                                                                                                                       <  3:228238/150‑1 (MQ=255)                                ggTCTGTGAAACAGTATAGAGGTTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGc                                                                                                                     <  2:334455/150‑1 (MQ=255)                                ggTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGc                                                                                                                     <  3:17496/150‑1 (MQ=255)                                 gTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCa                                                                                                                    <  3:253644/150‑1 (MQ=255)                                 gTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCa                                                                                                                    >  3:296286/1‑150 (MQ=255)                                 gTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCa                                                                                                                    >  1:221396/1‑150 (MQ=255)                                  tCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAg                                                                                                                   <  4:259011/150‑1 (MQ=255)                                     gtgAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGc                                                                                                                >  2:484347/1‑150 (MQ=255)                                      tgAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCa                                                                                                               >  4:299832/1‑150 (MQ=255)                                      tgAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCa                                                                                                               >  1:565405/1‑150 (MQ=255)                                        aaaCAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCAcc                                                                                                             >  1:124217/1‑150 (MQ=255)                                                taGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGtct                                                                                                     >  3:3913/1‑150 (MQ=255)                                                 agagATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGtctc                                                                                                    >  3:273386/1‑150 (MQ=255)                                                     aTTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGgtg                                                                                                <  2:50512/150‑1 (MQ=255)                                                     aTTCATCGGCACAAAGGCTTTGCTTTTTATAATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGgtg                                                                                                >  1:279754/1‑150 (MQ=255)                                                          tCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGtaata                                                                                           >  4:336427/1‑150 (MQ=255)                                                            ggCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAAc                                                                                         >  1:51517/1‑150 (MQ=255)                                                                  aaGGCTTTGCTTTTTGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgtttgt                                                                                   <  1:476529/150‑1 (MQ=255)                                                                       tttgctttTTGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTc                                                                              >  1:251321/1‑150 (MQ=255)                                                                       tttgctttTTGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTc                                                                              <  3:58637/150‑1 (MQ=255)                                                                         tgctttTTGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTccc                                                                            <  1:12527/150‑1 (MQ=255)                                                                            tttttGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCccgc                                                                         <  1:63974/150‑1 (MQ=255)                                                                              tttGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCccgccg                                                                       <  1:428226/150‑1 (MQ=255)                                                                              tttGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCccgccg                                                                       >  1:359528/1‑150 (MQ=255)                                                                                 gttaTTTATGCAAACCGTTTAGGCG‑TTCAAATAGCGCCGGTTTAATCGGTTCAACAGCGCAGGTCAGCCCCTTTTCGCCGTTTTCCGTGACAACTTAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTta                                                                    <  4:226667/147‑1 (MQ=255)                                                                                       tATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCg                                                              >  4:257568/1‑150 (MQ=255)                                                                                          tCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACt                                                           >  4:281985/1‑150 (MQ=255)                                                                                            aaaCC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGc                                                         >  3:373389/1‑150 (MQ=255)                                                                                            aaaCC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGc                                                         >  1:195340/1‑150 (MQ=255)                                                                                             aaCC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCa                                                        <  1:503543/150‑1 (MQ=255)                                                                                              aCC‑TTCAAGAGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAg                                                       >  1:395776/1‑150 (MQ=255)                                                                                                   tCAAGCCATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtttt                                                   >  4:292824/1‑150 (MQ=255)                                                                                                     gaGCGATTCAGATAGCGCCAGCTCAGTCGGTACACCTGCCATGGTCAGCCCTTTTTCGCCCTTGACCGCGACAACATAGCGCCGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtttttt                                                 <  4:145113/149‑1 (MQ=255)                                                                                                     aggCGATTCAGATAGCGCCAGATTAGTCGGTTCAACTGCGTTGGTCAGCTCTTTTTAGCCCTTGTCCGCGATAACATAACGCAGTCCACCTTCTTTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtttttt                                                 <  4:77754/148‑1 (MQ=255)                                                                                                     aaGCGATTCAGATAGGGCCAGCTAAATCGGTTCAACAGCGAAGGGCAGCCATTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGCCTGCAGttattt                                                 >  4:306686/1‑150 (MQ=255)                                                                                                      aGCGATTCAGATAGCGCCAGCATAATCGCATCAACAGCGATGGTCAGCCCCTTTTCGCCCTTGCACGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAACGACTGCAGTTTTTTc                                                <  2:292964/150‑1 (MQ=255)                                                                                                         gATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCTGCCTCTTTTCTCCGTTGTCCGCGACAACATAACGCAGTCCACCTTCTTTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTgg                                             <  2:344087/150‑1 (MQ=255)                                                                                                         gATTCAGAAGGCGCCAGCTTCATCGGTTCAACGGCGAAGGTCTGCGCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTgg                                             <  2:313582/150‑1 (MQ=255)                                                                                                             cgggtAGCGCTAGCTTAATCGGTTCAACAGCGAAGGTCAGCCTCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCCCCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCttt                                         <  4:33321/146‑1 (MQ=255)                                                                                                             cAGATAGCGCTAGCTTAATCGTTTCAACAACGGAGGTCGGCCTCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCttt                                         <  4:87081/150‑1 (MQ=255)                                                                                                                    cgcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGtttt                                  >  2:218649/1‑150 (MQ=255)                                                                                                                     ggcAGCTTAGTCGGTAGAACAGCGAGCGTAAGGCCCTTGTTTCCGTTGTCTGCTACAACATAACGCGGTGTTCCTTCTGTCCCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTg                                 <  2:179801/148‑1 (MQ=255)                                                                                                                     gcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTg                                 <  3:270008/150‑1 (MQ=255)                                                                                                                           ttAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAtt                           >  3:81516/1‑150 (MQ=255)                                                                                                                            tAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACCACATAACCCAGTGCACCTCCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGGTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTa                          <  3:103124/150‑1 (MQ=255)                                                                                                                            tAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTa                          <  3:249882/150‑1 (MQ=255)                                                                                                                             aaTCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTaa                         <  4:3913/150‑1 (MQ=255)                                                                                                                              aTCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAAt                        >  2:22402/1‑150 (MQ=255)                                                                                                                                 ggTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATggg                     <  3:299832/150‑1 (MQ=255)                                                                                                                                  gTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTGCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTTCCCGCCGTTAGTAGCGACTGCAGTTTTTTTTGGCTTTGTGTTGTGGTGATGAATgggg                    >  2:343874/1‑150 (MQ=255)                                                                                                                                  gTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATgggg                    >  3:10832/1‑150 (MQ=255)                                                                                                                                   ttCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                   <  1:17050/150‑1 (MQ=255)                                                                                                                                   ttCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                   >  3:372577/1‑150 (MQ=255)                                                                                                                                     cAACAGCGAGTGTCAGCCCCTTTTCGCCCTTGTCCGCGACAACATATCGCTGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGa                 <  4:41220/150‑1 (MQ=255)                                                                                                                                           cGAAGGTCAACCCCTTTTCGCCGTTTTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGATTTGTGTTTTGGTCATTAATGGGGTGAGTGTGc           >  3:139576/1‑150 (MQ=255)                                                                                                                                            gAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCg          >  1:512843/1‑150 (MQ=255)                                                                                                                                                    gCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCg  <  3:336427/150‑1 (MQ=255)                                                                                                                                                    gCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCg  <  1:159231/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      TAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACC‑TTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCG  >  NC_000913_3_hsa_tpiA/4111516‑4111811

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