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	85.7%	145.0 / 18.2	49	coding (40/1389 nt)	narU	nitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base . (6/1);  new base A (24/18);  total (30/19)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	NC_000913_3_hsa_tpiA	1,544,022	0	C	A	100.0%	164.4 / NA	49	L13F (TTG→TTT)	narU	nitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (30/19);  total (30/19)

TTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAAG  >  NC_000913_3_hsa_tpiA/1543879‑1544171                                                                                                                                                 |||                                                                                                                                                      ttGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑aa‑aaa                                                                                                                                                  <  2:685173/149‑2 (MQ=255)    agaTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑aa‑aaagat                                                                                                                                               >  2:346039/1‑145 (MQ=255)           cGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACCACTGACTGATATCCTGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑AACAAGTTAGCGACt                                                                                                                                         >  1:964562/1‑149 (MQ=255)            ggTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTa                                                                                                                                        <  1:346039/149‑1 (MQ=255)                   gcgcTAAATAGCATCCAGACACAGAAGGCAAGCAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTAttttt                                                                                                                                   <  1:976880/147‑1 (MQ=255)                   gcgcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTAttttt                                                                                                                                   >  2:976880/1‑147 (MQ=255)                           aaaaaaccAGAAACAGAAGGCAAGAAGTAGACAACTTACTGATATCCAAAGATTTCTTCTAGCAATATGTTTCCCTTTAATTTCCCAGAAGGCCGGATTTTCTGGTTAAAAAAAAAAAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGc                                                                                                                         <  2:704165/143‑1 (MQ=255)                                tCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGGGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGc                                                                                                                    >  1:190096/1‑149 (MQ=255)                                   aGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGAAGTGCCat                                                                                                                 >  1:276343/1‑149 (MQ=255)                                              gCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGGCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTca                                                                                                      <  1:451753/149‑1 (MQ=255)                                                aagaagTAGACAACTGACTGATATCCAGAGTTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTAATTTTGCAGTGCCATATTGTTCCTaaca                                                                                                    >  2:904640/1‑149 (MQ=255)                                                aagaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcaca                                                                                                    >  1:278188/1‑149 (MQ=255)                                                  gaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTAGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATg                                                                                                  <  1:490783/149‑1 (MQ=255)                                                      tAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcaca                                                                                              <  2:725697/149‑1 (MQ=255)                                                      tAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcaca                                                                                              >  1:653824/1‑149 (MQ=255)                                                             ctgaaacatatACAGAAAATTCATCTAGCAAAATTTATTCCTTTATTTTCCCAGAAAGCCGGATTTTCTGTATAAAAAAAAAAAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGt                                                                                       <  2:772466/142‑1 (MQ=255)                                                                agtgaTATCCAGAGATTTCGTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAAt                                                                                    <  2:259074/147‑1 (MQ=255)                                                                actgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACCTGCACACATTGGTAAt                                                                                    >  2:995396/1‑149 (MQ=255)                                                                     tatCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaa                                                                               >  2:259693/1‑149 (MQ=255)                                                                       tCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaa                                                                             >  1:516378/1‑149 (MQ=255)                                                                        ccAGAGATTTCTTCTAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa                                                                            <  1:639196/149‑1 (MQ=255)                                                                        ccAGAGATTTCTTCGAGCCATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCCCACATGCACACATTGGTAATTaaaaaaa                                                                            >  2:101025/1‑149 (MQ=255)                                                                        ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGACGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑AAGAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAGTGGTAATGaaaaaaa                                                                            >  1:151346/1‑149 (MQ=255)                                                                         cAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGt                                                                                                                      <  1:517935/106‑1 (MQ=255)                                                                         cAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGt                                                                                                                      >  2:517935/1‑106 (MQ=255)                                                                         cAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATCCACACATTGGTATGGaaaaaaaa                                                                           >  2:622001/1‑148 (MQ=255)                                                                          agggaTTTCTTCGAGGAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                                          <  2:151346/146‑1 (MQ=255)                                                                          agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                                          >  2:1025611/1‑149 (MQ=255)                                                                                   ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACg                                                                 >  2:143966/1‑149 (MQ=255)                                                                                         cAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTaa                                                           >  2:39002/1‑149 (MQ=255)                                                                                                ttttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCaa                                                    <  1:143966/149‑1 (MQ=255)                                                                                                   tCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTTTTTTTTTTTTTTGACATATTGATCCTAACATGCACACCTTTGTAAGGAAAAAAAGACAAAATACGACGTTAGgcgtattat                                                 >  2:1028574/1‑143 (MQ=255)                                                                                                      ttttttGTCCCGGAAGGCCGGATTTTCTGGTTTTTAGTTGCGAA‑AA‑AAGATAACGACTATTTTACTCATTTTGCAGTGCCATATTTTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                                              <  2:925558/149‑1 (MQ=255)                                                                                                      tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCATTAGccg                                              >  1:905112/1‑148 (MQ=255)                                                                                                      tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                                              >  1:420501/1‑149 (MQ=255)                                                                                                       ttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAg                                             <  1:39002/149‑1 (MQ=255)                                                                                                       ttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑AATAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAg                                             >  1:584651/1‑149 (MQ=255)                                                                                                            ttCCCATAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGttatt                                        <  2:294525/149‑1 (MQ=255)                                                                                                               ccAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCACTAGCCAGTTATTAGa                                     >  2:750410/1‑149 (MQ=255)                                                                                                                cAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGaa                                    >  2:885659/1‑149 (MQ=255)                                                                                                                cAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGaa                                    >  1:568688/1‑149 (MQ=255)                                                                                                                 aGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTa                                                                                      <  2:938397/98‑1 (MQ=255)                                                                                                                 aGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTa                                                                                      >  1:938397/1‑98 (MQ=255)                                                                                                                       ccGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGAAGTGTAAAAATTTCCAAAAACGAGCTAAGGCGCAATAGCCAGTTATTAGAATTAAgaa                             >  2:337969/1‑147 (MQ=255)                                                                                                                        cGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGAt                            >  2:821723/1‑149 (MQ=255)                                                                                                                         ggATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATg                           >  2:312360/1‑149 (MQ=255)                                                                                                                               tCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGCACATCCACAAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTgg                     <  1:259693/149‑1 (MQ=255)                                                                                                                                cTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTggg                    <  2:568688/149‑1 (MQ=255)                                                                                                                                   gTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGa                 >  1:314390/1‑149 (MQ=255)                                                                                                                                            tcgcgaaaa‑aaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGa         <  1:867703/142‑1 (MQ=255)                                                                                                                                            tcgcgaaaa‑aaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGa         <  1:885659/142‑1 (MQ=255)                                                                                                                                            tcgcgaaaa‑aaGATAACGACTATTTTTCTCATTTTGCAGGGCCACATCGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGa         <  1:622001/142‑1 (MQ=255)                                                                                                                                              gcgaaaa‑aaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGAtt       >  1:442508/6‑149 (MQ=255)                                                                                                                                              gcgaaaa‑aaGATAACGACTATTTTTCTCATTTTGCAGGGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGAtt       <  1:314240/144‑1 (MQ=255)                                                                                                                                               cgaaaa‑aaGATAGCGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTTGTAATGAAAAAGAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTg      >  2:817639/5‑149 (MQ=255)                                                                                                                                                gaaaa‑aaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGa     >  2:543231/4‑149 (MQ=255)                                                                                                                                                 aaaa‑aaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaa    <  1:821723/147‑1 (MQ=255)                                                                                                                                                  aaa‑aaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTaa                                >  1:352714/2‑120 (MQ=255)                                                                                                                                                  aaa‑aaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTaa                                <  2:352714/119‑1 (MQ=255)                                                                                                                                                  aaa‑aaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa   <  2:420501/148‑1 (MQ=255)                                                                                                                                                  aaa‑aaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa   >  2:997269/2‑149 (MQ=255)                                                                                                                                                  aaa‑aaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCAAATGAACACTTTGGGAATGAAAAAAAGACAATACACGAGGTAAGGCGCACTAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGAGCTGATTTaaa   >  2:949205/2‑149 (MQ=255)                                                                                                                                                   aa‑aaGATAACGAATAGTTTTCTTATTTTGCAGTGCAATACTGTCAGTAACATTGACACAGTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAAg  <  2:940410/149‑1 (MQ=255)                                                                                                                                                 |||                                                                                                                                                      TTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAAG  >  NC_000913_3_hsa_tpiA/1543879‑1544171

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