New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 = 1397797133 (0.590)39 (0.170) 4/156 36.1 0.37 intergenic (+101/+1797) RCAP_rcc01307/RCAP_rcc01309 hypothetical protein/IS66 family transposase
?NC_014034 1434665 = NA (NA)intergenic (‑444/‑514) RCAP_rcc01334/RCAP_rcc01335 hypothetical protein/Fis family GAF modulated sigma54 specific transcriptional regulator
Rejected: Position hash score below cutoff.

GTGGATCGAAACGACGCGCTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAA  >  NC_014034/1434567‑1434709
                                                                                                  |                                            
gTGGATCGAAACGACGCGCTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAggg                                            <  2:2913302/101‑1 (MQ=255)
 tGGATCGAAACGACGCGCTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGa                                           >  2:3856487/1‑101 (MQ=255)
 tGGATCGAAACGACGCGCTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGa                                           <  2:1549332/101‑1 (MQ=255)
 tGGATCGAAACGACGCGCTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGa                                           <  2:4976693/101‑1 (MQ=255)
  ggATCGAAACGACGCGCTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGAc                                          >  2:3318786/1‑101 (MQ=255)
  ggATCGAAACGACGCGCTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGAc                                          <  2:5358107/101‑1 (MQ=255)
   gATCGAAACGACGCGCTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACa                                         >  2:2333920/1‑101 (MQ=255)
   gATCGAAACGACGCGCTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACa                                         <  1:380202/101‑1 (MQ=255)
   gATCGAAACGACGCGCTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACa                                         >  1:2390294/1‑101 (MQ=255)
    aTCGAAACGACGCGCTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAg                                        <  2:922381/101‑1 (MQ=255)
    aTCGAAACGACGCGCTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAg                                        >  1:1913134/1‑101 (MQ=255)
       gAAACGACGCGCTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCt                                     >  2:3204082/1‑101 (MQ=255)
       gAAACGACGCGCTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCt                                     >  1:5062493/1‑101 (MQ=255)
        aaaCGACGCGCTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCtt                                    >  1:3200485/1‑101 (MQ=255)
         aaCGACGCGCTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCttt                                   <  2:3997968/101‑1 (MQ=255)
         aaCGACGCGCTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCttt                                   <  1:4844753/101‑1 (MQ=255)
          acggcgCGCTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTg                                  <  2:5370642/101‑1 (MQ=255)
          acgacgCGCTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGGGCCGAGGGACAGGCTTTg                                  >  2:1106905/1‑101 (MQ=255)
          acgacgCGCTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTg                                  <  1:750115/101‑1 (MQ=255)
           cgacgCGCTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGc                                 <  1:2188122/101‑1 (MQ=255)
           cgacgCGCTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGc                                 <  2:3946079/101‑1 (MQ=255)
           cgacgCGCTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGc                                 <  2:3162634/101‑1 (MQ=255)
           cgacgCGCTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGc                                 >  1:4676328/1‑101 (MQ=255)
            gacgCGCTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCa                                <  2:1808831/101‑1 (MQ=255)
            gacgCGCTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCa                                <  1:3552211/101‑1 (MQ=255)
             acgCGCTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAg                               <  2:910638/101‑1 (MQ=255)
             acgCGCTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAg                               >  2:4960275/1‑101 (MQ=255)
              cgcgCTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGt                              >  1:589437/1‑101 (MQ=255)
              cgcgCGGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGt                              <  1:2154226/101‑1 (MQ=255)
               gcgcTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTc                             <  1:4403868/101‑1 (MQ=255)
                cgcTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCg                            <  1:449929/101‑1 (MQ=255)
                cgcTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCg                            <  1:831675/101‑1 (MQ=255)
                cgcTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCg                            >  2:2786139/1‑101 (MQ=255)
                cgcTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCg                            >  2:578864/1‑101 (MQ=255)
                 gcTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGc                           >  2:3396345/1‑101 (MQ=255)
                  cTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCt                          >  1:4356605/1‑101 (MQ=255)
                  cTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCt                          <  1:1332817/101‑1 (MQ=255)
                  cTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCt                          <  1:453900/101‑1 (MQ=255)
                   tGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTc                         <  1:4538634/101‑1 (MQ=255)
                   tGAAGGAGGGTGATCGCGGTGGCGCCTTGTNGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTc                         <  1:452516/101‑1 (MQ=255)
                    gAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTcc                        <  2:1209382/101‑1 (MQ=255)
                    gAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTcc                        <  1:198878/101‑1 (MQ=255)
                    gAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCCCTcc                        >  2:2850474/1‑101 (MQ=255)
                     aaGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTccc                       <  2:5695271/101‑1 (MQ=255)
                     aaGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTccc                       <  1:3764385/101‑1 (MQ=255)
                      aggaggGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTcccc                      <  1:3811673/101‑1 (MQ=255)
                      aggaggGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTcccc                      <  1:1467059/101‑1 (MQ=255)
                        gaggGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTcccccc                    <  1:4683604/101‑1 (MQ=255)
                         aggGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCcg                   >  1:2394283/1‑101 (MQ=255)
                         aggGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGCCGCTccccccc                   >  2:4933451/1‑100 (MQ=255)
                          gggTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGACGGCCCGGCTTTGCAGTAGCTcccccacc                  >  2:1938258/1‑98 (MQ=17)
                           ggTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCcgcg                 >  2:5344472/1‑101 (MQ=255)
                           ggTGATCTCGGTGCCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCcgcg                 <  2:5324503/101‑1 (MQ=255)
                            gTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCgg                <  1:1598514/101‑1 (MQ=255)
                            gTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCACCCCGCgg                <  1:3487980/101‑1 (MQ=255)
                             tGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCggg               >  1:2664622/1‑101 (MQ=255)
                             tGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCggg               >  2:802733/1‑101 (MQ=255)
                              gATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCgggg              <  1:2334046/101‑1 (MQ=255)
                               aTCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCggggg             >  2:285710/1‑101 (MQ=255)
                               aTCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCggggg             >  2:430461/1‑101 (MQ=255)
                                tctcGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGc            >  1:3361677/1‑101 (MQ=40)
                                tctcGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGc            <  1:2776996/101‑1 (MQ=40)
                                tctcGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGc            >  1:744577/1‑101 (MQ=40)
                                tctcGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGc            <  2:1104038/101‑1 (MQ=40)
                                  tcGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGCGt          >  2:1539944/1‑101 (MQ=39)
                                   cGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGCggg         >  2:5658398/1‑99 (MQ=38)
                                    ggTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGCGTgg        >  1:1874022/1‑101 (MQ=37)
                                    ggTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGCGTgg        >  1:1925471/1‑101 (MQ=37)
                                    ggTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGCGTgg        <  1:2737474/101‑1 (MQ=37)
                                    ggTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGCGTgg        >  1:3561625/1‑101 (MQ=37)
                                    ggTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGCGTgg        <  2:1045801/101‑1 (MQ=37)
                                    ggTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGCGTgg        <  1:4548649/101‑1 (MQ=37)
                                    ggTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCCCGGGGGCGCgg        >  2:837618/1‑101 (MQ=255)
                                     gTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGCGTGGa       <  2:2394149/101‑1 (MQ=35)
                                     gTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGCGTGGa       <  1:4313697/101‑1 (MQ=35)
                                     gTGGCGCCTTGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGCGTGGa       <  2:3848059/101‑1 (MQ=25)
                                      tGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGCGTGGAt      <  1:3856729/101‑1 (MQ=34)
                                       ggCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGCGTGGATc     >  1:5533895/1‑101 (MQ=33)
                                       ggCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGCGTGGATc     <  2:2602616/101‑1 (MQ=33)
                                        gcgcCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGCGTGGATCg    >  2:5216212/1‑101 (MQ=32)
                                        gcgcCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGCGTGGATCg    <  2:1276064/101‑1 (MQ=32)
                                        gcgcCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGCGTGGATCg    <  2:359313/101‑1 (MQ=32)
                                         cgcCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGa   >  2:825359/1‑101 (MQ=31)
                                         cgcCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGa   >  1:1087582/1‑101 (MQ=31)
                                         cgcCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGCCTGGATCGa   >  2:3541562/1‑101 (MQ=14)
                                          gcCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaa  >  1:4554949/1‑101 (MQ=11)
                                                                                                  |                                            
GTGGATCGAAACGACGCGCTGAAGGAGGGTGATCTCGGTGGCGCCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAA  >  NC_014034/1434567‑1434709

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG