New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 = 1253880NA (NA)31 (0.150) 5/176 32.8 0.43 intergenic (+38/+53) RCAP_rcc01181/hisH hypothetical protein/imidazole glycerol phosphate synthase subunit HisH
?NC_014034 = 2577184 82 (0.350)intergenic (+106/‑323) fdx/purU ferredoxin VI/formyltetrahydrofolate deformylase
Rejected: Position hash score below cutoff.

CCCGCCGCGATGGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGAG....CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAAT  >  NC_014034/1253780‑1253955
                                                                                                    |                                                                                    
cccGCCGCGATGGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCt                                                                                      <  1:1233168/101‑1 (MQ=255)
cccGCCGCGATGGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCt                                                                                      <  1:1677177/101‑1 (MQ=255)
 ccgccgCGATGGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTc                                                                                     <  2:4488131/101‑1 (MQ=255)
 ccgccgCGATGGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTc                                                                                     <  2:3681146/101‑1 (MQ=255)
 ccgccgCGATGGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTc                                                                                     <  1:2811953/101‑1 (MQ=255)
 ccgccgCGATGGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTc                                                                                     <  2:2180426/101‑1 (MQ=255)
 ccgccgCGATGGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTc                                                                                     <  1:1947134/101‑1 (MQ=255)
    ccggGATGGGGCTGGGCAAGCGCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCggg                                                                                  <  2:2208844/101‑1 (MQ=255)
    ccgcGATGGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCggg                                                                                  >  1:335238/1‑101 (MQ=255)
    ccgcGATGGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCCCggg                                                                                  >  2:4048043/1‑101 (MQ=255)
     cgcgATGGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCgggg                                                                                 >  1:1480819/1‑101 (MQ=255)
      gcgATGGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCggggg                                                                                >  1:3426684/1‑101 (MQ=255)
      gcgATGGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCggggg                                                                                <  2:4178236/101‑1 (MQ=255)
      gcgATGGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCggggg                                                                                >  2:407227/1‑101 (MQ=255)
      gcgATGGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCggggg                                                                                >  1:2870002/1‑101 (MQ=255)
       cgATGGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGt                                                                               <  2:2546798/101‑1 (MQ=255)
       cgATGGGGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGt                                                                               <  2:1280675/101‑1 (MQ=255)
        gATGGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTg                                                                              <  2:4014192/101‑1 (MQ=255)
         aTGGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGc                                                                             <  2:3730303/101‑1 (MQ=255)
         aTGGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGc                                                                             <  2:226671/101‑1 (MQ=255)
         aTGGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGc                                                                             <  1:4096948/101‑1 (MQ=255)
         aTGGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGAGGTGc                                                                             <  2:2117473/101‑1 (MQ=255)
          tGGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCg                                                                            >  2:3120174/1‑101 (MQ=255)
          tGGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCg                                                                            >  2:3133959/1‑101 (MQ=255)
           ggTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCg.....g                                                                      <  2:2149643/101‑1 (MQ=255)
           ggTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCg.....g                                                                      <  2:3211487/101‑1 (MQ=255)
            gTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCg.....gg                                                                     <  1:1477081/101‑1 (MQ=255)
            gTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCg.....gg                                                                     <  2:1009168/101‑1 (MQ=255)
              gCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCg.....gggg                                                                   <  1:2491560/101‑3 (MQ=255)
              gCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCg.....gggg                                                                   <  2:3357661/101‑3 (MQ=255)
              gCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCg.....gggg                                                                   <  1:560487/101‑3 (MQ=255)
              gCTGGGCAAGCCCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCg.....gggg                                                                   <  1:4101229/101‑3 (MQ=255)
               cTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....c                                                              <  1:1894408/101‑1 (MQ=255)
               cTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....c                                                              <  1:4472112/101‑1 (MQ=255)
                tGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cg                                                             >  2:826320/1‑101 (MQ=255)
                cgggCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cg                                                             <  2:3855967/100‑1 (MQ=255)
                 gggCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cgg                                                            >  1:3535510/1‑100 (MQ=255)
                  ggCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cgga                                                           <  1:3904558/101‑3 (MQ=255)
                  ggCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cgga                                                           <  2:4385871/101‑3 (MQ=255)
                  ggCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cgga                                                           >  2:1550662/1‑99 (MQ=255)
                   gCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggag                                                          <  2:391647/101‑4 (MQ=255)
                    taaGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagc                                                         <  1:2121340/100‑5 (MQ=255)
                     aaGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagcc                                                        >  2:3299303/1‑96 (MQ=255)
                     aaGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagcc                                                        <  2:2834656/101‑6 (MQ=255)
                     aaGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagcc                                                        <  2:4255640/101‑6 (MQ=255)
                     aaGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagca                                                        >  2:2121190/1‑96 (MQ=255)
                      aGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagccc                                                       <  2:2022593/101‑7 (MQ=255)
                      aGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagccc                                                       >  2:1028516/1‑95 (MQ=255)
                        cTCGGCTTCCAGCCCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagccccc                                                     <  2:3792155/101‑9 (MQ=38)
                        cTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagccccc                                                     <  2:4254239/101‑9 (MQ=255)
                        cTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagccccc                                                     >  1:2540634/1‑93 (MQ=255)
                         tCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagcccccg                                                    >  2:3379154/1‑92 (MQ=38)
                             cTTCCAGCCCGGCTGCGCCTTGGAGGCGAAATAACCCGATCCCGCCCCGCTCCCATCTTCCCCAAAATCCCCCGGGGGTGCG.....GGGG....CGGCgcccccgccgg                                                <  2:2902753/101‑11 (MQ=19)
                                                     ggCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCg                                 >  1:783322/1‑101 (MQ=12)
                                                        gAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGaaa                              >  2:3503396/1‑101 (MQ=12)
                                                        gAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGaaa                              >  1:2234569/1‑101 (MQ=9)
                                                        gAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGaaa                              >  2:2113155/1‑101 (MQ=12)
                                                        gAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGGGCGGGGGCGGAGCCCCCGCCGGGTAGCCTCGGCGAAAGCCGAGGCGaaa                              >  2:2855155/1‑101 (MQ=9)
                                                         aaaTAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAAt                             >  1:3902527/1‑101 (MQ=9)
                                                           aTAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATcc                           >  2:3645734/1‑101 (MQ=9)
                                                            tAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCct                          >  1:3380305/1‑101 (MQ=9)
                                                            tAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCct                          >  1:4270036/1‑101 (MQ=9)
                                                            tAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCct                          >  1:4383723/1‑101 (MQ=9)
                                                              aaCCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCctct                        >  2:2745360/1‑101 (MQ=9)
                                                                ccGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCa                      >  2:128056/1‑101 (MQ=9)
                                                                 cGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCAc                     >  2:2252298/1‑101 (MQ=9)
                                                                 cGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGAGGAGCCCCCGCCGGGTCGCCGCGGCGAAAGCAGCGGCGAAATCCTCTCAc                     >  1:1160051/1‑101 (MQ=12)
                                                                  gATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACg                    >  1:2968515/1‑101 (MQ=9)
                                                                  gATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGGCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCCCg                    >  2:2784583/1‑101 (MQ=2)
                                                                   aTCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACgg                   >  2:2940972/1‑101 (MQ=9)
                                                                   aTCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACgg                   >  2:3473969/1‑101 (MQ=9)
                                                                    tCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGc                  >  1:776976/1‑101 (MQ=9)
                                                                    tCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGc                  >  1:4090564/1‑101 (MQ=9)
                                                                     cccgccccgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCt                 >  2:2918898/1‑101 (MQ=9)
                                                                        gccccgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTcc              >  2:621781/1‑101 (MQ=9)
                                                                        gccccgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTcc              >  2:1351082/1‑101 (MQ=9)
                                                                         ccccgcTCTCATCTTCCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCCCGGCTTCCa             >  1:1748467/1‑101 (MQ=9)
                                                                         ccacgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAGGCCGAGGCGTAATCCTCTCACGGCTTccc             >  1:2658816/4‑100 (MQ=2)
                                                                          cccgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAg            >  1:13246/1‑101 (MQ=9)
                                                                            cgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGtt          >  2:4562557/1‑101 (MQ=9)
                                                                            cgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGtt          >  1:589026/1‑101 (MQ=9)
                                                                             gcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTc         >  1:1072279/1‑101 (MQ=9)
                                                                              ctctCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCCCGGCTTCCAGTTCa        >  2:1374257/1‑101 (MQ=9)
                                                                               tctcATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAg       >  1:3239096/1‑101 (MQ=9)
                                                                               tctcATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAg       >  2:3098563/1‑101 (MQ=9)
                                                                                 tcATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGa     >  2:1546584/1‑101 (MQ=9)
                                                                                 tcATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGa     >  1:1852214/1‑101 (MQ=9)
                                                                                   aTCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGaaa   >  2:2896212/1‑101 (MQ=9)
                                                                                    tCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAAt  >  2:2188324/1‑101 (MQ=9)
                                                                                                    |                                                                                    
CCCGCCGCGATGGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGAG....CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAAT  >  NC_014034/1253780‑1253955

Base quality scores:  ATCG  < 3 ≤  ATCG  < 19 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG