Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_hsa_tpiA	1,544,022	1 bp→AA	coding (39/1389 nt)	narU ←	nitrate/nitrite transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_hsa_tpiA	1,544,021	1	.	A	86.2%	170.3 / 22.4	58	coding (40/1389 nt)	narU	nitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base . (7/1);  new base A (29/21);  total (36/22)
Fisher's exact test for biased strand distribution p-value = 1.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01
*	NC_000913_3_hsa_tpiA	1,544,022	0	C	A	98.3%	190.6 / NA	58	L13F (TTG→TTT)	narU	nitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (35/22);  total (36/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCG‑C‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAA  >  NC_000913_3_hsa_tpiA/1543879‑1544170                                                                                                                                                  |||                                                                                                                                                     ttGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCG‑A‑aaa‑aa                                                                                                                                                 >  2:319879/1‑148 (MQ=255)  tGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGTG‑A‑aaa‑aag                                                                                                                                                >  1:426778/1‑147 (MQ=255)   gagaTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCG‑A‑aaa‑aaga                                                                                                                                               <  1:472961/149‑4 (MQ=255)            ggTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTa                                                                                                                                       >  1:880865/1‑149 (MQ=255)                     gcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCa                                                                                                                              <  1:935865/149‑1 (MQ=255)                       tAAATAGCATCCAGACACAGAAAGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTGCCATTTCCAA‑A‑AAA‑AAATAACGACTATTTTTCTCCtt                                                                                                                            >  2:865130/1‑149 (MQ=255)                         aaTAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCGGGTTTCCA‑ATCGCGAA‑AAA‑AGAGAACGGCTATTTTTCTCAtttt                                                                                                                          >  1:397482/1‑149 (MQ=255)                              cATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGt                                                                                                                     <  2:324320/149‑1 (MQ=255)                                    gACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCata                                                                                                               >  2:549538/1‑149 (MQ=255)                                     acacAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTATTCGAGAAATATTTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCG‑AAAAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCatat                                                                                                              >  1:229590/1‑149 (MQ=255)                                       acaGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTg                                                                                                            >  1:350454/1‑149 (MQ=255)                                              gCAAGAAGTATACAACTTACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTCCTTCAAA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTca                                                                                                     <  1:1101473/149‑1 (MQ=255)                                              gCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAC‑AAA‑AGATAACGACTATTTTTTTTATTTTGCAGTGCCATATTATTCCTca                                                                                                     >  2:662942/1‑149 (MQ=255)                                                   aagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGc                                                                                                >  1:710080/1‑149 (MQ=255)                                                     gTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTc                                                                                                      >  2:556281/1‑141 (MQ=255)                                                     gTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTc                                                                                                      <  1:556281/141‑1 (MQ=255)                                                        gACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                                         >  1:232362/1‑135 (MQ=255)                                                        gACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                                         <  2:232362/135‑1 (MQ=255)                                                        gACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcacaca                                                                                           >  1:442806/1‑149 (MQ=255)                                                        gACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcacaca                                                                                           >  1:478218/1‑149 (MQ=255)                                                         aCAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTTCAGTGCCATATTGTTCCTCACATGCACACAt                                                                                          >  1:270063/1‑149 (MQ=255)                                                         aCAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAt                                                                                          >  1:350811/1‑149 (MQ=255)                                                           aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg                                                                                        >  2:1047552/1‑149 (MQ=255)                                                           aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg                                                                                        <  1:1047552/149‑1 (MQ=255)                                                           aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg                                                                                        >  1:75679/1‑149 (MQ=255)                                                           aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg                                                                                        <  1:533225/149‑1 (MQ=255)                                                              tgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCGATATTGTTCCTCACATGCACACATTGGTa                                                                                     >  1:1006247/1‑149 (MQ=255)                                                              tgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTa                                                                                     <  1:549538/149‑1 (MQ=255)                                                               gactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTaa                                                                                    >  1:911795/1‑149 (MQ=255)                                                                     tatCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaa                                                                              >  2:802106/1‑149 (MQ=255)                                                                       tCCAGAGATTTCTTCGAGCAATATGATCTGTAAAGTAGTGCCAGACGGCCGGACTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaa                                                                            >  1:871392/1‑149 (MQ=255)                                                                        ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa                                                                           <  1:763468/149‑1 (MQ=255)                                                                        ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa                                                                           <  1:319879/149‑1 (MQ=255)                                                                         cAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAg                                                                          >  1:220476/1‑149 (MQ=255)                                                                          agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCG‑AAAAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                                         >  1:420643/1‑149 (MQ=255)                                                                               tttCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaaaa                                                                    >  1:713238/1‑149 (MQ=255)                                                                                   ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACg                                                                >  2:758671/1‑149 (MQ=255)                                                                                          aaTATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCG‑A‑AAACAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGAACAAACACGAGGTAAg                                                         >  2:1098810/1‑149 (MQ=255)                                                                                           atatGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAGAACACGAGGTAAgg                                                        >  2:1148337/1‑149 (MQ=255)                                                                                               gTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCa                                                    >  1:457486/1‑149 (MQ=255)                                                                                                     cTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCata                                                                                                               >  2:1073773/1‑84 (MQ=255)                                                                                                     cTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCata                                                                                                               <  1:1073773/84‑1 (MQ=255)                                                                                                      tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                                             >  1:208933/1‑149 (MQ=255)                                                                                                      tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                                             >  2:893072/1‑149 (MQ=255)                                                                                                             tCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGttatta                                      <  2:880865/149‑1 (MQ=255)                                                                                                              cccAGAATGCCGGATTTTCTGGTTTCCA‑TTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAg                                     <  1:84325/149‑1 (MQ=255)                                                                                                               ccAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGca                                                                                               <  2:671606/90‑1 (MQ=255)                                                                                                               ccAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGca                                                                                               >  1:671606/1‑90 (MQ=255)                                                                                                                        cGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGt                                           <  2:1089985/133‑1 (MQ=255)                                                                                                                        cGGATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGt                                           >  1:1089985/1‑133 (MQ=255)                                                                                                                          gATTTTCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGa                         <  2:911795/149‑1 (MQ=255)                                                                                                                               tCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                                             >  2:363437/1‑124 (MQ=255)                                                                                                                               tCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                                             <  1:363437/124‑1 (MQ=255)                                                                                                                               tCTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTgg                    <  1:313787/149‑1 (MQ=255)                                                                                                                                cTGGTTTCCA‑GTCGCGAA‑AAA‑AGATAACGACTATTTTCCTCATTTTTCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTggg                   <  2:654173/149‑1 (MQ=255)                                                                                                                                     ttCCA‑GTCGCGAA‑AAA‑AGATAACTACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAACGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAg              <  1:986189/149‑1 (MQ=255)                                                                                                                                       ccA‑GTCGCGAA‑AAA‑AGATAACGACTATTTCTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTg            <  2:272168/149‑1 (MQ=255)                                                                                                                                        cA‑GTCGCG‑AAAAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGc           >  2:207466/1‑149 (MQ=255)                                                                                                                                              cgcgaaaaa‑aGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGAt       <  2:710080/143‑1 (MQ=255)                                                                                                                                                   aaaa‑aGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa  <  2:75679/148‑1 (MQ=255)                                                                                                                                                  |||                                                                                                                                                     TTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTCGCG‑C‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAA  >  NC_000913_3_hsa_tpiA/1543879‑1544170

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