Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_pgi805,3120GT21.3% 125.5 / 23.7 75E604* (GAG→TAG) ybhJaconitase family protein
Reads supporting (aligned to +/- strand):  ref base G (28/31);  new base T (8/8);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTG  >  NC_000913_3_hsa_pgi/805165‑805456
                                                                                                                                                   |                                                                                                                                                
cTGATTCCTTCCGGTGAAACGTCTTCTTATAGTTCAAATACTATTGGTCTGGGGGAGTTTACCCTGTCTCGCAGCGATCCAGGTTATGATTGCAGAAGTAAAACGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGa                                                                                                                                                 >  2:598187/1‑149 (MQ=255)
 tGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAg                                                                                                                                                >  1:484047/1‑149 (MQ=255)
    ttccttccGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTg                                                                                                                                             >  1:803958/1‑149 (MQ=255)
       cttccGGTGATACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGAGGTGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACa                                                                                                                                          <  2:1413786/149‑1 (MQ=255)
        ttccGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCTCGCTGACat                                                                                                                                         <  1:1292316/149‑2 (MQ=255)
        ttccGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACag                                                                                                                                         >  2:1087157/1‑149 (MQ=255)
        ttccGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACag                                                                                                                                         >  2:347311/1‑149 (MQ=255)
         tccGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATAAGCGTCTGGCGGGGAATGTCAGCGAGctgactga                                                                                                                                        >  2:982916/1‑149 (MQ=255)
         tccGGTGAAACCTCTTCTTATCGTTCAAAACCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCTCGCTGACatt                                                                                                                                        <  2:151450/149‑3 (MQ=255)
         tccGGCGAAACCTCTGCTTTTCCTTCAAATCCCATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACaga                                                                                                                                        <  2:1346378/149‑1 (MQ=255)
          ccGGTGAACCCTCTTCTTATCGTGCAAATCAGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCGGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCTCGCTGACattc                                                                                                                                       <  2:1058103/149‑4 (MQ=255)
          ccGGTGAAACCTCTTCTTATCGTTCAACTCCGATTGGTCTTGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACagag                                                                                                                                       <  1:833315/149‑1 (MQ=255)
          ccGGTGAAACCTCTTCTTATCGGTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCCCCGCGATCCCGGATATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCTCGCTGACattc                                                                                                                                       <  2:337036/149‑4 (MQ=255)
             gTGAAACCTCTTCTTCTCGTTCAAAGCCGCTCGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGCCTGCTTAGCGGGAAAATCAGCGTCTGGCGGGGAATGTCAGCTCGCTGACattcccc                                                                                                                                    <  2:1110677/149‑7 (MQ=255)
             gTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCTCGCTGACattcccc                                                                                                                                    <  2:1488690/149‑7 (MQ=255)
               taaaCCTCTTCTTATCGTTCAAATCCTATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATTTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGtt                                                                                                                                  <  1:1340135/148‑1 (MQ=255)
                       ttcttcTCGTGCAACTCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCa                                                                                                                          <  2:148754/149‑1 (MQ=255)
                        tcttATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCAt                                                                                                                         <  2:865597/149‑1 (MQ=255)
                          ttATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTa                                                                                                                       >  2:394589/1‑149 (MQ=255)
                                ttCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGa                                                                                                                 >  2:593119/1‑149 (MQ=255)
                                   aaaTCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTg                                                                                                              >  1:253015/1‑149 (MQ=255)
                                     aTCCGATTGGGCTGGCGGAGGTTACCCTGTCTCGCCGCGATCCCGGTTATGTTCGCAGAAGTACCGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGGGTTTGCGCGCATTAAGCAGATTGCt                                                                                                            <  1:1453157/149‑1 (MQ=255)
                                      tCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGGCTGGCGGGGAATGTCAGCGAGGTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTg                                                                                                           >  2:1320084/1‑149 (MQ=255)
                                        cGATTGGTCTGGCGGAGTTTACCCTGTGTCGCCGAGATCCCGTTTATGTTCGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTCGCGCGCATTAAGCAGATTGCTGGt                                                                                                         <  1:1050747/149‑1 (MQ=255)
                                             ggTCGGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGGCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGa                                                                                                    <  2:584855/149‑1 (MQ=255)
                                                 tGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCat                                                                                                <  1:1465258/149‑1 (MQ=255)
                                                   gCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAATCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCatat                                                                                              >  1:860523/1‑149 (MQ=255)
                                                       aGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGTCTGCTGAGCTGGAAAATCTGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGAt                                                                                          <  2:377534/149‑1 (MQ=255)
                                                        gTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATc                                                                                         <  1:472874/149‑1 (MQ=255)
                                                         tatACCCTTTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATcc                                                                                        <  2:1137768/147‑1 (MQ=255)
                                                              ccTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTTAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGc                                                                                   >  2:125739/1‑149 (MQ=255)
                                                              ccTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGCTTGCTGGTCAGGAGCATATTGATCCGCTGc                                                                                   >  1:388979/1‑149 (MQ=255)
                                                               cTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCa                                                                                  >  1:1321362/1‑149 (MQ=255)
                                                                tGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCaa                                                                                 >  1:651544/1‑149 (MQ=255)
                                                                tGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCaa                                                                                 >  1:1429383/1‑149 (MQ=255)
                                                                   ctcGCGGCGTGCCCGCGTCTGTTTGCAGAAGTTTGTGGACTTCTGAGCTGGAAATTCAGTGTCTGGTTGGGAATGTCAGCTTGCTTTAAGATGTGTTTGCTCGCAATAAGAATATTGCTTTTCAGGAGCATATTGATCCGCTGCAAACt                                                                              <  2:1061077/149‑1 (MQ=255)
                                                                    tcGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTg                                                                             >  1:1016484/1‑149 (MQ=255)
                                                                      gccgcGATCCCGGTTATGTTAGCAGAAGTAAAGCGACGGCTGAGCTGGAAACTCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGaa                                                                           <  1:982916/149‑1 (MQ=255)
                                                                         gcgATCCCGGTTATGTTAGCAGCAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTAAGCTAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAAtt                                                                        <  1:684972/149‑1 (MQ=255)
                                                                         gcgATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAAtt                                                                        <  2:1321362/149‑1 (MQ=255)
                                                                         gcgATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAAtt                                                                        >  2:136588/1‑149 (MQ=255)
                                                                          cgATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTg                                                                       >  1:882232/1‑149 (MQ=255)
                                                                              cccGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTATGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgcag                                                                   <  1:1439767/149‑1 (MQ=255)
                                                                                 gggtATGTTAGCAGAAGTAACGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCAt                                                                <  2:1151443/146‑1 (MQ=255)
                                                                                    tATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGt                                                             <  1:593119/149‑1 (MQ=255)
                                                                                      tGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCt                                                           <  1:136588/149‑1 (MQ=255)
                                                                                         tttcAGACGTAAAGCGACTGCGGAGCTTGTACACCCTCGGCTGACGTGCTATGTCAGCGAGCTGACAGAGGTCTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATg                                                        <  1:669562/146‑1 (MQ=255)
                                                                                         tAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATg                                                        <  2:253015/149‑1 (MQ=255)
                                                                                           gCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGGCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGAATATTGATCCGCTGCTAACTGAAATTGGAAGCATGTTCTATGCg                                                      >  2:1123412/1‑149 (MQ=255)
                                                                                             aGAAGTAACGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGt                                                    <  1:1274800/149‑1 (MQ=255)
                                                                                              gAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTg                                                   >  1:484930/1‑149 (MQ=255)
                                                                                              gAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCGGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTg                                                   <  2:484047/149‑1 (MQ=255)
                                                                                                aGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCTCGCATTAAGCAGATTGCTGGTCAGGAGCATTTTGATCCGCTGCATACTGAAATTGGCAGCATGGTCTATGCGGTGaa                                                 <  1:1155539/149‑1 (MQ=255)
                                                                                                 gTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGaaa                                                <  1:1087157/149‑1 (MQ=255)
                                                                                                 gTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGaaa                                                >  2:367939/1‑149 (MQ=255)
                                                                                                 gTAAAGCGACTGCTGAGCAGGAAAATCAGCGTCTGGCGGGGAATCTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGAATGCTGGTCAGGAGCATATTTATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGaaa                                                >  2:768389/1‑149 (MQ=255)
                                                                                                             ctgAGCAGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGt                                    >  1:74742/1‑149 (MQ=255)
                                                                                                                 gCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTg                                >  1:426753/1‑149 (MQ=255)
                                                                                                                  cTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTgc                               <  2:803958/149‑1 (MQ=255)
                                                                                                                    ggAAAATCAGCGTCTGGCGGGGAATGTCAGCGATCTGACAGAGGTGTTTGCGCGCATTAAGCAGTTTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTTTGCGGTGAAACCAGGCGATGGTTCTgcgc                             <  1:347311/149‑1 (MQ=255)
                                                                                                                           cAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGAATGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTGTATGCGGTTAAACAATGCGATGGTTCTGCGCGTGTACa                      >  2:994078/1‑149 (MQ=255)
                                                                                                                             gCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAgg                    >  1:409697/1‑149 (MQ=255)
                                                                                                                                 cTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAggcggc                <  2:1468395/149‑1 (MQ=255)
                                                                                                                                  tGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGgcggcg               >  1:1247991/1‑149 (MQ=255)
                                                                                                                                       acctctGTCAGCTCGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTg          >  1:775040/6‑149 (MQ=255)
                                                                                                                                        ggAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGc         <  2:1016484/149‑1 (MQ=255)
                                                                                                                                        cctctGTCAGCTCGCTGACAGAGGTGTTTGCGAGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGCTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGc         >  2:1173506/5‑149 (MQ=255)
                                                                                                                                        cctctGTCAGCTAGCTGAAATATGTGTATGCGCGCATTAAGCAGATTTCTTGACAGTAGCATATTTAACCGCTGCAAACTTAAATTGGAAGAATGGTATATGCGGTGAAACCAGGCTATGGTTATGCGCGTGAACCGGCGCCTAGCTGc         >  2:733945/5‑149 (MQ=255)
                                                                                                                                         ctctGTCAGCTCGCTTACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGCTTCTGCGCGTGAACAGGCGGCGAGCTGcc        >  2:1475700/4‑149 (MQ=255)
                                                                                                                                         ctctGTCAGCTCGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAATCTGAAATTTGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGcc        >  2:94065/4‑149 (MQ=255)
                                                                                                                                         ctctGTCAGCTCGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCAGGGTCGATGCGGTGAAGCCAGGCGATGGTTCTGCGCGTGCACCGGCGGCGAGCTGcc        >  2:724304/4‑149 (MQ=255)
                                                                                                                                         ctcagTCAGCTCGCTGCCAGAGGTGTTTGCGCGCATTAAGCAGATTGCACGTCAGGAGCATATTCATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGAGAATCCAGGCGATCGTTCTGCGCGTGAACAGGCGTCGAGCTGcc        >  1:276362/5‑149 (MQ=255)
                                                                                                                                          tctGTCAGCTCGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCa       >  1:780970/3‑149 (MQ=255)
                                                                                                                                               cAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCgtg  <  2:651544/149‑1 (MQ=255)
                                                                                                                                               cAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCGGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTGCTGCGCGTGAACAGGCGGCGAGCTGCCAGCgtg  <  2:1429383/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTG  >  NC_000913_3_hsa_pgi/805165‑805456

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: ‑