Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA4,296,0420CT25.0% 78.6 / 26.8 57intergenic (+266/+376)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base C (19/23);  new base T (9/5);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC  >  NC_000913_3_hsa_tpiA/4295898‑4296176
                                                                                                                                                 |                                                                                                                                       
gATAAGGCGTTTAGGCCGCAGCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTTCGTGTGGCGTTTTACGCCGCATCCGACATCAATGCCt                                                                                                                                     <  2:655305/149‑1 (MQ=11)
gATAAGGCGTTTACGCCGGATACGGCATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCt                                                                                                                                     <  2:390037/149‑1 (MQ=255)
   aaGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                                                                                                  <  2:185912/149‑1 (MQ=17)
     ggCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGc                                                                                                                                >  1:124655/1‑149 (MQ=32)
       cttaaaTGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGa                                                                                                                              >  2:592747/6‑149 (MQ=12)
            atcccccATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCtt                                                                                                                         <  1:564796/143‑1 (MQ=17)
             gactgCATCCGACATGAATGCGTGATGCGCCTCTTAACGCGTCTTATCAGGCACACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGGCGTACAAGGGGTTTCCGCCGCATACTACATCAATGCCTGATGCGACGCTTg                                                                                                                        <  2:526948/145‑1 (MQ=11)
              gCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGc                                                                                                                       >  1:247790/1‑149 (MQ=21)
                cGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGAc                                                                                                                             <  2:382994/141‑1 (MQ=17)
                cGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGAc                                                                                                                             >  1:382994/1‑141 (MQ=32)
                   aTCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGGAATAGCCCGATCTAGCGTGATTTTAGAGGTCGGATAAGGCGGTTACGCAGCATACGACTCCAATGCCTTATGCGACGCttcgcgcct                                                                                                                  >  2:60670/1‑142 (MQ=11)
                   aTCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                                                                                                  >  1:284560/1‑149 (MQ=25)
                    tCCGACATCAACGCCTGATGTGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTAATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTc                                                                                                                 <  2:534912/149‑1 (MQ=255)
                    tCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCATGATTTAGCGTGATTTTGTAGGACGGATAAGGCGTTTACGCCGCGTCCGACATCAATGCCGGATGCGGCGCTTGCCGCGTc                                                                                                                 >  1:33478/1‑149 (MQ=11)
                    gtcGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGATTACGCCAGACACCGCGATAGCCTGACTTAGCGTAATTTTGTTGGTCGGATAAGGCCTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTc                                                                                                                 <  2:261296/147‑1 (MQ=11)
                       gACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAc                                                                                                                                         >  1:261879/1‑122 (MQ=32)
                       gACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCTGCGTCTTa                                                                                                              >  2:438283/1‑149 (MQ=21)
                        aCATCAACGCCTGTTGCTACGCTTATCGCGTCTTATCAGGCCTATGCCAGATAGAGCAATAGCCTGCTTTAGCGTGATATAGTAGGTACAATACTGCGTCTCTGCCCCACCCGATATCGACGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                            <  2:291323/149‑1 (MQ=2)
                          aTCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCa                                                                                                          >  1:235451/1‑149 (MQ=21)
                          aTCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCa                                                                                                          >  2:354044/1‑149 (MQ=38)
                             aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAt                                                                                                                                  >  1:33254/1‑123 (MQ=32)
                              aCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                      >  2:156822/1‑149 (MQ=37)
                                  cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                                                                                                  >  2:550245/1‑147 (MQ=17)
                                         gACGCTTAACCCGTCTCATCAGGCCTACGCGAGAGAGCGCGTTAGCCTGATGTAGCGGGATTTTGTAGGTCGGATAAGGCGTTTACGCCGTATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACa                                                                                           <  2:247961/149‑1 (MQ=11)
                                             cTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAgcgc                                                                                       <  1:412154/149‑1 (MQ=39)
                                                  cgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                                                                                                  <  2:41416/118‑1 (MQ=37)
                                                  cgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                                                                                                  >  1:41416/1‑118 (MQ=255)
                                                   gcgTCTTATCAGGCCTACGCCAGAGGGCGCAATGGCCTGATTTAGGGTGATTTTGTAGGGCGGATAAGGCGTTTACGCCGCACCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc                                                                                 <  2:81509/149‑1 (MQ=11)
                                                          aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGAGATTTTGTAAGTCCGATAAGGCGCTTACGCGGCATCCGACATCAAAGCCTGATGCGACACTTGCCGCGTCTTATCAGGCCTACGTCAATCATCGCAAAAGCCTGAttt                                                                          >  2:156048/1‑149 (MQ=11)
                                                           tCAGGCCTACGCCAGACAGCGCGATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa                                                                         <  1:429348/149‑1 (MQ=255)
                                                             aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTACACAAGCGTTTTGAGACAGCGCAATAGCCTGATTTAGc                                                                       <  2:488279/149‑1 (MQ=11)
                                                              ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACTCTTGCCGCGTCCTATCAGGCCAACGCCAGACAGCGCCATAGCCTGCTTGAGCg                                                                      >  2:670663/1‑149 (MQ=11)
                                                              ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCg                                                                      >  1:360735/1‑149 (MQ=39)
                                                                      ccAGACCGCGAAATGGCCTGATTCAGAGTTAATTTGTGGGGCGGGTCAGGGGTTTACGCCGCGTCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTg                                                              <  2:217795/149‑1 (MQ=11)
                                                                       cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGt                                                             <  2:360735/149‑1 (MQ=38)
                                                                            aGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTc                                                        >  2:360321/1‑149 (MQ=38)
                                                                                  aTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATaa                                                  <  1:308267/149‑1 (MQ=38)
                                                                                   tAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTATGTCGGATaaa                                                 >  2:563858/1‑148 (MQ=38)
                                                                                    aGCCTGATGTAGCGTGATTTTGTAGGTTGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGGTGTGGTAGGCCGGATAAgg                                                <  2:490528/149‑1 (MQ=11)
                                                                                      ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                           >  2:261472/1‑88 (MQ=18)
                                                                                      ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                           <  1:261472/88‑1 (MQ=255)
                                                                                       cTGATTTCGCGTGAATCGGTGGGGCGGACGAGGCGTTTACGCCGCATCCGGCATGAACGCCTGATGCGACGCTTGCCGCGTCGTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGt                                             <  1:372712/149‑1 (MQ=11)
                                                                                          aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTc                                                        >  2:443995/1‑135 (MQ=17)
                                                                                          aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTc                                                        <  1:443995/135‑1 (MQ=40)
                                                                                                      tttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                       <  1:555329/76‑1 (MQ=255)
                                                                                                      tttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                       >  2:555329/1‑76 (MQ=255)
                                                                                                         gTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGCCGGATCGGGCGATAATCCCGCATCCAACAt                           >  1:154176/1‑149 (MQ=11)
                                                                                                           aGGTCGGATAAGGCGTTTACGCCGCATCCGACATAAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCCATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATCAGGCGTTTACGCCGCATCCGACATcc                         >  1:63016/1‑148 (MQ=255)
                                                                                                            ggTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAACGCGTTTATGCCGCATCCCACCTCaa                        >  2:187076/1‑149 (MQ=255)
                                                                                                             gTCGGATAAGGCGTTTACGGCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTCGCGTGATTTTGTAGGTCGGATAAGGCGTTTCTGCCCCATCCGACATCCAc                       >  2:572014/1‑149 (MQ=11)
                                                                                                                   tAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGa                 >  1:185912/1‑149 (MQ=38)
                                                                                                                     agccGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTAGGTGCAATCTGGTATGCCGCATCCGACATCAACGCCTGATg               <  1:670663/146‑1 (MQ=11)
                                                                                                                      ggCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTc                                                                                                                 >  2:230827/1‑51 (MQ=255)
                                                                                                                      ggCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTc                                                                                                                 <  1:230827/51‑1 (MQ=255)
                                                                                                                                gccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAt                                                                  >  2:279374/1‑88 (MQ=40)
                                                                                                                                gccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAt                                                                  <  1:279374/88‑1 (MQ=40)
                                                                                                                                  cgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCGGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGc                                               <  1:606486/105‑1 (MQ=38)
                                                                                                                                   gcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAAcgc  <  1:156822/149‑1 (MQ=35)
                                                                                                                                                 |                                                                                                                                       
GATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC  >  NC_000913_3_hsa_tpiA/4295898‑4296176

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