Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_hsa_tpiA 389,925 (T)6→7 intergenic (‑198/‑326) hemB ← / → yaiT 5‑aminolevulinate dehydratase (porphobilinogen synthase)/pseudogene, autotransporter family;putative structure; Not classified; interrupted by IS3; putative flagellin structural protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA389,9191.T96.0% 261.7 / 3.7 75intergenic (‑192/‑332)hemB/yaiT5‑aminolevulinate dehydratase (porphobilinogen synthase)/pseudogene, autotransporter family;putative structure; Not classified; interrupted by IS3; putative flagellin structural protein
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base T (34/38);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 6.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTA‑T‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGAAATACT  >  NC_000913_3_hsa_tpiA/389770‑390063
                                                                                                                                                      |                                                                                                                                                  
gAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGttatta                                                                                                                                                     >  3:169422/1‑150 (MQ=255)
       ggTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTA‑t‑t‑ttttt                                                                                                                                           <  3:386386/150‑2 (MQ=255)
        gTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTA‑t‑t‑tttttg                                                                                                                                          <  2:807708/150‑3 (MQ=255)
                  tttACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGtt                                                                                                                                 <  3:95756/150‑1 (MQ=255)
                     aCGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTgcg                                                                                                                              <  2:221189/150‑1 (MQ=255)
                       gAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTTATTTTTTTTTATATCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCgct                                                                                                                            >  2:582610/1‑150 (MQ=255)
                       gAAAGTTGTCGCGATGATGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAGCAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCgct                                                                                                                            >  3:525654/1‑150 (MQ=255)
                        aaaGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCgctg                                                                                                                           <  2:463443/150‑1 (MQ=255)
                        aaaGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCgctg                                                                                                                           <  2:674635/150‑1 (MQ=255)
                          aGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCgctgct                                                                                                                         <  2:332258/150‑1 (MQ=255)
                           gTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCta                                                                                                                        <  1:540801/150‑1 (MQ=255)
                            ttGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCtat                                                                                                                       <  4:452174/150‑1 (MQ=255)
                              gTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATaa                                                                                                                     >  2:569962/1‑150 (MQ=255)
                                 gcgATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAAtt                                                                                                                  >  1:34347/1‑150 (MQ=255)
                                         gACAAGAAGAGAATGGAATAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCAGTTATAATATTTTTAATCTTTATTTTTAGCTGTTAGTTATTATT‑T‑TTTTGGTGTTTGGTTGCGCTGCTATAAATTGAATAAtt                                                                                                          <  4:131574/150‑1 (MQ=255)
                                           cAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGaata                                                                                                             >  3:447344/1‑145 (MQ=255)
                                           cAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGaata                                                                                                             <  4:447344/145‑1 (MQ=255)
                                           cAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTaa                                                                                                        >  3:29606/1‑150 (MQ=255)
                                             agaagaGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTaaaa                                                                                                      <  3:262426/150‑1 (MQ=255)
                                              gaagaGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTgcg                                                                                                                              <  3:507192/125‑1 (MQ=255)
                                              gaagaGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTgcg                                                                                                                              >  4:507192/1‑125 (MQ=255)
                                                    aaTGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTttgtt                                                                                               <  4:151847/150‑1 (MQ=255)
                                                    aaTGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTttgtt                                                                                               <  2:815549/150‑1 (MQ=255)
                                                      tGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGc                                                                                             >  1:62180/1‑150 (MQ=255)
                                                        gAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCat                                                                                           <  1:296888/150‑1 (MQ=255)
                                                        gAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCat                                                                                           >  1:515505/1‑150 (MQ=255)
                                                         aaGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCata                                                                                          >  2:599705/1‑150 (MQ=255)
                                                         aaGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATACACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCatt                                                                                          >  2:571809/1‑149 (MQ=255)
                                                            agagGCCATGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTGTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTc                                                                                       <  2:594024/150‑1 (MQ=255)
                                                               ggCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATa                                                                                    <  2:719217/150‑1 (MQ=255)
                                                                gCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAg                                                                                   >  1:666142/1‑150 (MQ=255)
                                                                   aGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATa                                                                                    <  2:102436/146‑1 (MQ=255)
                                                                   aGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATa                                                                                    >  1:102436/1‑146 (MQ=255)
                                                                     gACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACttt                                                                                                  >  2:459991/1‑130 (MQ=255)
                                                                     gACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACttt                                                                                                  <  1:459991/130‑1 (MQ=255)
                                                                        aTCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTa                                                                           >  2:333700/1‑150 (MQ=255)
                                                                           tttttAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATc                                                                         <  1:372399/149‑1 (MQ=255)
                                                                           tttttAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATc                                                                         >  2:372399/1‑149 (MQ=255)
                                                                           tttttAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCt                                                                        <  3:352562/150‑1 (MQ=255)
                                                                                attAGAAAAAGCAATAGTTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGAAGTTTGGTTGCGTTGCTATATATTGAATAATTAAAATTTTGTTGAATATTCATAGGTTATTTATCTGgttt                                                                   <  2:284062/148‑1 (MQ=255)
                                                                                aaTAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgttt                                                                   <  1:207833/150‑1 (MQ=255)
                                                                                    gAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGTATTTTGTTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgt                                                               >  2:789225/1‑150 (MQ=255)
                                                                                      aaCAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgttt                                                             <  3:504383/150‑1 (MQ=255)
                                                                                      aaCAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgttt                                                             <  3:504554/150‑1 (MQ=255)
                                                                                          gCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGa                                                         >  1:564743/1‑150 (MQ=255)
                                                                                           cAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGat                                                        >  1:518533/1‑150 (MQ=255)
                                                                                            aataatTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTTTTTTTTTTTTACGCGGAGATCGATTGATGAATTTAAACTTTCTTGCTTGTTCATAGGGTATGTATCTGgtttgtgattgagtgt                                                       >  4:187368/1‑145 (MQ=255)
                                                                                             ataatTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGCTTGTGatat                                                      <  1:582610/150‑1 (MQ=255)
                                                                                               aatTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGAtattt                                                    <  3:379258/150‑1 (MQ=255)
                                                                                                  tttATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGAtatttatt                                                 <  4:169422/150‑1 (MQ=255)
                                                                                                   ttATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATttatta                                                >  1:436027/1‑150 (MQ=255)
                                                                                                       atTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCg                                            >  2:604746/1‑150 (MQ=255)
                                                                                                            cTGAAATTAGTTTTATTGTTAAGTTATAGCTGTTGATTATTA‑T‑TGTTTTAGAGTTTGTTTGCGGTGATATAAATTGAATGATTAACAATTTGTTTCATATTCATAGAATACTTATCTGTTTTGTGTTTGTGATATTTATTATGCGTATGc                                       <  2:569485/150‑1 (MQ=255)
                                                                                                            cTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTg                                                                                                                                    >  1:1896/1‑57 (MQ=255)
                                                                                                            cTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTg                                                                                                                                    <  2:1896/57‑1 (MQ=255)
                                                                                                            cTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGc                                       <  1:781200/150‑1 (MQ=255)
                                                                                                            cTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGATATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGc                                       <  1:752722/150‑1 (MQ=255)
                                                                                                             tGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGTATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCt                                      <  4:455672/150‑1 (MQ=255)
                                                                                                             tGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTA‑T‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCtt                                     >  2:546390/1‑150 (MQ=255)
                                                                                                               aaaaTATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCtt                                     <  2:430252/149‑1 (MQ=255)
                                                                                                               aaaaTATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCtt                                     >  1:430252/1‑149 (MQ=255)
                                                                                                               aaaaTATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTc                                    >  2:158535/1‑150 (MQ=255)
                                                                                                                aaaTATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCa                                   >  3:267804/1‑150 (MQ=255)
                                                                                                                      ttttattttttattTATAGCTGTTGGTTATTATT‑T‑TTTTGGATTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACa                             >  4:381383/1‑150 (MQ=255)
                                                                                                                      ttttAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACa                             >  2:200904/1‑150 (MQ=255)
                                                                                                                       tttAATCTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACaa                            <  1:425701/150‑1 (MQ=255)
                                                                                                                       tttAATCTTTATTTATAGCTGTTGGATATTA‑TGT‑TTTTGGAGTTTGTTTGCGCTGCTATACATTGAATATTTAAAACTTTGTTGCATATTCATAGGATATTTATCTTGTTTGTGATTGTGATAGTTATTATTCGTATGCTTCAAAAACaa                            >  4:261988/1‑150 (MQ=255)
                                                                                                                             cTTTATTTATAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTAt                      <  1:447888/150‑1 (MQ=255)
                                                                                                                               ttatttatAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCt                    <  1:191663/150‑1 (MQ=255)
                                                                                                                               ttatttatAGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCt                    >  1:842585/1‑150 (MQ=255)
                                                                                                                                     ataGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACa              >  3:145778/1‑150 (MQ=255)
                                                                                                                                     ataGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACa              >  2:627025/1‑150 (MQ=255)
                                                                                                                                     ataGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACa              >  2:304672/1‑150 (MQ=255)
                                                                                                                                      taGCTGTTGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGGTATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACAt             >  3:41891/1‑150 (MQ=255)
                                                                                                                                          tgtttgttATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAATTTTAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTGGTGATATTTATTATGCGTATGCTTCAAAAGCAAAATTATCTGCAACATGTTg         <  2:62180/150‑1 (MQ=255)
                                                                                                                                             tGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGGTATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGaaa      <  4:41891/150‑1 (MQ=255)
                                                                                                                                             tGGTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGaaa      <  2:842585/150‑1 (MQ=255)
                                                                                                                                              ggTTATTATT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGAAAt     >  2:471711/1‑150 (MQ=255)
                                                                                                                                                 tattatT‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGAAATACt  <  2:176120/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
GAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTA‑T‑T‑TTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGAAATACT  >  NC_000913_3_hsa_tpiA/389770‑390063

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