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.

GTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTT  >  NC_000913_3_hsa_tpiA/389778‑390055
                                                                                                                                                |                                                                                                                                            
gTAGGGATGTTTTATGAAAGTTGTCGCGATGTTTACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAAGAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAttttttt                                                                                                                                     <  2:361346/149‑2 (MQ=255)
         ttttACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTgg                                                                                                                            >  2:242151/1‑149 (MQ=255)
         ttttACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGCTATTATTTTTTTGGGGTTTgt                                                                                                                            >  2:10265/1‑148 (MQ=255)
              cGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTgcg                                                                                                                       >  1:103646/1‑149 (MQ=255)
                aaaGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGTTTGCgct                                                                                                                     >  2:13011/1‑149 (MQ=255)
                                 gACAAGAAGAGAATGGAAGAGAGGCTAGGACAGCTTTTTAATATAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGaatcac                                                                                                    <  2:173356/149‑4 (MQ=255)
                                     agaagaGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTaaa                                                                                               >  2:517736/1‑149 (MQ=255)
                                         gagaATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACtt                                                                                           >  1:63292/1‑149 (MQ=255)
                                            aaTGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTttgt                                                                                        >  2:236034/1‑149 (MQ=255)
                                                gAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTCGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCa                                                                                    <  1:581908/149‑1 (MQ=255)
                                                         ccAGGACATATATTTAATAGAAAAATCAATAATTTTAGATTCACTGAAAATATTTTTAATCTTCATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCaccc                                                                           <  2:611149/149‑4 (MQ=255)
                                                           aGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGa                                                                         >  1:440644/1‑149 (MQ=255)
                                                           aGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGa                                                                         >  1:97803/1‑149 (MQ=255)
                                                                      ttAATAGTAACAGCAATAATTATATATGCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTTGTTATTATTTTTTTGGAGTTTGCTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTgg                                                              <  2:309983/149‑1 (MQ=255)
                                                                       tAATAGAGACATCAATAATTTTATATTCACTGAAATTATTTTTAATCTTTATTTATAGCTGTTTGTTATTATTTTTTTGGAGTTTGGTTGCACTGCTATAAATTTAATCAATCACAAAGATCGGTATATTCATAGGATATTTATCTGgt                                                             <  2:305717/149‑1 (MQ=255)
                                                                        aaTAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTGTAGCTGTGGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAGACTTTGTTGCATATTCATAGGATATTTCTGTGgtt                                                            >  1:22429/1‑149 (MQ=255)
                                                                         aTAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTAGATCTTTTTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgttt                                                           >  1:125384/1‑149 (MQ=255)
                                                                         aTAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgttt                                                           >  1:434651/1‑149 (MQ=255)
                                                                         aTAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGCTTCTTATTCTTTTGGAGTTTGGTTGCGCGGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgttt                                                           >  1:190389/1‑149 (MQ=255)
                                                                             aaaCAGCAATAATTTTATATTGCCTGAAAATATTTTTAATCTTTATTTAGAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgt                                                       <  1:87218/149‑1 (MQ=255)
                                                                              aaaaGCAATAATTTTATATTCACTGAAAATATTTTTAATCTCTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgtt                                                      <  2:413953/146‑1 (MQ=255)
                                                                                       aatTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGAtatt                                             <  2:440644/149‑1 (MQ=255)
                                                                                         ttttATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGAtattta                                           <  1:543929/149‑1 (MQ=255)
                                                                                          tttATATTAACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGATGCTATAAATTGAATAATTAAAACTTTATTGCATTATCATAGGTTATTTGTCTTGTTTTTGTATGTGttatttat                                          >  2:390900/1‑149 (MQ=255)
                                                                                                       aaaaTATTTTTAATCTTTATTTATAGCTGTGGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAAGTTTTTTTCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATTCGTATGCtt                             >  2:566488/1‑149 (MQ=255)
                                                                                                       aaaaTAGTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGGTGCGATAAATTTAATAATTTAAAGGTTGTTGCATATTTGTAGGATGTTTATGTGGTTTGTGTTTGTGATATTTATTATGCGTATGCtt                             <  2:386063/149‑1 (MQ=255)
                                                                                                           tatTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCaaa                         <  1:242151/149‑1 (MQ=255)
                                                                                                                 tattCTTTGTTTATAGGTGTTGTTTATTATTTTTTTTGAGTTTGTTCGTGTTGCTATAAATTGAAATATTAAAACTTCGTTGCATAGTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACaaa                   <  2:392286/146‑1 (MQ=255)
                                                                                                                 tAATTTGTATTTAGAGCTGTTGGTTATTGTTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTTAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGCTATTTATTATGCGTATGCTTCAAAAACaaa                   <  2:392381/149‑1 (MQ=255)
                                                                                                                  aaTCATTATTTATAGCTGTTTGTTATTATTTTTTGGGAGTTGGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATTTTCATAGGATATTTAGCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACaaaa                  <  2:597122/149‑1 (MQ=255)
                                                                                                                    tCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAATCAAACtt                >  2:191827/1‑149 (MQ=255)
                                                                                                                    tCTTTATTTATAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAAtt                >  1:428884/1‑149 (MQ=255)
                                                                                                                      tttatttaTAGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTAt              <  1:517736/149‑1 (MQ=255)
                                                                                                                        tatttatGGCGGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCt            <  1:177306/149‑1 (MQ=255)
                                                                                                                               aGCTGTTGGTTATTATTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACaaa                   >  2:644609/1‑135 (MQ=255)
                                                                                                                               aGCTGTTGGTTATTATCTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACaaa                   <  1:644609/135‑1 (MQ=255)
                                                                                                                                  tGTTGGTTATTATTTTTTTGCAGTTTGGTTGCGGTGCTATCCTTTGAATATTTAAAACTTTACTACATCACAAATCGTGATTCACTCAATTGGTTTTCTTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGtt  <  1:318232/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                            
GTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTT  >  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: