Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_pgi4,296,0880CT36.2% 46.7 / 63.3 69intergenic (+266/+376)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base C (24/20);  new base T (24/1);  total (48/21)
Fisher's exact test for biased strand distribution p-value = 2.65e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCC  >  NC_000913_3_hsa_pgi/4295943‑4296235
                                                                                                                                                 |                                                                                                                                                   
ggATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGcc                                                                                                                                                  >  1:3163/1‑149 (MQ=17)
 gctacggCGTTTATGCCGCATCCGTCATGAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCt                                                                                                                                                 <  2:1030339/144‑1 (MQ=12)
 gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCAGCATCCGACATCAATGCCt                                                                                                                                                 >  2:24277/1‑149 (MQ=12)
  aTAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTg                                                                                                                                                >  1:1279108/1‑149 (MQ=32)
    aaGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                                                                                                              >  1:1625037/1‑149 (MQ=17)
    aaGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                                                                                                              >  2:383908/1‑149 (MQ=17)
    aaGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGGGTTATCAGGCCTACGCCAGACAGCGCAATAGCATGATTTAGCGTGATTTTGTAGTTAGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAt                                                                                                                                              >  2:942717/1‑149 (MQ=11)
     aGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCTACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATg                                                                                                                                             <  1:942717/149‑1 (MQ=17)
      ggCGTTGATGTCGGATGCGGCATAAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGc                                                                                                                                            >  1:94191/1‑149 (MQ=11)
       gCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTGTGTAGGTCGGATACGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCg                                                                                                                                           >  2:520253/1‑149 (MQ=255)
       gCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCg                                                                                                                                           >  1:646281/1‑149 (MQ=17)
       gCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCg                                                                                                                                           >  1:60309/1‑149 (MQ=17)
       gCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCATACTCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGGCGGATAAGTCGTTTACGACGCATCCGAAATAACCTCCTGATgtt                                                                                                                                           >  2:875933/1‑147 (MQ=11)
        ggTTTCTGCCGCCTCCTAGATCAACGCCTGATGCGACGCTTCACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGa                                                                                                                                          <  2:392825/148‑1 (MQ=11)
             aTGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCtt                                                                                                                                     >  1:569491/1‑149 (MQ=32)
              tGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTg                                                                                                                                    <  2:1279108/149‑1 (MQ=32)
                ccGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTGTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGcc                                                                                                                                  >  1:1332727/1‑149 (MQ=21)
                 cGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTCGCcg                                                                                                                                 >  1:1246007/1‑149 (MQ=18)
                  gCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCcgc                                                                                                                                >  1:1651998/1‑149 (MQ=25)
                  gCATCCGACAACAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTTATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCTTGATGCGACGCTTGCcgc                                                                                                                                <  2:1376653/149‑1 (MQ=255)
                    aTCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCTATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                                                                                                              <  1:1335014/149‑1 (MQ=25)
                     tCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTc                                                                                                                             >  1:1137048/1‑149 (MQ=25)
                     cacGACATCAACGCCTTATGCTACGCTTCCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCCTGATCTTGTAGGTCGGCTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTc                                                                                                                             <  2:1063819/147‑1 (MQ=11)
                        gACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAGTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCcgcgagtgc                                                                                                                          >  2:1388281/1‑144 (MQ=255)
                        gACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                                                          >  1:367945/1‑149 (MQ=25)
                        gACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                                                          >  2:513149/1‑149 (MQ=25)
                        gACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGACTACGCCAGACAGCGGAATAGCCTGATTAAGCGTGATATTGTAGGTCGGATAAGCCGTATACGCCGCATCCGCGATCAACCCCTGATGCGACGGTGGCCGCGTCGTa                                                                                                                          >  1:901286/1‑149 (MQ=11)
                           aTCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGAATTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCAGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCa                                                                                                                       >  1:192357/1‑149 (MQ=255)
                             cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                                                                                                     >  2:524346/1‑149 (MQ=25)
                              aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                    >  2:1507518/1‑149 (MQ=25)
                              aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                    >  1:829368/1‑149 (MQ=25)
                               aCGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                   <  2:587468/149‑1 (MQ=255)
                                 gCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACGGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTa                                                                                                                 >  1:1080136/1‑149 (MQ=25)
                                 gCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTa                                                                                                                 >  2:370335/1‑149 (MQ=25)
                                 gCCTGATAAGACGCGTTAAGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACTCTTGCCGCGTCTTATCAGGCCTa                                                                                                                 >  2:1143868/1‑149 (MQ=11)
                                  ccTGAAGCGACGCTTAACGCGTACTATCAGGCCTACGCCCGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGCTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                                                                                                                <  1:994611/149‑1 (MQ=11)
                                   cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                                                                                                               >  2:971406/1‑147 (MQ=18)
                                   cTGATAAGACGCGTTAAGCGTCGCATCAGGCGTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                                                                                                               >  1:1406674/1‑147 (MQ=11)
                                    tGATGCGACGCTTGACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGAc                                                                                                                                         <  1:760197/122‑1 (MQ=17)
                                    tGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatct                                                                                                              >  1:1553520/1‑146 (MQ=17)
                                     gataagACGCGTTAAGCGTCGCATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTAGGATAAGGCGGTTACGCCGCATCCGGAATCAACGCCTGATGCGACGCGTGCCGCGTCTTATAAGGCCTACGcg                                                                                                             >  1:825020/6‑148 (MQ=11)
                                      aTGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCa                                                                                                            <  1:640042/149‑1 (MQ=37)
                                        gCGACTCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGa                                                                                                          <  1:1171421/149‑1 (MQ=255)
                                          gACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAAGAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACa                                                                                                        >  1:768567/1‑149 (MQ=37)
                                                   cgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTg                                                                                                                                                <  1:259700/100‑1 (MQ=33)
                                                   cgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTg                                                                                                                                                >  2:259700/1‑100 (MQ=33)
                                                           aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAAGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGAGATAAAAGGCCTACTCCAGCCAGCGCAATAGCGGGAttt                                                                                       >  2:105030/1‑149 (MQ=11)
                                                            tCAGTCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGAGATTTTGTAGGTCGGCTAAGGCGTTTAGGCCGCTTCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa                                                                                      <  2:478159/149‑1 (MQ=11)
                                                            tCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa                                                                                      <  2:1332727/149‑1 (MQ=37)
                                                                   tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCt                                                                                                                                                 >  1:703012/1‑83 (MQ=33)
                                                                   tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCt                                                                                                                                                 <  2:703012/83‑1 (MQ=33)
                                                                   tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAt                                                                               >  1:1277460/1‑149 (MQ=37)
                                                                      gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAgc                                                                                                      >  2:220744/1‑123 (MQ=17)
                                                                         cgACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTAGGATAATGCTTTTACGCCGCATCAGACATCAACGCCTGATGCGACGCTTTCCGCGTCTTATCAGGCCTACGCCAGAAAGCGCAATAGCCTGATTTAGCGTGATTTTGTa                                                                         <  1:540628/148‑1 (MQ=11)
                                                                                  aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACCGCGCAGTACCCTGAAGTAGCGTGATTGTGAAGGTGGGATa                                                                >  2:322310/1‑149 (MQ=11)
                                                                                       ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGATTTTACGCCGAATCCGACATCAACGCCTGATGCGACGCTTTCCGCGTCTTATCAGGACTAAGAAAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGct                                                           >  1:999281/1‑148 (MQ=11)
                                                                                                  gTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCa                                                >  1:597329/1‑149 (MQ=37)
                                                                                                            atgTCGGATGCGGCGTAAACGCCTTATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCAACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCAAGCGACATCa                                      >  1:1526677/3‑149 (MQ=11)
                                                                                                              gTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAAc                                    >  2:617779/1‑149 (MQ=37)
                                                                                                               tCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCt                                                                                            <  2:109126/92‑1 (MQ=40)
                                                                                                                       ggCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGGCAGCCAGAGCAAGAGCCTGATGTAGCGTGAGTTGGGAGGTCGGATCAGGCGTGTATGCCGCCTCAGAGCGCAACGCCTGATGc                           >  2:1462877/1‑149 (MQ=11)
                                                                                                                              aCGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                                                                                                                >  2:1359939/1‑57 (MQ=255)
                                                                                                                                cccgcTTCCGACATCAACGCCATCTGCGACGCTCGCCGCGTCTTATGAGGAGTACGCCATCCAGCGCAATCGCCTGCTTTAGGGTGATTTTGTATGTCGGATAATGCGTTTATGCCGCACCCGACATCAACGCCTGATGCGACGCTTaa                  <  2:901286/148‑1 (MQ=11)
                                                                                                                                         gACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGAt                                                                                         <  2:936675/69‑1 (MQ=40)
                                                                                                                                         gACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGAt                                                                                         >  1:936675/1‑69 (MQ=40)
                                                                                                                                            aTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTc                                                                     >  1:47986/1‑86 (MQ=38)
                                                                                                                                            aTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTc                                                                     <  2:47986/86‑1 (MQ=38)
                                                                                                                                               aaCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGc   >  1:36411/1‑149 (MQ=32)
                                                                                                                                                aCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATACGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGcc  >  1:1127964/1‑149 (MQ=9)
                                                                                                                                                 |                                                                                                                                                   
GGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCC  >  NC_000913_3_hsa_pgi/4295943‑4296235

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