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.T92.6% 407.7 / 25.4 121intergenic (‑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 . (3/6);  new base T (65/47);  total (68/53)
Fisher's exact test for biased strand distribution p-value = 1.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 131 total aligned reads displayed.

GAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGAAATACT  >  NC_000913_3_hsa_tpiA/389770‑390063
                                                                                                                                                      |                                                                                                                                                  
gAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGttatta                                                                                                                                                     >  2:214817/1‑150 (MQ=255)
     aGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAttttt                                                                                                                                              >  1:23917/1‑150 (MQ=255)
      gggTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAtttttt                                                                                                                                             <  2:24598/150‑1 (MQ=255)
       ggTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAttttttt                                                                                                                                            >  3:779091/1‑149 (MQ=255)
             gATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGtt                                                                                                                                      >  3:104943/1‑150 (MQ=255)
             gATGTTTGACGAAAGTTGTTGCGATATTGATAAGAAGAGAATGTAAGAGAGGCCAGAACATCTTTTTAATAGATACAAAAATAATTTTGTATTCAGTGAAAATATTTTTAATCTTTATTTATAGCTGTTGCTTATTATCTTTTTCGAGtt                                                                                                                                      <  4:161948/150‑1 (MQ=255)
               ggggtgaCGAAAGTGGTCGGGATGTTGACAAGAAGAGAAGGGAAGAGGGGCCAGGAGATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTgg                                                                                                                                   <  2:985150/144‑1 (MQ=255)
                 ttttACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGt                                                                                                                                  >  2:81742/1‑150 (MQ=255)
                   ttACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTg                                                                                                                                <  2:1069085/150‑1 (MQ=255)
                     aCGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTgcg                                                                                                                              >  2:1141202/1‑150 (MQ=255)
                       gAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTTCgct                                                                                                                            >  1:318221/1‑150 (MQ=255)
                        aaaGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCgctg                                                                                                                           >  1:802938/1‑150 (MQ=255)
                         aaGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCgctgc                                                                                                                          <  4:587527/150‑1 (MQ=255)
                           gTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCta                                                                                                                        >  2:460711/1‑150 (MQ=255)
                               tCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATaaa                                                                                                                    >  2:79017/1‑150 (MQ=255)
                                  cgATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTg                                                                                                                 <  1:570804/150‑1 (MQ=255)
                                   gATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGa                                                                                                                <  3:209915/150‑1 (MQ=255)
                                    aTGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGaa                                                                                                               <  3:23470/150‑1 (MQ=255)
                                    aTGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGaa                                                                                                               <  2:1092800/150‑1 (MQ=255)
                                          aCAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTa                                                                                                         >  1:892667/1‑150 (MQ=255)
                                            aagaagAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTaaa                                                                                                       >  3:387384/1‑150 (MQ=255)
                                            aagaagAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTaaa                                                                                                       >  2:125827/1‑150 (MQ=255)
                                              gaagaGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAAc                                                                                                     <  2:802938/150‑1 (MQ=255)
                                                 gagaATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACttt                                                                                                  <  2:892667/150‑1 (MQ=255)
                                                  agaATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCCCTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGATATAAATTGAATAATTAAAACTttg                                                                                                 >  2:1019830/1‑150 (MQ=255)
                                                    aaTGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTttgtt                                                                                               <  4:486596/150‑1 (MQ=255)
                                                    aaTGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTttgtt                                                                                               >  1:932321/1‑150 (MQ=255)
                                                     aTGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTttgttg                                                                                              >  2:59952/1‑150 (MQ=255)
                                                       ggAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCa                                                                                            >  3:479331/1‑150 (MQ=255)
                                                        gAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTTTTTTTTTTTTACAACAGCAGTTAATTATTTTTTAGGCGTTTGGTCGCTCTAGTCTAAATCGAATAATTTCACCTTCGTTGCat                                                                                           >  2:440780/1‑150 (MQ=255)
                                                        gAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCat                                                                                           >  3:193059/1‑150 (MQ=255)
                                                         aaGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAttttttt                                                                                                                                            >  3:508616/1‑99 (MQ=255)
                                                         aaGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAttttttt                                                                                                                                            <  4:508616/100‑2 (MQ=255)
                                                           gagagGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACttt                                                                                                  <  1:226866/140‑1 (MQ=255)
                                                           gagagGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACttt                                                                                                  >  2:226866/1‑140 (MQ=255)
                                                                gCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAg                                                                                   <  4:779091/150‑1 (MQ=255)
                                                                gCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAg                                                                                   >  3:483607/1‑150 (MQ=255)
                                                                        aTCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAATTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTAt                                                                          <  3:221519/150‑1 (MQ=255)
                                                                        aTCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGATGTTGGTTATTATTTTTTTGGAGTTTGGTTGCTATGCTATAAATTGAATAATTAAAACTTTGTTGCATATTAATAGGATAtttt                                                                           >  4:712042/1‑149 (MQ=255)
                                                                         tCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTAt                                                                          >  4:262333/1‑150 (MQ=255)
                                                                            ttttAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTg                                                                       >  1:759709/1‑150 (MQ=255)
                                                                                aaTAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgttt                                                                   <  3:429411/150‑1 (MQ=255)
                                                                                  tAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgt                                                                 <  2:908531/150‑1 (MQ=255)
                                                                                  tAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgt                                                                 >  1:750670/1‑150 (MQ=255)
                                                                                  tAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtg                                                                <  1:179014/150‑1 (MQ=255)
                                                                                   aGAAACAGCAATATTTTTCTATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGATTTTGTTTTCTCTGCTATAAATTTAATATTTAAAATTTTTTTGCATCTTTATAGGATATTTTTCTGgtttgtg                                                                >  4:507879/1‑150 (MQ=255)
                                                                                   aGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtg                                                                >  1:258850/1‑150 (MQ=255)
                                                                                      aaCAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgttt                                                             <  2:940542/150‑1 (MQ=255)
                                                                                        cAGCAATAATTTTATATTCACTGTAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgtttgt                                                           <  4:106725/150‑1 (MQ=255)
                                                                                        cAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgtttgt                                                           <  1:125827/150‑1 (MQ=255)
                                                                                          gCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTTTGa                                                         >  4:736592/1‑150 (MQ=255)
                                                                                           cAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGat                                                        >  1:15898/1‑150 (MQ=255)
                                                                                           cAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGat                                                        >  2:125799/1‑150 (MQ=255)
                                                                                            aataatTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTTTTTTTTTTTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGTTTTGTGTTTGTGatt                                                       >  4:328006/1‑149 (MQ=255)
                                                                                            aataatTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGata                                                       >  2:225672/1‑150 (MQ=255)
                                                                                            aataatTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTGGGTTATTATTTTTTTGGAGTTTTGTTGCGCTGCTAGAAATTGAATAATTAAAACTTTGGTGCATATGCATAGGAGATTTATTGGGTTTGTGTTTGTGatt                                                       >  1:319020/1‑149 (MQ=255)
                                                                                                 ttttATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGAtatttat                                                  >  4:433567/1‑150 (MQ=255)
                                                                                                 ttttATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGAtatttat                                                  <  4:390017/150‑1 (MQ=255)
                                                                                                   ttATATTCACTGAAAATATTTTTAATCTTTATTGATAGCTGTTGGTAATTATTTTTTTGTGGTTTGGTTGCGCTGCTATTAATTGAATAATTAAAGCTTTGTTGCATATTCATATGATATTTGTCTGGTTTGTGGTTGTAATATttatta                                                >  4:117371/1‑150 (MQ=255)
                                                                                                    tataTTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTtattat                                               <  1:1171187/150‑1 (MQ=255)
                                                                                                     atatTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATg                                              <  2:750670/150‑1 (MQ=255)
                                                                                                         tCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTTTTTTTTTGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATTCGtt                                          >  4:61386/1‑149 (MQ=255)
                                                                                                         tCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTa                                          >  1:219404/1‑150 (MQ=255)
                                                                                                           aCTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATg                                        <  4:483607/150‑1 (MQ=255)
                                                                                                             tGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTTTTTTTTTTTTTCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTTTTTTTTGTTTGTGATATTTTTtattcttattct                                      >  4:598949/1‑150 (MQ=255)
                                                                                                             tGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGTTTTTTGTTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATTCt                                      >  4:150455/1‑150 (MQ=255)
                                                                                                              gAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGTGTTTGGTTGTGCTGCTAAAAATTGAATAATTAAAACTTTGTTGCATATTCATACGTTGTTTATCTGGTTTGTGTTTGTGTTATTTATTATACGTATGAtt                                     >  4:696377/1‑150 (MQ=255)
                                                                                                               aaaaTATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTc                                    <  1:65916/150‑1 (MQ=255)
                                                                                                               aaaaTATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTc                                    >  2:815391/1‑150 (MQ=255)
                                                                                                                aaaTATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCa                                   <  1:1010007/150‑1 (MQ=255)
                                                                                                                aaaTATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCaa                                  >  1:1077736/1‑150 (MQ=255)
                                                                                                                    atTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCaaaaa                               <  3:194722/150‑1 (MQ=255)
                                                                                                                    atTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATCAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCaaaaa                               >  3:332506/1‑150 (MQ=255)
                                                                                                                     tttttAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAAc                              <  3:504166/150‑1 (MQ=255)
                                                                                                                      ttttAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACa                             >  4:19264/1‑150 (MQ=255)
                                                                                                                      ttgtAATCTTTATTTATTGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTATGTGATATTTATTATGCGTATGCTTCAAAAACa                             <  4:558943/147‑1 (MQ=255)
                                                                                                                       tttAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACaa                            <  3:393884/150‑1 (MQ=255)
                                                                                                                       tttAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTaaaaaaaaaa                           <  4:303338/150‑1 (MQ=255)
                                                                                                                        ttAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgt                                                                 <  1:371800/112‑1 (MQ=255)
                                                                                                                        ttAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgt                                                                 >  2:371800/1‑112 (MQ=255)
                                                                                                                        ttAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACaaa                           <  4:466108/150‑1 (MQ=255)
                                                                                                                        ttAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACaaa                           <  2:458477/150‑1 (MQ=255)
                                                                                                                             cTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGAtattt                                                    >  1:1221604/1‑120 (MQ=255)
                                                                                                                             cTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGAtattt                                                    <  2:1221604/120‑1 (MQ=255)
                                                                                                                             cTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTAt                      >  3:663222/1‑150 (MQ=255)
                                                                                                                             cTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTAt                      <  3:19264/150‑1 (MQ=255)
                                                                                                                              tttatttaTAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATc                     >  1:325633/1‑150 (MQ=255)
                                                                                                                               ttatttatAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCt                    >  4:323339/1‑150 (MQ=255)
                                                                                                                               ttatttatAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCt                    >  1:77962/1‑150 (MQ=255)
                                                                                                                                   ttatAGCTGTTGGTCATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCaa                <  1:242266/150‑1 (MQ=255)
                                                                                                                                       aGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATg            <  3:580615/150‑1 (MQ=255)
                                                                                                                                        gCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGt           <  3:193178/150‑1 (MQ=255)
                                                                                                                                         cTGTTGGTTATTATTTTTTTTGTTTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGtt          >  2:243709/1‑150 (MQ=255)
                                                                                                                                         cTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAATCTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAATAACGAAATTATCTGCAACATGtt          <  4:772460/150‑1 (MQ=255)
                                                                                                                                          tGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTg         <  2:77962/150‑1 (MQ=255)
                                                                                                                                          tGTTGGTTATTATTTTTTTGGAGTTTGGTCGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGAATATTCATAGGGTATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTg         <  4:196771/150‑1 (MQ=255)
                                                                                                                                            ttGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGAAAAAAAAAAAAAAACAAAATTATCTGCAACATGTTGaa       <  2:915899/150‑1 (MQ=255)
                                                                                                                                               gTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGAAATa    >  1:572471/1‑150 (MQ=255)
                                                                                                                                                ttattatTTTTTTGGAGTTTGGTTGCGCCGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATTTGGTTGGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGGTGAAATAc   >  4:775361/1‑150 (MQ=255)
                                                                                                                                                 tattatTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGAAATACt  <  1:225672/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
GAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGAAATACT  >  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: