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	98.3%	189.1 / NA	58	L13F (TTG→TTT)	narU	nitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (26/31);  minor base G (1/0);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTC‑GC‑G‑C‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCT  >  NC_000913_3_hsa_tpiA/1543884‑1544162                                                                                                                                               |                                                                                                                                                 aTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGCAGGCCGGATTTTCTGGTTTCCA‑GTC‑GC‑G‑A‑aa‑a‑aagctaa                                                                                                                                     >  1:533026/1‑143 (MQ=255) aTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATa‑a                                                                                                                                     >  1:957586/1‑149 (MQ=255) aTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGGAGACAACTGAGTGATATCCAGAGATTTCTTCGAGCCATATGTTTTCCTTTATTTTCCCGGAAGCCCGGTTTTTCTGGTTTCCA‑GTC‑GC‑G‑G‑AA‑ACAGATa‑a                                                                                                                                     >  1:162612/1‑149 (MQ=255)  ttAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGTTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑Ac                                                                                                                                    >  2:137409/1‑149 (MQ=255)  ttAACGGTAACTGCGCTAAATAGCATCCAGACAAAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑Ac                                                                                                                                    <  2:4151/149‑1 (MQ=255)    aaCGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGa                                                                                                                                  >  1:566319/1‑149 (MQ=255)        gTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTAt                                                                                                                              >  1:819591/1‑149 (MQ=255)             tGCGCTAAATAGCATCCAGACACCAAAGGCAAGAAGTAGACAACTGACTAATACCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCTGATTTTCTGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTtc                                                                                                                         <  1:355556/149‑1 (MQ=255)             tGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTtc                                                                                                                         <  2:819591/149‑1 (MQ=255)               cgcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTtctc                                                                                                                       <  1:148119/149‑1 (MQ=255)                  tAAATAGCATCCAGGAACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTTCTCAtt                                                                                                                    <  2:284227/149‑1 (MQ=255)                       aGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCa                                                                                                               <  1:790802/149‑1 (MQ=255)                                   caGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGt                                                                                                   >  1:860262/1‑149 (MQ=255)                                      aaGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGAATTTTCTGGTTTCCA‑TTC‑GC‑G‑A‑AA‑ATAGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTcc                                                                                                >  2:915746/1‑149 (MQ=255)                                         gCAAGAAGTAGACAACTTACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTca                                                                                             <  1:555911/149‑1 (MQ=255)                                           aagaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcaca                                                                                           <  1:425432/149‑1 (MQ=255)                                                gTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGTTTTTCTGGTTTCCA‑GTC‑GC‑G‑AAAA‑A‑AGATA‑ACGACTATTTTTCTCATTTTTCAGTACCATATTGTTCCTCACATTcac                                                                                      >  1:480234/1‑149 (MQ=255)                                                 tAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcaca                                                                                     >  2:29557/1‑149 (MQ=255)                                                      aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg                                                                                >  2:829129/1‑149 (MQ=255)                                                      aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCAGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg                                                                                <  1:956975/149‑1 (MQ=255)                                                              gaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGTTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaa                                                                        <  1:137409/149‑1 (MQ=255)                                                                  tCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTaa                                                                            <  2:3855/141‑1 (MQ=255)                                                                  tCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTaa                                                                            >  1:3855/1‑141 (MQ=255)                                                                  tCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaa                                                                    >  2:1021680/1‑149 (MQ=255)                                                                     agagTTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                                 <  2:957586/149‑1 (MQ=255)                                                                       agaTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGATTAAATGCA‑GG‑A‑A‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACa                                                               <  1:763086/149‑1 (MQ=255)                                                                        aaTTACTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCAGGGTTTCTGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTGTCTCATGTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaa                                                              >  2:789338/2‑149 (MQ=255)                                                                         aTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCCCATGGGCACAAGGGTAGTTGAAAGAAAACaaa                                                             >  2:430279/1‑149 (MQ=255)                                                                               tcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa                                                       <  1:953678/149‑1 (MQ=255)                                                                               tcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa                                                       <  1:29557/149‑1 (MQ=255)                                                                                           ttttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCaa                                           >  1:17149/1‑149 (MQ=255)                                                                                           ttttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCaa                                           >  2:144580/1‑149 (MQ=255)                                                                                            tttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAAt                                          >  1:633307/1‑149 (MQ=255)                                                                                            tttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAATCGCAAt                                          >  2:343892/1‑149 (MQ=255)                                                                                               ccTTTCTTCTCCCAGACGGCTGGATTTTCTGGTTTCCA‑CTC‑GC‑C‑A‑AATA‑AGATA‑ACGAATCTTTTTCTCATGTTGCAGTGCCATTTTGTTCCCCACTTGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGc                                       <  1:430279/149‑1 (MQ=255)                                                                                               ccTTTATTGTCCCAGGAGGCCGGATTTTCTGGTTTCCG‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGc                                       <  2:186953/149‑1 (MQ=255)                                                                                                ctttctttTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGcc                                      <  2:930573/149‑1 (MQ=255)                                                                                                 tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                                 <  2:468879/89‑1 (MQ=255)                                                                                                 tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                                 >  1:468879/1‑89 (MQ=255)                                                                                                 tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATg                                                                                         <  1:905852/97‑1 (MQ=255)                                                                                                 tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTC‑GC‑G‑AAAA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATg                                                                                         >  2:905852/1‑97 (MQ=255)                                                                                                  cgagtatCCCAGAGGGCCGGATCTACTAGTTTACA‑GTC‑TCGG‑A‑AA‑A‑AGATA‑AGGACCATTTTTCACATTTTGCAGTGCCATATTGTTCCTCACATGCTAACATTGGAAATGAAAAAAAGACAAGACACGAGGTAAGGCGCAATAGCCAg                                    <  2:994459/143‑1 (MQ=255)                                                                                                          ccAGAAGGGCGGATTTTCTGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGa                            <  1:437501/149‑1 (MQ=255)                                                                                                               aGGCAGGAGTTTGTGGTTGCCA‑GTCTGG‑A‑A‑AA‑A‑AGATA‑ACGAATATTTTTCTCATTTTGCAGTGCCATATTGTGCCTCAGATGCACACATTGGTAATGAAAAAAAGACAAAAGACGAGCTACGGTGCAATAGCCAGTTATTAGAATTaa                       <  2:178432/149‑1 (MQ=255)                                                                                                                ggCCGGATTTTCTGGTTTCCA‑GTC‑GC‑G‑AAAA‑A‑AGATA‑ATCACTATTTTTCTCATTTTGCAGTGCCTTATTGCGCCGCACATCCCGGATTTTGTAATTAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAg                      <  1:1021680/149‑1 (MQ=255)                                                                                                                  ccGGATTTTCTGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAACAAGGAGAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGa                    >  1:808662/1‑149 (MQ=255)                                                                                                                  ccGGATTTTCTGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGa                    >  1:785584/1‑149 (MQ=255)                                                                                                                  acGGAGTTTCTGGTTTCCA‑GTC‑GC‑G‑AAAA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACCTTGGTAATGAAAAAAAGACAAAACACGACGTAAGGCACAATAGCCAATTATTAGAATTAAggc                    >  2:958120/2‑148 (MQ=255)                                                                                                                    ggATTTTCTGGTTTCCA‑GTC‑GC‑G‑A‑AA‑A‑AAAAATCCCCCAATCGTTCGGATGTTGCAGGACCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGACAGTTATTAGAATTAGGGATg                  <  2:480234/149‑1 (MQ=255)                                                                                                                        tttCTGGTTTCCA‑GTC‑GA‑GAA‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAAtt              <  1:781823/149‑1 (MQ=255)                                                                                                                          tCTGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTgg            <  1:829129/149‑1 (MQ=255)                                                                                                                           cTGGGTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTggg           <  2:860262/149‑1 (MQ=255)                                                                                                                            tGGTTTCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGt          >  2:579804/1‑149 (MQ=255)                                                                                                                                ttCCA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTGACATGAAGGATGAATTGGGTGaat      >  2:396046/1‑148 (MQ=255)                                                                                                                                  ccA‑GTC‑GC‑GAA‑AA‑A‑AGATA‑ACGACCATTTTTCTCATTTTGCAGTGCCATATTGTCCCCCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTg    <  1:166493/149‑1 (MQ=255)                                                                                                                                   cA‑GTC‑GC‑G‑AAAA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCCTCCTGCTCACCACATGCCCACATTGGTAATGCAACAAACACAAAACACCAAGTAAGGCGCAATAGCCAGATATTAGAATTAAGGGTGAATTGGGTGAAGtca   >  2:870004/1‑147 (MQ=255)                                                                                                                                    a‑GTC‑GC‑G‑AAAA‑A‑AGATA‑CCGATTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCt  <  1:652091/149‑1 (MQ=255)                                                                                                                                    a‑GTC‑GC‑G‑AAAA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTTGGTGAAGTGCt  <  1:653284/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                 ATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCA‑GTC‑GC‑G‑C‑AA‑A‑AGATA‑ACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCT  >  NC_000913_3_hsa_tpiA/1543884‑1544162

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