Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_hsa_tpiA	4,604,698	C→A	A186E (GCA→GAA)	bglJ →	bgl operon transcriptional activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_hsa_tpiA	4,604,698	0	C	A	100.0%	208.3 / NA	66	A186E (GCA→GAA)	bglJ	bgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (34/32);  total (34/32)

GACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTC  >  NC_000913_3_hsa_tpiA/4604556‑4604843                                                                                                                                               |                                                                                                                                                   gaccgacTGAACAATCAGTGGTACAATAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAAAACCAAAACAATCCGTGAACATaa                                                                                                                                              >  2:713833/1‑149 (MQ=255)  accgacTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATaaa                                                                                                                                             >  2:701183/1‑149 (MQ=255)      acTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGGCGAGCAGCTTAAACGCAATATCAAGACGATCCGGAAACATAAATTTa                                                                                                                                         >  1:685569/1‑149 (MQ=255)      acTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTa                                                                                                                                         >  2:652744/1‑149 (MQ=255)       cTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTaa                                                                                                                                        <  1:162273/149‑1 (MQ=255)          aaCAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAAt                                                                                                                                            <  2:328739/142‑1 (MQ=255)          aaCAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAAt                                                                                                                                            >  1:328739/1‑142 (MQ=255)          aaCAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAAtgt                                                                                                                                     >  2:385726/1‑149 (MQ=255)             aaccattGGAACATTACCCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGAt                                                                                                                                  <  1:510601/143‑1 (MQ=255)               ttaGTGGTACAGTAACCAAAGCCGGACGTTCAGCCCGACGGAGAGAGAAATATTGCGCTTGATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGt                                                                                                                                <  2:829336/147‑1 (MQ=255)                aaGTGGTATATTAACTAATGCCGGACGTTATGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTGCTCAATGTCGCAAATTGTCGTGATGTTTAATTGCAATATCAAAATGATCGGTGATTATAGATTTAATGTGATGTc                                                                                                                               <  2:253969/148‑1 (MQ=255)                    ggTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGaaa                                                                                                                           <  2:337883/149‑1 (MQ=255)                     gTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTGGAAAc                                                                                                                          >  1:213180/1‑149 (MQ=255)                     gTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTGTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAAc                                                                                                                          >  1:688606/1‑149 (MQ=255)                          ttAACCAAAGCCGGACGTTAAGCCCCACGGAGAGAGAAATATTGCGGTTTATGTCGCGTGGCTACTGAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATATATTTAATGTGATGTCGAACCTGGGc                                                                                                                     >  1:895630/1‑149 (MQ=255)                               cAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAg                                                                                                                >  1:300130/1‑149 (MQ=255)                                   gCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTAATGTCGAAACTGGGCGTCAGTT‑Ct                                                                                                           >  1:816162/1‑149 (MQ=255)                                     cGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAACCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGa                                                                                                         <  1:33608/149‑1 (MQ=255)                                     cGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAACCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGa                                                                                                         <  1:33772/149‑1 (MQ=255)                                     cGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGa                                                                                                         <  2:588760/149‑1 (MQ=255)                                       gACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGACg                                                                                                       >  1:18554/1‑149 (MQ=255)                                           tgaaGCCCGACGGGGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGACGCAgg                                                                                                   <  1:314791/147‑1 (MQ=255)                                                      ggAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGAGGCAGGGTTGTTGGAgg                                                                                        >  2:561117/1‑149 (MQ=255)                                                        agagagaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAGTATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGACGCAGGGTTGTTGGAGGcc                                                                                      >  1:508269/1‑149 (MQ=255)                                                        agagagaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGGTGTCGAAACTGGGCGTCAGTT‑CTGACGCAGGGTTGTTGGAGGcc                                                                                      <  2:865968/149‑1 (MQ=255)                                                              aaaTATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGCCAGTT‑CTGACGCAGGGTTGTTGGAAGCCGCAGat                                                                                >  2:823994/1‑149 (MQ=255)                                                               aaTATTGCTTTTTCGGTGGCGTGGCTACTCAATAACACAAATTTCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGACGCAGGGTTGTTGGAGGCCGCAGata                                                                               <  1:678058/149‑1 (MQ=255)                                                                atatTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGtt                                                                                                              <  2:717068/118‑1 (MQ=255)                                                                atatTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGtt                                                                                                              >  1:717068/1‑118 (MQ=255)                                                                   ttGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGACGCAGGGTTGTTGGGGGCCGCCGATAGTCt                                                                           >  1:756497/1‑149 (MQ=255)                                                                       gcTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTATACGCAATATCAAAACGAGCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGAGTCAGTT‑CTGACGCAGGGTTGTTGGAGGGCACACATATTCTGTta                                                                       <  2:895630/149‑1 (MQ=255)                                                                       gcTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑TTGACGCAGGGTTGTTGAAGGCCGCAGATATTCtctca                                                                       >  2:100306/1‑145 (MQ=255)                                                                        cTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGACGCAGGgttgt                                                                                              <  1:671073/125‑1 (MQ=255)                                                                        cTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGACGCAGGgttgt                                                                                              >  2:671073/1‑125 (MQ=255)                                                                             tGTCGCTTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAAAGCAATATCAAAACGATCCGTGAACATAAATTTACTGTGATGTCGAAACTGGGCGTCAGTT‑CTGACGCATGGTTGTTGGAGGCCGCAGATATTCTGTtatgtatg                                                                 <  2:621514/149‑1 (MQ=255)                                                                              gTCGCGTGGTTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGc                                                                <  1:69484/149‑1 (MQ=255)                                                                                         aCTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACCTAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGaa                                                     <  2:107520/149‑1 (MQ=255)                                                                                         aCTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATAAAAACAATGCGTGAACATCAATTCAAAGTGATGTCGAAACTGGGCGTCAGTT‑CTGACGCAGCGTTGTTGGATGCCGCAGATGTTCTGTTTTGTCTGCGGAATTCCgta                                                     >  2:908340/1‑147 (MQ=255)                                                                                           tCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAAc                                                   >  1:36092/1‑149 (MQ=255)                                                                                              aTGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTTGGCGTCAGTT‑CTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAg                                                <  1:385726/149‑1 (MQ=255)                                                                                               tGACAGAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGAGCCGTGAACATAAGTTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGACGCAGGGGTGGTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGt                                               <  2:519917/149‑1 (MQ=255)                                                                                                  cacaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAAt                                            >  2:685536/1‑149 (MQ=255)                                                                                                   acaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAAtg                                           >  1:655330/1‑149 (MQ=255)                                                                                                      aaTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATgtgt                                        <  2:300130/149‑1 (MQ=255)                                                                                                                cagcTTAAACGCAATATCAAAACGGTCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTa                              <  2:32784/149‑1 (MQ=255)                                                                                                                cagcTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTa                              <  2:272021/149‑1 (MQ=255)                                                                                                                   cTTAAAAGCAAGATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCTGAACATGTAGTGTGTTGCATCCCTATTa                           >  2:466757/1‑149 (MQ=255)                                                                                                                    ttAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCTTATaaa                          >  2:684411/1‑149 (MQ=255)                                                                                                                    ttAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTaa                          >  1:413491/1‑149 (MQ=255)                                                                                                                      aaaCGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATc                        >  1:605366/1‑149 (MQ=255)                                                                                                                      aaaCGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGACGCAGGGTTGTTGCAGGCCGCCGATCCTCCGTGGATTATGCGGCATTGCCAAACAAGTAATGTGTTGCATCCCTATTAATc                        >  1:187203/1‑149 (MQ=255)                                                                                                                       aaCGAAAGATCAAAACGAAGCGTGAACATACGTTTAATGTGATGTCGTAACTGGGCGTCAGTT‑CTGACGCAGGGTTGTTGGAGGCTGTATATATTCTGTTATGTATGCGGCATTGCGATACAATTTATGTGTTGCATCCCTATTAATcc                       <  2:853683/149‑1 (MQ=255)                                                                                                                        aCGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGACGCAggg                                                                                                  <  2:365814/73‑1 (MQ=255)                                                                                                                        aCGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGACGCAggg                                                                                                  >  1:365814/1‑73 (MQ=255)                                                                                                                         cGCAATATCAAACCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGc                     <  1:795306/149‑1 (MQ=255)                                                                                                                          gCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCa                    <  1:701183/149‑1 (MQ=255)                                                                                                                             atatCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGACGGAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCatga                 <  2:18554/149‑1 (MQ=255)                                                                                                                                  aaaaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCg            >  2:729952/1‑149 (MQ=255)                                                                                                                                  aaaaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGACCAAGTGGTGGTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCg            >  2:361699/1‑149 (MQ=255)                                                                                                                                   aaaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGAGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCgg           <  1:652744/149‑1 (MQ=255)                                                                                                                                    aaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTACTG‑TGTTGGATTGTTGGACGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCggg          <  1:683460/149‑1 (MQ=255)                                                                                                                                      cGATTCGTGAACATAAATTTAATGAGATGTCGAAACTGGGCGTTAGTT‑TTGACGCAGGGTTGTTGGAGGCCCCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGtt        <  1:522837/149‑1 (MQ=255)                                                                                                                                         tCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTAc     >  1:976413/1‑149 (MQ=255)                                                                                                                                            gTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTa                                              >  1:238323/1‑105 (MQ=255)                                                                                                                                            gTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCAACCCTATTAATCCGCATGATGCCGGGTTTACTTc  >  1:123169/1‑149 (MQ=255)                                                                                                                                            gTGAACATAAATTCAGTGGGCGTTCGAACCTGGGCGACAGTT‑CGGACGCACGGTTGTTGGCGGCCGCAGATATTCTGTTATGTATGCGCCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGACAACAGGTTTACATc  <  2:390946/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                   GACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTT‑CTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTC  >  NC_000913_3_hsa_tpiA/4604556‑4604843

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