New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 1253913 =NA (NA)36 (0.180) 6/164 23.3 0.34 intergenic (+71/+20) RCAP_rcc01181/hisH hypothetical protein/imidazole glycerol phosphate synthase subunit HisH
?NC_014034 = 2590221 141 (0.580)intergenic (‑120/‑136) glnE/RCAP_rcc02415 glutamate‑ammonia‑ligase adenylyltransferase/hypothetical protein
Rejected: Position hash score below cutoff.

CGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGAG....CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGAACTGCGTGCCGAC  >  NC_014034/1253835‑1254013
                                                                                       |                                                                                                    
cGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGaa                                                                                         >  1:4041066/1‑101 (MQ=9)
 gAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGaaa                                                                                        >  1:1003715/1‑101 (MQ=9)
 gAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGaaa                                                                                        >  1:1635325/1‑101 (MQ=9)
  aaaTAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAAt                                                                                       >  2:4313225/1‑101 (MQ=9)
  aaaTAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAAt                                                                                       >  1:3249004/1‑101 (MQ=9)
   aaTAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATc                                                                                      >  1:4081788/1‑101 (MQ=9)
     tAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCct                                                                                    >  1:1872906/1‑101 (MQ=9)
      aaaCCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGGAAGCCGAGGCGAAATCctc                                                                                   >  1:1289581/1‑101 (MQ=9)
      aaaCCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCctc                                                                                   >  2:206579/1‑101 (MQ=9)
       aaCCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCctct                                                                                  >  2:848520/1‑101 (MQ=9)
        aCCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCtctc                                                                                 >  2:4833764/1‑101 (MQ=9)
          cGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAAGCGAAATCCTCTCAc                                                                               >  2:4119555/1‑101 (MQ=9)
              cccgccccgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCt                                                                           >  1:2239900/1‑101 (MQ=9)
              cccgccccgcTCTCATCTTGCCCTAAATACCTCGGGGGCGCGGGGGCGGAGACCCCGCCGGTTCGCCTCGGCGAAAACCGAAGCGAAATCCTCTCACTGCt                                                                           >  1:4191660/1‑101 (MQ=12)
                 gcccggcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTcc                                                                        >  2:576820/1‑101 (MQ=9)
                   cccgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAg                                                                      >  1:1252084/1‑101 (MQ=9)
                   cccgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGCCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAg                                                                      >  1:335838/1‑101 (MQ=9)
                    ccgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGt                                                                     >  1:3117854/1‑101 (MQ=9)
                      gcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTc                                                                   >  2:2439186/1‑101 (MQ=9)
                      gcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTc                                                                   >  1:1231573/1‑101 (MQ=9)
                        tctcATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAg                                                                 >  1:3654273/1‑101 (MQ=9)
                        tctcATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAg                                                                 >  1:4018356/1‑101 (MQ=9)
                        tctcATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAg                                                                 >  1:70412/1‑101 (MQ=9)
                         ctcATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAgg                                                                >  1:3076505/1‑101 (MQ=9)
                         ctcATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAgg                                                                >  1:258326/1‑101 (MQ=9)
                         ctcATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAgg                                                                >  2:780262/1‑101 (MQ=9)
                          tcATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGa                                                               >  2:2751823/1‑101 (MQ=9)
                          tcATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGa                                                               >  2:2965275/1‑101 (MQ=9)
                           cATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGaa                                                              >  1:4376166/1‑101 (MQ=9)
                            aTCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGGCGCCTCGGCGAAAGCCGAGGGGAAATCCTCTCACGGCTTCCCGTTCAGGaaa                                                             >  1:152997/1‑101 (MQ=9)
                             tCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAAt                                                            >  2:3355303/1‑101 (MQ=9)
                             tCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAAt                                                            >  2:3233721/1‑101 (MQ=9)
                             tCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAAt                                                            >  1:614744/1‑101 (MQ=9)
                             tCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAAt                                                            >  1:4722541/1‑101 (MQ=9)
                                                         gggccccg....cgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGAt                            <  2:3701300/94‑1 (MQ=255)
                                                            cggagcccccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGa                             >  2:4934030/10‑101 (MQ=255)
                                                             ggagcccccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGAt                            >  1:385677/9‑101 (MQ=255)
                                                             cccccccccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGAt                            <  2:1251963/93‑1 (MQ=255)
                                                              gagcccccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGAtt                           >  1:2593273/8‑101 (MQ=255)
                                                               gccccccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGAttt                          <  1:3634218/95‑1 (MQ=255)
                                                               cccccgccccGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGAttt                          <  2:1804185/93‑1 (MQ=255)
                                                               cccccccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGAttt                          <  1:4834091/95‑1 (MQ=255)
                                                               cccccccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGAttt                          <  2:4863526/95‑1 (MQ=255)
                                                               agcccccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGAttt                          >  1:1373676/7‑101 (MQ=255)
                                                               agcccccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGAttt                          >  2:188010/7‑101 (MQ=255)
                                                               agcacccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGAttt                          >  2:2105458/7‑101 (MQ=255)
                                                                ggccgtcgcGGGCCGCCTCGGCAAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTtc                         <  1:1040260/93‑1 (MQ=255)
                                                                gcccccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTtc                         >  1:260268/6‑101 (MQ=255)
                                                                ccgccgcccGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTtc                         <  2:2055303/94‑1 (MQ=255)
                                                                 ccgcccccGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTtct                        <  1:2071836/95‑1 (MQ=255)
                                                                  cccggcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTtctc                       <  1:62400/97‑1 (MQ=255)
                                                                   cccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTcc                      <  2:2239745/99‑1 (MQ=255)
                                                                   ccccccGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTcc                      <  2:2567642/97‑1 (MQ=255)
                                                                    cgcccGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCg                     <  2:4151678/98‑1 (MQ=255)
                                                                     cgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCgg                    <  2:515483/101‑1 (MQ=255)
                                                                       ccGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGGTTTCTCCGGGt                  <  2:335785/101‑1 (MQ=255)
                                                                       ccGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGt                  >  1:1803416/1‑101 (MQ=255)
                                                                       ccGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGt                  <  2:3248753/101‑1 (MQ=255)
                                                                       ccGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGt                  <  2:2386267/101‑1 (MQ=255)
                                                                        cGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTg                 >  2:289752/1‑101 (MQ=255)
                                                                        cGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTg                 <  2:2272790/101‑1 (MQ=255)
                                                                         gggTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTgg                >  2:2337890/1‑101 (MQ=255)
                                                                         gggTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTgg                >  2:4860093/1‑101 (MQ=255)
                                                                          ggTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGa               >  2:584343/1‑101 (MQ=255)
                                                                          ggTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGa               <  1:732574/101‑1 (MQ=255)
                                                                           gTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGaa              <  2:3897773/101‑1 (MQ=255)
                                                                           gTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGaa              <  1:4403056/101‑1 (MQ=255)
                                                                           gTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGaa              >  1:2527967/1‑101 (MQ=255)
                                                                            tCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGAAc             >  1:2278225/1‑101 (MQ=255)
                                                                            tCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGAAc             <  2:1872769/101‑1 (MQ=255)
                                                                            tCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGAAc             >  2:4619878/1‑101 (MQ=255)
                                                                             cGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGAACt            <  1:4863406/101‑1 (MQ=255)
                                                                             cGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGAACt            >  2:115580/1‑101 (MQ=255)
                                                                             cGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGAACt            >  2:3282935/1‑101 (MQ=255)
                                                                             cGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGAACt            >  2:1542784/1‑101 (MQ=255)
                                                                             cGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGAACt            >  2:2627540/1‑101 (MQ=255)
                                                                             cGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGAACt            >  2:4684557/1‑101 (MQ=255)
                                                                               ccgcgtcGAAAGCCGAGGCTAAAGCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGAACTGc          <  2:810197/95‑1 (MQ=255)
                                                                               ccTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGAACTGc          <  1:2446582/101‑1 (MQ=255)
                                                                               ccTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGAACTGc          <  2:2482783/101‑1 (MQ=255)
                                                                                 tCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGAACTGCGt        <  2:2496925/101‑1 (MQ=255)
                                                                                 tCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGAACTGCGt        <  1:1093445/101‑1 (MQ=255)
                                                                                 tCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGAACTGCGt        <  2:3654001/101‑1 (MQ=255)
                                                                                  cGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGAACTGCGTg       >  2:399136/1‑101 (MQ=255)
                                                                                  cGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGAACTGCGTg       >  2:2659258/1‑101 (MQ=255)
                                                                                  cGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGAACTGCGTg       >  2:3546257/1‑101 (MQ=255)
                                                                                  cGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGAACTGCGTg       >  2:2091134/1‑101 (MQ=255)
                                                                                   ggCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGAACTGCGTGc      >  2:280822/1‑101 (MQ=255)
                                                                                   ggCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGAACTGCGTGc      <  2:4507316/101‑1 (MQ=255)
                                                                                   ggCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGAACTGCGTGc      <  2:4040763/101‑1 (MQ=255)
                                                                                   ggCGAAAGCCGAAGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGAACTGCGTGc      >  1:2730628/1‑101 (MQ=255)
                                                                                    gCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGAACTGCGTGcc     <  1:2338054/101‑1 (MQ=255)
                                                                                     cGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGAACTGCGTGCcg    >  2:4794469/1‑101 (MQ=255)
                                                                                     cGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGAACTGCGTGCcg    >  1:1264335/1‑101 (MQ=255)
                                                                                     cGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGAACTGCGTGCcg    <  2:2636362/101‑1 (MQ=255)
                                                                                     cGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGAACTGCGTGCcg    <  2:3068820/101‑1 (MQ=255)
                                                                                      gAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGAACTGCGTGCcga   >  1:1431874/1‑101 (MQ=255)
                                                                                       aaaGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGAACTGCGTGCcgac  >  2:4711514/1‑101 (MQ=255)
                                                                                       aaaGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGAACTGCGTGCcgac  >  2:5029044/1‑101 (MQ=255)
                                                                                       aaaGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGAACTGCGTGCcgac  >  2:254932/1‑101 (MQ=255)
                                                                                       |                                                                                                    
CGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGAG....CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGTGGAACTGCGTGCCGAC  >  NC_014034/1253835‑1254013

Base quality scores:  ATCG  < 3 ≤  ATCG  < 22 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG