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.T83.3% 87.6 / 13.8 36intergenic (‑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/4);  new base T (16/14);  total (18/18)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

GTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑A‑T‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTT  >  NC_000913_3_hsa_tpiA/389778‑390055
                                                                                                                                                |                                                                                                                                            
gTAGGGATGTTTTATGAAAGTTGTCGCGATGTTTACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAAGAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑A‑t‑t‑ttttt                                                                                                                                     <  2:361346/149‑2 (MQ=255)
         ttttACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTgg                                                                                                                            >  2:242151/1‑149 (MQ=255)
         ttttACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGCTATT‑ATT‑T‑TTTT‑GGGGTTTgt                                                                                                                            >  2:10265/1‑148 (MQ=255)
              cGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTgcg                                                                                                                       >  1:103646/1‑149 (MQ=255)
                aaaGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGTTTGCgct                                                                                                                     >  2:13011/1‑149 (MQ=255)
                                 gACAAGAAGAGAATGGAAGAGAGGCTAGGACAGCTTTTTAATATAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGaatcac                                                                                                    <  2:173356/149‑4 (MQ=255)
                                     agaagaGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTaaa                                                                                               >  2:517736/1‑149 (MQ=255)
                                         gagaATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACtt                                                                                           >  1:63292/1‑149 (MQ=255)
                                            aaTGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTttgt                                                                                        >  2:236034/1‑149 (MQ=255)
                                                gAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑A‑T‑T‑TTTTCGGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCa                                                                                    <  1:581908/149‑1 (MQ=255)
                                                         ccAGGACATATATTTAATAGAAAAATCAATAATTTTAGATTCACTGAAAATATTTTTAATCTTCATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCaccc                                                                           <  2:611149/149‑4 (MQ=255)
                                                           aGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGa                                                                         >  1:440644/1‑149 (MQ=255)
                                                           aGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGa                                                                         >  1:97803/1‑149 (MQ=255)
                                                                      ttAATAGTAACAGCAATAATTATATATGCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑TGTTATT‑ATT‑T‑TTTT‑GGAGTTTGCTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTgg                                                              <  2:309983/149‑1 (MQ=255)
                                                                       tAATAGAGACATCAATAATTTTATATTCACTGAAATTATTTTTAATCTTTATTTATAGCTGTT‑TGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCACTGCTATAAATTTAATCAATCACAAAGAT‑CGGTATATTCATAGGATATTTATCTGgt                                                             <  2:305717/149‑1 (MQ=255)
                                                                        aaTAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTGTAGCTGTG‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAGACTTTGTTGCATATTCATAGGATATTTCTGTGgtt                                                            >  1:22429/1‑149 (MQ=255)
                                                                         aTAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTAGATCTTT‑T‑T‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgttt                                                           >  1:125384/1‑149 (MQ=255)
                                                                         aTAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgttt                                                           >  1:434651/1‑149 (MQ=255)
                                                                         aTAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GCTTCTT‑A‑T‑TCTTTT‑GGAGTTTGGTTGCGCGGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgttt                                                           >  1:190389/1‑149 (MQ=255)
                                                                             aaaCAGCAATAATTTTATATTGCCTGAAAATATTTTTAATCTTTATTTAGAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgt                                                       <  1:87218/149‑1 (MQ=255)
                                                                              aaaaGCAATAATTTTATATTCACTGAAAATATTTTTAATCTCTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgtt                                                      <  2:413953/146‑1 (MQ=255)
                                                                                       aatTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGAtatt                                             <  2:440644/149‑1 (MQ=255)
                                                                                         ttttATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGAtattta                                           <  1:543929/149‑1 (MQ=255)
                                                                                          tttATATTAACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGATGCTATAAATTGAAT‑AATTAAAACTTTATTGCATTATCATAGGTTATTTGTCTTGTTTTTGTATGTGttatttat                                          >  2:390900/1‑149 (MQ=255)
                                                                                                       aaaaTATTTTTAATCTTTATTTATAGCTGTG‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAAGTTTTTTTCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATTCGTATGCtt                             >  2:566488/1‑149 (MQ=255)
                                                                                                       aaaaTAGTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGGTGCGATAAATTTAAT‑AATTTAAAGGTTGTTGCATATTTGTAGGATGTTTATGTGGTTTGTGTTTGTGATATTTATTATGCGTATGCtt                             <  2:386063/149‑1 (MQ=255)
                                                                                                           tatTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCaaa                         <  1:242151/149‑1 (MQ=255)
                                                                                                                 tattCTTTGTTTATAGGTGTT‑GTTTATT‑ATT‑T‑TTTT‑TGAGTTTGTTCGTGTTGCTATAAATTGAAA‑TATTAAAACTTCGTTGCATAGTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACaaa                   <  2:392286/146‑1 (MQ=255)
                                                                                                                 tAATTTGTATTTAGAGCTGTT‑GGTTATTGT‑T‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTTAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGCTATTTATTATGCGTATGCTTCAAAAACaaa                   <  2:392381/149‑1 (MQ=255)
                                                                                                                  aaTCATTATTTATAGCTGTT‑TGTTATT‑A‑T‑T‑TTTTGGGAGTTGGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATTTTCATAGGATATTTAGCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACaaaa                  <  2:597122/149‑1 (MQ=255)
                                                                                                                    tCTTTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAATCAAACtt                >  2:191827/1‑149 (MQ=255)
                                                                                                                    tCTTTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAAtt                >  1:428884/1‑149 (MQ=255)
                                                                                                                      tttatttaTAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTAt              <  1:517736/149‑1 (MQ=255)
                                                                                                                        tatttatGGCGGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCt            <  1:177306/149‑1 (MQ=255)
                                                                                                                               aGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACaaa                   >  2:644609/1‑135 (MQ=255)
                                                                                                                               aGCTGTT‑GGTTATT‑A‑TCT‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACaaa                   <  1:644609/135‑1 (MQ=255)
                                                                                                                                  tGTT‑GGTTATT‑ATT‑T‑TTTT‑GCAGTTTGGTTGCGGTGCTATCCTTTGAAT‑ATTTAAAACTTTACTACATCACAAATCGTGATTCACTCAATTGGTTTTCTTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGtt  <  1:318232/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                            
GTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑A‑T‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTT  >  NC_000913_3_hsa_tpiA/389778‑390055

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