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.T98.7% 279.6 / ‑3.0 78intergenic (‑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 . (0/1);  new base T (37/40);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTA‑T‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGAAATAC  >  NC_000913_3_hsa_tpiA/389773‑390062
                                                                                                                                                   |                                                                                                                                                
gCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGAAATCTTTTTAATCGAAACAGCACTAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTCTAGCTGTTGGTTATTA‑t‑tt                                                                                                                                              >  1:1353615/1‑150 (MQ=255)
 cAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTA‑t‑ttt                                                                                                                                             <  3:505086/150‑1 (MQ=255)
    ggTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTA‑t‑tttttt                                                                                                                                          <  2:866419/150‑2 (MQ=255)
    ggTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTA‑t‑tttttt                                                                                                                                          <  2:244236/150‑2 (MQ=255)
         ggATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGt                                                                                                                                      >  4:317553/1‑150 (MQ=255)
         ggATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGt                                                                                                                                      >  1:1189107/1‑150 (MQ=255)
          gATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGtt                                                                                                                                     <  1:1171440/150‑1 (MQ=255)
                  aCGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTgcg                                                                                                                             >  1:1168870/1‑150 (MQ=255)
                    gAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCgct                                                                                                                           >  4:825046/1‑150 (MQ=255)
                    gAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCgct                                                                                                                           >  4:706755/1‑150 (MQ=255)
                     aaaGTTGTCGCGATGTTGACAAGAAGAGAATGGGAAAGAGGCCAGGATAAATTTTTAATATAAACAGAAATAATTTTATATTCGTTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCgctg                                                                                                                          <  3:243118/150‑1 (MQ=255)
                       aGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCgctgct                                                                                                                        <  4:645043/150‑1 (MQ=255)
                             cgcgATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAAt                                                                                                                  <  3:825046/150‑1 (MQ=255)
                              gcgATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGttt                                                                                                                                    >  1:47970/1‑131 (MQ=255)
                              gcgATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGttt                                                                                                                                    <  2:47970/131‑1 (MQ=255)
                                 aTGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGaa                                                                                                              <  2:197845/150‑1 (MQ=255)
                                 aTGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGaa                                                                                                              <  4:357519/150‑1 (MQ=255)
                                   gTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGaata                                                                                                            <  3:706755/150‑1 (MQ=255)
                                    ttGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGaataa                                                                                                           <  4:180746/150‑1 (MQ=255)
                                      gACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAAtt                                                                                                         >  4:104536/1‑150 (MQ=255)
                                      gACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAAtt                                                                                                         <  2:545027/150‑1 (MQ=255)
                                       aCAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTa                                                                                                        >  2:992540/1‑150 (MQ=255)
                                         aagaagAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTaaa                                                                                                      <  1:1172844/150‑1 (MQ=255)
                                           gaagaGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTTTCTTTTGATTGCGCTGCTGTTAATTGATGAATTAAAAc                                                                                                    >  4:455650/1‑150 (MQ=255)
                                           gaagaGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGaa                                                                                                              <  1:72031/140‑1 (MQ=255)
                                           gaagaGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGaa                                                                                                              >  2:72031/1‑140 (MQ=255)
                                              gagaATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCgctgc                                                                                                                         <  1:627813/126‑1 (MQ=255)
                                              gagaATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCgctgc                                                                                                                         >  2:627813/1‑126 (MQ=255)
                                              gagaATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACttt                                                                                                 <  2:1333593/150‑1 (MQ=255)
                                               agaATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTCAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTttg                                                                                                >  1:1161852/1‑150 (MQ=255)
                                                 aaTGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTttgtt                                                                                              >  1:1140285/1‑150 (MQ=255)
                                                 aaTGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTATATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTttctt                                                                                              >  4:627524/1‑150 (MQ=255)
                                                  aTGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCgctgc                                                                                                                         <  4:570112/122‑1 (MQ=255)
                                                  aTGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCgctgc                                                                                                                         >  3:570112/1‑122 (MQ=255)
                                                   tGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGc                                                                                            >  4:712770/1‑150 (MQ=255)
                                                   tGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGc                                                                                            <  2:250793/150‑1 (MQ=255)
                                                   tGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGc                                                                                            <  1:969692/150‑1 (MQ=255)
                                                   tGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGc                                                                                            <  3:512213/150‑1 (MQ=255)
                                                    ggAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCa                                                                                           <  2:1091544/150‑1 (MQ=255)
                                                         agagGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTc                                                                                      <  2:764813/150‑1 (MQ=255)
                                                         agagGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTc                                                                                      >  1:929970/1‑150 (MQ=255)
                                                            ggCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATa                                                                                   <  1:1240245/150‑1 (MQ=255)
                                                            ggCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATa                                                                                   <  1:189421/150‑1 (MQ=255)
                                                            ggCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATa                                                                                   <  3:616005/150‑1 (MQ=255)
                                                              ccAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAgg                                                                                 >  1:993272/1‑150 (MQ=255)
                                                                 ggACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGata                                                                              <  2:929970/150‑1 (MQ=255)
                                                                    cATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATAttt                                                                           <  2:775767/150‑1 (MQ=255)
                                                                    cATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATAttt                                                                           >  4:463785/1‑150 (MQ=255)
                                                                     aTCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTTTTTTTTTTTTGCGATGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTa                                                                          >  2:489743/1‑150 (MQ=255)
                                                                     aTCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTa                                                                          >  2:1009961/1‑150 (MQ=255)
                                                                         ttttAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAg                                                                                  <  2:411641/138‑1 (MQ=255)
                                                                         ttttAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAg                                                                                  >  1:411641/1‑138 (MQ=255)
                                                                             aaTAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgttt                                                                  <  1:1008977/150‑1 (MQ=255)
                                                                               tAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTTTTTTTTTTTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATATGTTTTTTATCTGgtttgt                                                                >  2:1089523/1‑150 (MQ=255)
                                                                                 gAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgt                                                              >  3:479261/1‑150 (MQ=255)
                                                                                   aaCAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGACTAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgttt                                                            <  1:992540/150‑1 (MQ=255)
                                                                                     cAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgtttgt                                                          <  4:479261/150‑1 (MQ=255)
                                                                                     cAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgtttgt                                                          <  3:104536/150‑1 (MQ=255)
                                                                                       gCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTTTTTTTTGTTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTTTTTGTGa                                                        >  4:704697/1‑150 (MQ=255)
                                                                                       gCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGc                                                                                            >  2:921254/1‑114 (MQ=255)
                                                                                       gCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGc                                                                                            <  1:921254/114‑1 (MQ=255)
                                                                                       gCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGa                                                        >  4:685283/1‑150 (MQ=255)
                                                                                        cAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTA‑TATTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGat                                                       <  1:609074/150‑1 (MQ=255)
                                                                                              ttttATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGAtattt                                                   <  3:672392/148‑1 (MQ=255)
                                                                                              ttttATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGAtattt                                                   >  4:672392/1‑148 (MQ=255)
                                                                                                 tataTTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTtattat                                              <  4:128828/150‑1 (MQ=255)
                                                                                                       cACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTAt                                        <  3:463785/150‑1 (MQ=255)
                                                                                                        aCTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATg                                       >  4:621910/1‑150 (MQ=255)
                                                                                                         cTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgt                                                              >  1:245734/1‑126 (MQ=255)
                                                                                                         cTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgt                                                              <  2:245734/126‑1 (MQ=255)
                                                                                                         cTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGc                                      <  2:392994/150‑1 (MQ=255)
                                                                                                               atatTTTTAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCaaa                                <  2:1189107/150‑1 (MQ=255)
                                                                                                                   aattAATCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGTTTTTTGGTTGCGTTGGTATAAATTGAATATTTATCAGTTTGTTGTAGATTAGTGGGAAATTTCACTGCATTGTGTTTGTTCTGTTGTTTATGCATATTCTTCAAAAACa                            >  2:1310492/3‑150 (MQ=255)
                                                                                                                        aTCTTTATTTATAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAAtt                       >  3:213843/1‑150 (MQ=255)
                                                                                                                           tttatttaTAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgt                                                              >  3:642993/1‑108 (MQ=255)
                                                                                                                           tttatttaTAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgt                                                              <  4:642993/108‑1 (MQ=255)
                                                                                                                           tttatttaTAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATc                    >  4:353722/1‑150 (MQ=255)
                                                                                                                             tatttatAGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTg                  <  2:1245238/150‑1 (MQ=255)
                                                                                                                               tttatAGCTGTTGGTTATTATT‑TTTTTTTTTTTTTTTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCTTATGCTTTAAAAACAATATTATCTGCa                >  2:462178/1‑150 (MQ=255)
                                                                                                                                  ataGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTTCTATAAATTGAATAATTAAAACTTTGTTACATATTCATAGGATAttt                                                                           <  2:839098/88‑1 (MQ=255)
                                                                                                                                  ataGCTGTTGGTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATAttt                                                                           >  1:839098/1‑88 (MQ=255)
                                                                                                                                            gTTATTATT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGAAATa   >  3:540995/1‑150 (MQ=255)
                                                                                                                                             ttattatT‑TTTTTGGAGTTTGGTTTCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGAAATAc  >  3:763460/1‑150 (MQ=255)
                                                                                                                                             ttattatT‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGAAATAc  <  1:1009961/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTA‑T‑TTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGAAATAC  >  NC_000913_3_hsa_tpiA/389773‑390062

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