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

GAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGAAATACT  >  NC_000913_3_hsa_tpiA/389770‑390063
                                                                                                                                                      |                                                                                                                                                  
gAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGttatta                                                                                                                                                     >  3:169422/1‑150 (MQ=255)
       ggTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAttttttt                                                                                                                                           <  3:386386/150‑2 (MQ=255)
        gTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAtttttttg                                                                                                                                          <  2:807708/150‑3 (MQ=255)
                  tttACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGtt                                                                                                                                 <  3:95756/150‑1 (MQ=255)
                     aCGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTgcg                                                                                                                              <  2:221189/150‑1 (MQ=255)
                       gAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTTATTTTTTTTTATATCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCgct                                                                                                                            >  2:582610/1‑150 (MQ=255)
                       gAAAGTTGTCGCGATGATGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAGCAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCgct                                                                                                                            >  3:525654/1‑150 (MQ=255)
                        aaaGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCgctg                                                                                                                           <  2:463443/150‑1 (MQ=255)
                        aaaGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCgctg                                                                                                                           <  2:674635/150‑1 (MQ=255)
                          aGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCgctgct                                                                                                                         <  2:332258/150‑1 (MQ=255)
                           gTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCta                                                                                                                        <  1:540801/150‑1 (MQ=255)
                            ttGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCtat                                                                                                                       <  4:452174/150‑1 (MQ=255)
                              gTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATaa                                                                                                                     >  2:569962/1‑150 (MQ=255)
                                 gcgATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAAtt                                                                                                                  >  1:34347/1‑150 (MQ=255)
                                         gACAAGAAGAGAATGGAATAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCAGTTATAATATTTTTAATCTTTATTTTTAGCTGTTAGTTATTATTTTTTTGGTGTTTGGTTGCGCTGCTATAAATTGAATAAtt                                                                                                          <  4:131574/150‑1 (MQ=255)
                                           cAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGaata                                                                                                             >  3:447344/1‑145 (MQ=255)
                                           cAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGaata                                                                                                             <  4:447344/145‑1 (MQ=255)
                                           cAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTaa                                                                                                        >  3:29606/1‑150 (MQ=255)
                                             agaagaGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTaaaa                                                                                                      <  3:262426/150‑1 (MQ=255)
                                              gaagaGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTgcg                                                                                                                              <  3:507192/125‑1 (MQ=255)
                                              gaagaGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTgcg                                                                                                                              >  4:507192/1‑125 (MQ=255)
                                                    aaTGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTttgtt                                                                                               <  4:151847/150‑1 (MQ=255)
                                                    aaTGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTttgtt                                                                                               <  2:815549/150‑1 (MQ=255)
                                                      tGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGc                                                                                             >  1:62180/1‑150 (MQ=255)
                                                        gAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCat                                                                                           <  1:296888/150‑1 (MQ=255)
                                                        gAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCat                                                                                           >  1:515505/1‑150 (MQ=255)
                                                         aaGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCata                                                                                          >  2:599705/1‑150 (MQ=255)
                                                         aaGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATACACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCatt                                                                                          >  2:571809/1‑149 (MQ=255)
                                                            agagGCCATGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTGTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTc                                                                                       <  2:594024/150‑1 (MQ=255)
                                                               ggCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATa                                                                                    <  2:719217/150‑1 (MQ=255)
                                                                gCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAg                                                                                   >  1:666142/1‑150 (MQ=255)
                                                                   aGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATa                                                                                    <  2:102436/146‑1 (MQ=255)
                                                                   aGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATa                                                                                    >  1:102436/1‑146 (MQ=255)
                                                                     gACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACttt                                                                                                  >  2:459991/1‑130 (MQ=255)
                                                                     gACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACttt                                                                                                  <  1:459991/130‑1 (MQ=255)
                                                                        aTCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTa                                                                           >  2:333700/1‑150 (MQ=255)
                                                                           tttttAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATc                                                                         <  1:372399/149‑1 (MQ=255)
                                                                           tttttAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATc                                                                         >  2:372399/1‑149 (MQ=255)
                                                                           tttttAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCt                                                                        <  3:352562/150‑1 (MQ=255)
                                                                                attAGAAAAAGCAATAGTTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGAAGTTTGGTTGCGTTGCTATATATTGAATAATTAAAATTTTGTTGAATATTCATAGGTTATTTATCTGgttt                                                                   <  2:284062/148‑1 (MQ=255)
                                                                                aaTAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgttt                                                                   <  1:207833/150‑1 (MQ=255)
                                                                                    gAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGTATTTTGTTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgt                                                               >  2:789225/1‑150 (MQ=255)
                                                                                      aaCAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgttt                                                             <  3:504383/150‑1 (MQ=255)
                                                                                      aaCAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgttt                                                             <  3:504554/150‑1 (MQ=255)
                                                                                          gCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGa                                                         >  1:564743/1‑150 (MQ=255)
                                                                                           cAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGat                                                        >  1:518533/1‑150 (MQ=255)
                                                                                            aataatTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTTTTTTTTTTTTACGCGGAGATCGATTGATGAATTTAAACTTTCTTGCTTGTTCATAGGGTATGTATCTGgtttgtgattgagtgt                                                       >  4:187368/1‑145 (MQ=255)
                                                                                             ataatTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGCTTGTGatat                                                      <  1:582610/150‑1 (MQ=255)
                                                                                               aatTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGAtattt                                                    <  3:379258/150‑1 (MQ=255)
                                                                                                  tttATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGAtatttatt                                                 <  4:169422/150‑1 (MQ=255)
                                                                                                   ttATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATttatta                                                >  1:436027/1‑150 (MQ=255)
                                                                                                       atTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCg                                            >  2:604746/1‑150 (MQ=255)
                                                                                                            cTGAAATTAGTTTTATTGTTAAGTTATAGCTGTTGATTATTATTGTTTTAGAGTTTGTTTGCGGTGATATAAATTGAATGATTAACAATTTGTTTCATATTCATAGAATACTTATCTGTTTTGTGTTTGTGATATTTATTATGCGTATGc                                       <  2:569485/150‑1 (MQ=255)
                                                                                                            cTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTg                                                                                                                                    >  1:1896/1‑57 (MQ=255)
                                                                                                            cTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTg                                                                                                                                    <  2:1896/57‑1 (MQ=255)
                                                                                                            cTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGc                                       <  1:781200/150‑1 (MQ=255)
                                                                                                            cTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGATATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGc                                       <  1:752722/150‑1 (MQ=255)
                                                                                                             tGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGTATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCt                                      <  4:455672/150‑1 (MQ=255)
                                                                                                             tGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCtt                                     >  2:546390/1‑150 (MQ=255)
                                                                                                               aaaaTATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCtt                                     <  2:430252/149‑1 (MQ=255)
                                                                                                               aaaaTATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCtt                                     >  1:430252/1‑149 (MQ=255)
                                                                                                               aaaaTATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTc                                    >  2:158535/1‑150 (MQ=255)
                                                                                                                aaaTATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCa                                   >  3:267804/1‑150 (MQ=255)
                                                                                                                      ttttattttttattTATAGCTGTTGGTTATTATTTTTTTGGATTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACa                             >  4:381383/1‑150 (MQ=255)
                                                                                                                      ttttAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACa                             >  2:200904/1‑150 (MQ=255)
                                                                                                                       tttAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACaa                            <  1:425701/150‑1 (MQ=255)
                                                                                                                       tttAATCTTTATTTATAGCTGTTGGATATTATGTTTTTGGAGTTTGTTTGCGCTGCTATACATTGAATATTTAAAACTTTGTTGCATATTCATAGGATATTTATCTTGTTTGTGATTGTGATAGTTATTATTCGTATGCTTCAAAAACaa                            >  4:261988/1‑150 (MQ=255)
                                                                                                                             cTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTAt                      <  1:447888/150‑1 (MQ=255)
                                                                                                                               ttatttatAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCt                    <  1:191663/150‑1 (MQ=255)
                                                                                                                               ttatttatAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCt                    >  1:842585/1‑150 (MQ=255)
                                                                                                                                     ataGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACa              >  3:145778/1‑150 (MQ=255)
                                                                                                                                     ataGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACa              >  2:627025/1‑150 (MQ=255)
                                                                                                                                     ataGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACa              >  2:304672/1‑150 (MQ=255)
                                                                                                                                      taGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGGTATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACAt             >  3:41891/1‑150 (MQ=255)
                                                                                                                                          tgtttgttATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAATTTTAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTGGTGATATTTATTATGCGTATGCTTCAAAAGCAAAATTATCTGCAACATGTTg         <  2:62180/150‑1 (MQ=255)
                                                                                                                                             tGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGGTATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGaaa      <  4:41891/150‑1 (MQ=255)
                                                                                                                                             tGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGaaa      <  2:842585/150‑1 (MQ=255)
                                                                                                                                              ggTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGAAAt     >  2:471711/1‑150 (MQ=255)
                                                                                                                                                 tattatTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGAAATACt  <  2:176120/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: