Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA4,296,0420CT26.7% 29.8 / 16.6 30intergenic (+266/+376)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base C (11/11);  new base T (3/5);  total (14/16)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.80e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAG  >  NC_000913_3_hsa_tpiA/4295895‑4296186
                                                                                                                                                   |                                                                                                                                                
tCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACg                                                                                                                                                 >  1:228549/1‑149 (MQ=32)
    aTAAGGCGTCTCTGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGCCCTACGCCAGCCAGCGCAATAGCCTGATGTATCGCCATTTTGTAGGCCGGATAAGTCGTGTACGCCGCATCCGACACCAACGCCTg                                                                                                                                             >  2:132988/1‑149 (MQ=11)
      aaGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAt                                                                                                                                           <  1:326980/149‑1 (MQ=32)
                   cGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGCCCCTTGCcg                                                                                                                              >  2:414115/1‑149 (MQ=255)
                      ttcagACATCAGCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCATGATTTAGCGTGATTTTGTAGGTCGGATACGGCGTTTCCGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGt                                                                                                                           <  2:321269/145‑1 (MQ=11)
                       cccggtgtCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGGCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTc                                                                                                                          <  1:416536/142‑1 (MQ=255)
                                aaCGCCGGATGCGACGCTTAACGCGTCTTATCAGGCTGACGCGAGACATGCCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                 <  2:297580/149‑1 (MQ=11)
                                   gCCTGATGCGACGCGTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTa                                                                                                              <  1:265448/149‑1 (MQ=21)
                                               gCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                 <  1:300183/134‑1 (MQ=21)
                                               gCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                 >  2:300183/1‑134 (MQ=21)
                                                     cgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGTCGCTTCTTATCAGGCCTACGCCAGACAGCGCAATAg                                                                                            <  1:414844/149‑1 (MQ=25)
                                                             aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGc                                                                                                                                    <  2:231681/101‑1 (MQ=18)
                                                             aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGc                                                                                                                                    >  1:231681/1‑101 (MQ=33)
                                                                 ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGCCAGCGCATTATCCTGCTTTAGCg                                                                                >  2:188381/1‑149 (MQ=255)
                                                                    cTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCCGACCTACGCCAGACCGCGCAAAAGCCTTATTTCGCGTGa                                                                             >  2:328964/1‑149 (MQ=11)
                                                                      aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGGTTTAGCGTgttg                                                                           >  2:286304/1‑146 (MQ=255)
                                                                            gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTTTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAg                                                                     <  1:302259/149‑1 (MQ=255)
                                                                                    aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAACCATCGCAATAGCCTGGTTTAGCGTGATTTTGTAGGTCGGATa                                                             >  2:322471/1‑149 (MQ=255)
                                                                                      tAGCCTGATTTAGGCTGCTTGTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCGGAGTTAGCGTGATTTTGTAGGTCGCATAAg                                                           <  1:407922/149‑1 (MQ=11)
                                                                                         ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                                                                                                  >  2:134793/1‑91 (MQ=16)
                                                                                         ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                                                                                                  <  1:134793/91‑1 (MQ=255)
                                                                                             aTTTTGCGTGATTTTGTAGGTGGGAGAAGGCGTCTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCGAGACGGCGCAATAGCCTGATTACGCGTGATTTTGTGGGTCGGGTAAGCCgataa                                                    >  2:319927/1‑145 (MQ=11)
                                                                                                           tGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGTTGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACa                                      <  2:228549/149‑1 (MQ=37)
                                                                                                             tAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATc                                    >  1:307250/1‑149 (MQ=39)
                                                                                                                                  gccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTCTCAGGCCTACGCCAGACCGCGCACTAGCCTGCTTTAGCGTGATTTCGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCCGATGCGACGCtttg               >  1:404980/1‑147 (MQ=11)
                                                                                                                                   ccgcgaccGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAAc              <  1:18240/143‑1 (MQ=21)
                                                                                                                                   ccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTg                            <  2:179818/135‑1 (MQ=35)
                                                                                                                                   ccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTg                            >  1:179818/1‑135 (MQ=37)
                                                                                                                                              aTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCa   <  1:286304/149‑1 (MQ=34)
                                                                                                                                               tCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAg  <  1:129297/149‑1 (MQ=33)
                                                                                                                                                   |                                                                                                                                                
TCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAG  >  NC_000913_3_hsa_tpiA/4295895‑4296186

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: