Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_hsa_tpiA 389,919:1 +T 100% intergenic (‑192/‑332) 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.5% 249.7 / ‑2.7 68intergenic (‑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 (33/34);  total (33/35)
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.39e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

CAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTA‑‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGAA  >  NC_000913_3_hsa_tpiA/389774‑390058
                                                                                                                                                  |                                                                                                                                            
cAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATTGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTA‑‑tttt                                                                                                                                         <  1:230233/150‑1 (MQ=255)
   ggTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTA‑‑ttttttt                                                                                                                                      <  1:91821/150‑2 (MQ=255)
   ggTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTA‑‑ttttttt                                                                                                                                      <  1:23246/150‑2 (MQ=255)
   ggTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTA‑‑ttttttt                                                                                                                                      >  2:350452/1‑149 (MQ=255)
   ggTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTA‑‑ttttttt                                                                                                                                      >  3:64896/1‑149 (MQ=255)
             ttttACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTgtt                                                                                                                             >  2:97436/1‑148 (MQ=255)
             ttttACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGt                                                                                                                             >  2:98539/1‑150 (MQ=255)
              tttACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAATCAGCAATAATTTTATATTCTCTGAAATTATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGtt                                                                                                                            <  2:537582/150‑1 (MQ=255)
               ttACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTg                                                                                                                           >  3:271402/1‑150 (MQ=255)
               ttACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTg                                                                                                                           >  2:181225/1‑150 (MQ=255)
                 aCGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTgcg                                                                                                                         <  2:401176/150‑1 (MQ=255)
                 aCGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTgcg                                                                                                                         <  3:289146/150‑1 (MQ=255)
                  cGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTgcgc                                                                                                                        <  3:193444/150‑1 (MQ=255)
                   gAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCgct                                                                                                                       >  1:460883/1‑150 (MQ=255)
                    aaaGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCgctg                                                                                                                      >  2:271402/1‑150 (MQ=255)
                          gTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATaa                                                                                                                <  1:44584/150‑1 (MQ=255)
                           tCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATaaa                                                                                                               <  3:315283/150‑1 (MQ=255)
                                   ttGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGaataa                                                                                                       >  2:181431/1‑150 (MQ=255)
                                      aCAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTa                                                                                                    >  2:172453/1‑150 (MQ=255)
                                       cAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTaa                                                                                                   <  1:181431/150‑1 (MQ=255)
                                        aagaagAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCCCTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGTTGTTTTTTTGCGCTGCTATAAATTGAATAATTaca                                                                                                  >  2:275571/1‑148 (MQ=255)
                                         agaagaGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTaaaa                                                                                                 <  2:139502/150‑1 (MQ=255)
                                         agaagaGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTaaaa                                                                                                 >  3:386271/1‑150 (MQ=255)
                                         agaagaGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTaaaa                                                                                                 <  4:50633/150‑1 (MQ=255)
                                           aagaGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACt                                                                                               <  1:270384/150‑1 (MQ=255)
                                              agaATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTttg                                                                                            <  1:406802/150‑1 (MQ=255)
                                              agaATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTttg                                                                                            >  2:317987/1‑150 (MQ=255)
                                                aaTGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAATTttgtt                                                                                          >  4:204107/1‑150 (MQ=255)
                                                    gAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCat                                                                                      <  2:534949/150‑1 (MQ=255)
                                                            gCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAAGAATTAAAACTTTGTTGCATATTCATAg                                                                              <  1:420384/150‑1 (MQ=255)
                                                             ccAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAgg                                                                             >  3:108661/1‑150 (MQ=255)
                                                                 gACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTTTTTTTTTTTTTCGCTGCTATAAATTGAATAATTAAAACTTTTTTGCATATTCATAGGatat                                                                         >  2:418946/1‑150 (MQ=255)
                                                                  aCATCTGTTGGATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATAtt                                                                        <  2:168305/150‑1 (MQ=255)
                                                                       tttttAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCt                                                                   >  1:342652/1‑150 (MQ=255)
                                                                         tttAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTgg                                                                 >  1:253296/1‑150 (MQ=255)
                                                                         tttAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTgg                                                                 >  3:230397/1‑150 (MQ=255)
                                                                            aaTAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgttt                                                              <  1:62506/150‑1 (MQ=255)
                                                                              tAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTTTTTTTTTTTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTTTCTTttttgt                                                            >  4:207791/1‑150 (MQ=255)
                                                                                 aaaCAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGTTTTTTGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgtt                                                         >  4:301301/1‑150 (MQ=255)
                                                                                 aaaCAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAGACTTTGTTGCATATTCATAGGATATTTCTCTGgtttgtgtt                                                         >  1:530416/1‑150 (MQ=255)
                                                                                 aaaCAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgtt                                                         >  4:288464/1‑150 (MQ=255)
                                                                                  aaCAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgttt                                                        <  2:475131/150‑1 (MQ=255)
                                                                                     aGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGtttgtgtttgtg                                                     <  1:236918/150‑1 (MQ=255)
                                                                                      gCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGa                                                    <  3:229479/150‑1 (MQ=255)
                                                                                       cAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGat                                                   <  2:184629/150‑1 (MQ=255)
                                                                                       cAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTA‑‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGata                                                  <  3:268478/150‑1 (MQ=255)
                                                                                        aataatTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGata                                                  >  2:506951/1‑150 (MQ=255)
                                                                                          taatTTTATATTCACTGAAAATATTTTTATTCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGAtatt                                                <  4:357418/150‑1 (MQ=255)
                                                                                          taatTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTTTTTTTTTTTTTCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTCGTTTGTGTTTTTTAtatt                                                >  4:251062/1‑150 (MQ=255)
                                                                                          taatTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGAtatt                                                >  3:384839/1‑150 (MQ=255)
                                                                                               ttATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATttatta                                           <  3:136031/150‑1 (MQ=255)
                                                                                                   atTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCg                                       <  1:476748/150‑1 (MQ=255)
                                                                                                      cACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTAt                                    <  4:271402/150‑1 (MQ=255)
                                                                                                      cACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTAt                                    <  1:98539/150‑1 (MQ=255)
                                                                                                         tGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGAGTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCt                                 >  3:379329/1‑150 (MQ=255)
                                                                                                          gAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCtt                                >  1:570741/1‑150 (MQ=255)
                                                                                                               tatTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGTAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGTTATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCaaaa                           >  2:119125/1‑150 (MQ=255)
                                                                                                               tatTTTTAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCaaaa                           <  1:227823/150‑1 (MQ=255)
                                                                                                                   tttAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACaa                       >  2:326727/1‑150 (MQ=255)
                                                                                                                    ttAATCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACaaa                      <  1:575548/150‑1 (MQ=255)
                                                                                                                        tCTTTATTTATAGCTGTTGGTTATTAT‑TTTTTTGGTTTTTTTTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTa                  >  2:503488/1‑150 (MQ=255)
                                                                                                                           ttatttatAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCt               >  3:295520/1‑150 (MQ=255)
                                                                                                                           ttatttatAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCt               >  3:245891/1‑150 (MQ=255)
                                                                                                                           ttatttatAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCt               >  4:66238/1‑150 (MQ=255)
                                                                                                                               ttatAGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCaa           <  2:510921/150‑1 (MQ=255)
                                                                                                                                tataGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAAc          <  4:26393/150‑1 (MQ=255)
                                                                                                                                 ataGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACa         <  2:290523/150‑1 (MQ=255)
                                                                                                                                   aGCTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCATATGCTTAAAAAACAAAATTATCTGCAACATg       <  2:166359/150‑1 (MQ=255)
                                                                                                                                    gCTGTTGGTTATTATTTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATg       <  2:14831/150‑1 (MQ=255)
                                                                                                                                     cTGTTGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATTTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGtt     <  3:196659/150‑1 (MQ=255)
                                                                                                                                       ttttgTTATTAT‑TTTTTTGGAGTTTATTTTCTTTTCTATAAATTAAATAATTAAAACTTTTTTGCATATTCATAAGATATTTATCTGCTTTGTGTTTGTGATATTTATTATGAAAAAAAAAAAAAAACAAAATTATCTGCAACATGTTGa   <  4:108661/146‑1 (MQ=255)
                                                                                                                                        ttGGTTATTAT‑TTTTTTTTTTTTTTTTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTTTTTGTTTTTTTCAAAAACAAAATTATCTGCAACATGTTGaa  >  4:251124/1‑150 (MQ=255)
                                                                                                                                        ttGGTTATTAT‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGaa  <  1:97436/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
CAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTA‑‑TTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGAA  >  NC_000913_3_hsa_tpiA/389774‑390058

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