Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_hsa_tpiA	1,544,022	C→A	L13F (TTG→TTT)	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,022	0	C	A	96.8%	188.3 / NA	62	L13F (TTG→TTT)	narU	nitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (33/27);  total (34/28)
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

TGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑C‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAG  >  NC_000913_3_hsa_tpiA/1543880‑1544158                                                                                                                                                 |                                                                                                                                              tGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑A‑a‑aa‑aag                                                                                                                                      >  2:1089874/1‑147 (MQ=255)     gttAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑AAA‑AA‑A‑GATaa                                                                                                                                  <  2:701624/148‑1 (MQ=255)      ttAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAAc                                                                                                                                 >  2:851923/1‑149 (MQ=255)        aaCGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGa                                                                                                                               <  1:641447/149‑1 (MQ=255)         aCGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑A‑A‑AA‑AGGATAACGAc                                                                                                                              >  1:246619/1‑149 (MQ=255)            gTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTAt                                                                                                                           >  1:889920/1‑149 (MQ=255)               aCTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTGCCTTTATTTTCCCAAAAGTCCGGATTTTCTGGTTTCCAGTCGCCA‑C‑A‑AA‑A‑GCTAACGACTAtttt                                                                                                                        >  1:551220/1‑149 (MQ=255)                cTGCGCTAAATAGAATCCAGACACAGAAGGCAAGAAGTAGACAAATGACTGATATCCAGAGATTTCTTCGAGCAAGATGTTTTCCTTTATTTTCTCAGAAGGCCGGATTTTCTTGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTAttttt                                                                                                                       <  1:547324/149‑1 (MQ=255)                 tGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑AAA‑AA‑A‑GATAACGACTAtttttt                                                                                                                      >  2:90432/1‑148 (MQ=255)                 tGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCAAGAAGGCCGGATTTTCTGTTTTCCAGTCG‑CGAA‑A‑AA‑A‑TATAACGACTATTTTtc                                                                                                                      >  2:852486/1‑149 (MQ=255)                         aTAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTg                                                                                                              <  1:720697/149‑1 (MQ=255)                             cATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGt                                                                                                          <  1:470667/149‑1 (MQ=255)                                ccAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGcc                                                                                                       >  1:585965/1‑149 (MQ=255)                                    acacAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCatat                                                                                                   >  1:699609/1‑149 (MQ=255)                                       aaGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGt                                                                                                <  2:706229/148‑1 (MQ=255)                                          caGGCAATAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTcc                                                                                             <  2:1239860/148‑1 (MQ=255)                                             gCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑AAA‑AA‑A‑GATAAGGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTca                                                                                          >  1:351405/1‑149 (MQ=255)                                               aagaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcaca                                                                                        <  1:952673/149‑1 (MQ=255)                                               aagaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcaca                                                                                        >  1:130387/1‑149 (MQ=255)                                               aagaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGCTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcaca                                                                                        >  1:446951/1‑149 (MQ=255)                                                  aagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑AAA‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGc                                                                                     >  1:214481/1‑149 (MQ=255)                                                       gACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcacaca                                                                                >  2:1117308/1‑149 (MQ=255)                                                       gACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcacaca                                                                                >  1:1117744/1‑149 (MQ=255)                                                        aCAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCCGTGCCATATTGTTCCTCACATGCACACAt                                                                               >  1:212989/1‑149 (MQ=255)                                                          aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CC‑C‑AGAA‑A‑GATCAAGATGTTTGTTGTTATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg                                                                             <  1:294347/149‑1 (MQ=255)                                                          aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTAACAGAAGACAGGCTTTTCTGCTTTCCAGGCC‑AGAA‑A‑AA‑A‑GATAACGACAATTATTTCCATTTTGCAGCGCCAAATTTTTCCTCACATGCACACATTg                                                                             >  2:294347/1‑149 (MQ=255)                                                              gactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTaa                                                                         <  2:1117744/149‑1 (MQ=255)                                                                 gcatatCCCGGGATTTCCTCCAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGTATTTTCTGGTTTCCAGCCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGa                                                                      <  1:118455/147‑1 (MQ=255)                                                                   atatCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaa                                                                    >  2:694606/1‑149 (MQ=255)                                                                      tCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATTaaaaaa                                                                 >  2:122497/1‑149 (MQ=255)                                                                      tCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaa                                                                 >  1:706219/1‑149 (MQ=255)                                                                       ccAGATATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa                                                                <  1:261301/149‑1 (MQ=255)                                                                       ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa                                                                <  1:66971/149‑1 (MQ=255)                                                                         agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTAttt                                                                                                                         <  1:90241/90‑1 (MQ=255)                                                                         agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                              <  1:92899/149‑1 (MQ=255)                                                                         agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑AAA‑AA‑A‑GATAACGACTAttt                                                                                                                         >  2:90241/1‑90 (MQ=255)                                                                             aTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaaa                                                          <  1:737513/149‑1 (MQ=255)                                                                                  ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACg                                                     >  1:1015291/1‑149 (MQ=255)                                                                                       gCAATATGTTATCCTTTATTTACCCAAAAGACCGGCTTTTCTGTTTTCCATTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTGCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑a                                               >  2:1045574/1‑149 (MQ=255)                                                                                         aaTATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAg                                             >  2:602707/1‑149 (MQ=255)                                                                                              gTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑AAA‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCa                                        <  1:324909/149‑1 (MQ=255)                                                                                               ttttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑AAA‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATTCACACATGTGTAATGCAAACAAGACAAAACCCGAGGA‑AAGGCGCaa                                       >  2:1130845/1‑149 (MQ=255)                                                                                               ttttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑AAA‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCaa                                       >  2:488745/1‑149 (MQ=255)                                                                                               ttttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTGTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCaa                                       <  1:544442/149‑1 (MQ=255)                                                                                                   ccTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCAATAGc                                   <  1:177699/149‑1 (MQ=255)                                                                                                     tttATTTTCCCAGAAGGCCGTATTTTCTGGTTTCCAGTCG‑CG‑A‑A‑AA‑ATGATAACGACTATTTTTCTCCTTTTGCAGTGCCATATTGTTCCTCACATGCACTCATTGGTCATGAAAAAATGACAAGACACGAAGT‑ACAGCGCAACAgtcg                                 >  2:312262/1‑146 (MQ=255)                                                                                                     tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCAATAGccg                                 >  1:566497/1‑148 (MQ=255)                                                                                                     tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCAATAGCCa                                 >  2:1065993/1‑149 (MQ=255)                                                                                                         ttttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCAATAGCCAGtta                             >  1:1011879/1‑149 (MQ=255)                                                                                                           ttCCCAGGAGGCCGGATTTTCTGGTTTCCAGTCG‑AG‑AAA‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCAATAGCCAGttatt                           <  1:532564/149‑1 (MQ=255)                                                                                                              ccAGAAGGCCGGAGTTTCTGGTTACCAGAAG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCAATAGCCAGTTATTAGa                        <  2:319947/149‑1 (MQ=255)                                                                                                                 gAAGGCGGGTTTTTCTGGTTTTGAATCG‑TGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGGGGT‑AAGGCGCAATAGCCAGTTATTAGAAtt                     <  2:24504/149‑1 (MQ=255)                                                                                                                   aGGCCTGATTGTCTTGTTTCAAGGCT‑CG‑A‑A‑AAGA‑GATAACGACTATTTTTCTCAATTTGCAGTGCCATATTGTTCCTCATATGCACACAATGGTAAAGAAAAAAAGACAAAACACGAGGT‑AAGGCGCAATAGCCAGTTATTAGAATTaa                   <  2:601697/149‑1 (MQ=255)                                                                                                                      ccGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg                                                                             >  2:597364/1‑89 (MQ=255)                                                                                                                      ccGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg                                                                             <  1:597364/89‑1 (MQ=255)                                                                                                                         gATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCAATAGCCAGTTATTAGAATTAAGGATGa             >  1:817975/1‑149 (MQ=255)                                                                                                                             ttCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTg         <  2:343702/149‑1 (MQ=255)                                                                                                                             ttCTGGTTTCAAGTCT‑CG‑A‑ATAA‑A‑GATAACGGCTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTTACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCAATCGCCAGTTATTAGAATTGAGGATGACTTg         <  2:615720/149‑1 (MQ=255)                                                                                                                                  gTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAA‑AAGCGGTGATGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGa    >  1:1029902/1‑149 (MQ=255)                                                                                                                                   tttCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAACATGATTAAACACGAGAT‑AAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGaa   >  1:603090/1‑149 (MQ=255)                                                                                                                                   tttCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGaa   <  1:1141241/149‑1 (MQ=255)                                                                                                                                    ttCCAGTCG‑CG‑AAA‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGGAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAg  <  2:889920/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                              TGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑C‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAG  >  NC_000913_3_hsa_tpiA/1543880‑1544158

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