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

CGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCA  >  NC_000913_3_hsa_tpiA/4295896‑4296185
                                                                                                                                                  |                                                                                                                                               
cgggcaaGGCATTTCTACCCCATCATACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGc                                                                                                                                               <  2:625007/144‑1 (MQ=11)
cGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGc                                                                                                                                               <  1:429544/149‑1 (MQ=17)
cGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGc                                                                                                                                               <  2:191379/149‑1 (MQ=32)
  gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTCACGCTTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCt                                                                                                                                             >  2:3730/1‑149 (MQ=255)
  gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGCCATCAACGCCt                                                                                                                                             >  1:336400/1‑149 (MQ=17)
  gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCt                                                                                                                                             <  1:35273/149‑1 (MQ=17)
  gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGcgt                                                                                                                                             >  1:205531/1‑147 (MQ=17)
       ggCGTTTATGCCGCACCCTACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGc                                                                                                                                        <  1:152814/149‑1 (MQ=17)
        gCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATGCGACATCAATGCCTGATGCg                                                                                                                                       >  1:44771/1‑149 (MQ=17)
        gCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCg                                                                                                                                       >  1:448242/1‑149 (MQ=17)
         gatttACGCCGCATCCGACATCAACGCCTGATGCGACGATTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGa                                                                                                                                      <  2:747535/147‑1 (MQ=12)
          gTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGAc                                                                                                                                     <  1:80125/149‑1 (MQ=32)
          gTTTATGCCGCATCCGACAGCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGAc                                                                                                                                     <  1:183719/149‑1 (MQ=17)
                gCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTAGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGc                                                                                                                               >  2:719971/1‑149 (MQ=18)
                 ccGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGcc                                                                                                                              >  2:130507/1‑149 (MQ=34)
                 ccGCATCAGGCACTAACGCCTGATGCGACGCTTAACGCGTCTTATCAGTGCTACGCTAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGcc                                                                                                                              <  2:205094/149‑1 (MQ=11)
                  cGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATCGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACCGTTGCcg                                                                                                                             >  2:360114/1‑149 (MQ=16)
                  cGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCcg                                                                                                                             >  1:609996/1‑149 (MQ=35)
                   gCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCACGCCCACGCCAGACAGCGCAGTCGCTTTATTTAGCGTTGTTTTTTAGGTCTGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCcgc                                                                                                                            <  1:558086/149‑1 (MQ=11)
                       ccGACATGAAAGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCt                                                                                                                        <  2:419807/149‑1 (MQ=255)
                         gACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                                                      >  1:756752/1‑149 (MQ=25)
                         gACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                                                      >  1:474964/1‑149 (MQ=25)
                          aCATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                                     >  1:734684/1‑149 (MQ=21)
                             tCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                                                  >  1:158276/1‑149 (MQ=39)
                             tCAACGCCTGATGCGACCCTTAGCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                                                  <  1:90282/149‑1 (MQ=255)
                              cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                                                                                                 <  1:377685/149‑1 (MQ=25)
                               caCGGCTGATGGGACGCTTGACGTGTCTTATCAGGCCTAAGTCCGGCAGGGAAATGGGCTGATTTAGGGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                <  2:121341/148‑1 (MQ=11)
                                aCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                               >  2:132382/1‑149 (MQ=37)
                                 cGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCt                                                                                                              <  2:609996/149‑1 (MQ=39)
                                  gCCTGATCCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTa                                                                                                             <  1:34842/149‑1 (MQ=255)
                                  gCCTGATAAGACGCGTTAAGCGTCGCATCAGGCCTCCGCCAGCCAGCGCCATAGCCTGATTCAGCCTATTTCTGTAGGGCGGCTGAGGCGTTTACGCCGCATCCGACATCAATGCCTGTGGCGACCCTTGCCGCGTCTTATCAGGCCtc                                                                                                             >  2:747552/1‑148 (MQ=2)
                                   ccTTATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                                                                                                            <  1:387807/149‑1 (MQ=255)
                                   ccTGATGGGAGGCTTACGGCGTCTTAGCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                                                                                                            <  2:336400/149‑2 (MQ=11)
                                    cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                                                                                                           >  1:321700/1‑147 (MQ=18)
                                    cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                                                                                                           <  2:734684/149‑3 (MQ=17)
                                    cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                                                                                                           <  1:310361/149‑3 (MQ=17)
                                       aTGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCa                                                                                                        <  1:288019/149‑1 (MQ=255)
                                       aTGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                                                                                                           >  2:466824/1‑144 (MQ=18)
                                       aTGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                                                                                                           <  1:466824/146‑3 (MQ=17)
                                           gACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACa                                                                                                    >  1:733050/1‑149 (MQ=38)
                                            aCGCTTACCGCGTCTTATCAGGTCTACGCCATACAGAGCAATAGCATGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAg                                                                                                   <  2:274621/149‑1 (MQ=11)
                                              gCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGAcgtcg                                                                                                 >  1:349000/1‑145 (MQ=37)
                                                ttAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGCCGCTTGCCGCGTCTTATCAGGCCTACGCCAGACCGCGCa                                                                                               >  1:296035/1‑149 (MQ=255)
                                                 tgccgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCaa                                                                                              <  1:422851/146‑1 (MQ=37)
                                                  aaCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                                                                                                 <  2:267485/129‑1 (MQ=25)
                                                  aaCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                                                                                                 >  1:267485/1‑129 (MQ=40)
                                                   aCGCGTCTTATCAGGCCTACGACAGACGGCGCAATGGCATGATTTAGCGTGATTTTGTAGGTCGGACAAGGCGTTTACGCCGCATCCGGCATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATa                                                                                            <  2:9759/149‑1 (MQ=11)
                                                             tCAGGCCTACGTCCGACAGCGCAATGACCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCGTCCGACATCAACGCCTGATGGTACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa                                                                                  <  1:525864/149‑1 (MQ=11)
                                                             tCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCATTAGCCTGCTTTa                                                                                  >  1:40531/1‑149 (MQ=255)
                                                             tCAGGCCTACGCCAGACAGCGCAAGAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCGGCGTCTTATGAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa                                                                                  >  1:382/1‑149 (MQ=255)
                                                               aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGc                                                                                <  1:132382/149‑1 (MQ=39)
                                                                ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATATCCTGAGTTAGCg                                                                               >  1:46589/1‑149 (MQ=255)
                                                                   cTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATa                                                                                            >  1:62766/1‑133 (MQ=40)
                                                                   cTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATa                                                                                            <  2:62766/133‑1 (MQ=18)
                                                                    tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa                                                                                  >  2:129417/1‑142 (MQ=38)
                                                                    tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa                                                                                  <  1:129417/142‑1 (MQ=39)
                                                                       gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAtttt                                                                        >  2:672379/1‑149 (MQ=38)
                                                                          aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTa                                                                     <  2:733050/149‑1 (MQ=37)
                                                                              aGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACCCTAGACATCGCAAAAGCTTTATTTATCGTGATTTTCTAGATc                                                                 >  2:177703/1‑149 (MQ=14)
                                                                                 aaaaTAGACTGATTTAGCGTCAATTTGTGGGGCGTATAAGGCGTTTACGCCGCATCAGACATAAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGa                                                              <  2:702164/147‑1 (MQ=11)
                                                                                        ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTa                                                                     >  1:582249/1‑135 (MQ=38)
                                                                                        ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTa                                                                     <  2:582249/135‑1 (MQ=39)
                                                                                        ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATCAGGCg                                                       >  1:272619/1‑149 (MQ=37)
                                                                                         cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGCCGGCGCAATAGCCTGCTTTACCTTAACTCTCGAGATCGGATAATGCGt                                                      >  2:130559/1‑149 (MQ=11)
                                                                                            aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCCCAGCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTa                                                   <  1:207771/149‑1 (MQ=255)
                                                                                                       ttttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCTGCATGCg                                        >  2:682725/1‑149 (MQ=255)
                                                                                                               gTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAAc                                <  2:205531/149‑1 (MQ=37)
                                                                                                                tCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCa                                                                                                                   <  1:568299/65‑1 (MQ=255)
                                                                                                                tCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCa                                                                                                                   >  2:568299/1‑65 (MQ=255)
                                                                                                                 gggATAGGGCGTGTCCGCCGGACGCGAGATAAACGCCTGATGCGACGCTTGCCGCGTCGTATCAGGCCTGCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGc                              <  1:321389/148‑1 (MQ=11)
                                                                                                                    aTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACCCCAGCCAGCGCAGTAGCCTGATTTAGCGGGATTTTGTAGTTCGGATAAGGCGGTTATGCCGCATCCGACATCAATGCCTg                           >  2:11419/1‑149 (MQ=14)
                                                                                                                        ggCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                                                      <  1:657898/54‑1 (MQ=255)
                                                                                                                        ggCGTATACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                                                      >  2:657898/1‑54 (MQ=11)
                                                                                                                         gCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACCCCTGTTGCg                      >  1:516813/1‑149 (MQ=255)
                                                                                                                           ggttACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGGGATTTTGTAGGTCGGATAAGGGGTTTATGCCGCATCCGACCTCAACGCCTGATGCGAc                    >  1:35412/3‑149 (MQ=255)
                                                                                                                             atACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGTGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGc                  <  1:688173/148‑1 (MQ=255)
                                                                                                                                 gccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTaa              >  1:416984/1‑149 (MQ=35)
                                                                                                                                 gccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTaa              >  1:51449/1‑149 (MQ=34)
                                                                                                                                         cGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGTGATTTTGTAGGTCGGATAATGCGTTTATTCCGCATCCGACCTTAACGCCTCATGCGGTGCTTGCCGCGTCtt      >  2:92433/1‑149 (MQ=11)
                                                                                                                                         cGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATc                                   <  1:371656/120‑1 (MQ=17)
                                                                                                                                         cGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATc                                   >  2:371656/1‑120 (MQ=18)
                                                                                                                                          gACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTa     >  2:206854/1‑149 (MQ=35)
                                                                                                                                             aTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCa  >  1:357437/1‑149 (MQ=34)
                                                                                                                                                 aCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa                                                                                  >  2:643270/1‑65 (MQ=33)
                                                                                                                                                 aCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa                                                                                  >  2:641990/1‑65 (MQ=33)
                                                                                                                                                 aCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa                                                                                  <  1:641990/65‑1 (MQ=33)
                                                                                                                                                 aCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa                                                                                  <  1:643270/65‑1 (MQ=33)
                                                                                                                                                  |                                                                                                                                               
CGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCA  >  NC_000913_3_hsa_tpiA/4295896‑4296185

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: