New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 = 125586524 (0.140)60 (0.350) 27/146 14.7 0.39 intergenic (+35/+207) RCAP_rcc01184/RCAP_rcc01185 CsbD family protein/methyl‑accepting chemotaxis sensory transducer
?NC_014034 = 2486558 162 (0.680)intergenic (‑186/+8) RCAP_rcc02312/soxG hypothetical protein/sarcosine oxidase subunit gamma
Rejected: Position hash score below cutoff.

CAGAAATACGGCCGCACCAAGGAGGCGGCGAAGAAGGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GG.....CTG.TC.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..TGGTAC..CCCTGCGGGACTCCCCCCGGGATATTTCCGCCAA  >  NC_014034/1255765‑1255938
                                                                                                     |                                                                                    
cagAAATACGGCCGCACCAAGGAGGCGGCGAAGAAGGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.gg                                                                                      >  2:2557052/1‑101 (MQ=255)
cagAAATACGGCCGCACCAAGGAGGCGGCGAAGAAGGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.gg                                                                                      >  1:2407627/1‑101 (MQ=255)
  gAAATACGGCCGCACCAAGGAGGCGGCGAAGAAGGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.gggc                                                                                    >  2:2938298/1‑99 (MQ=255)
   aaaTACGGCCGCACCAAGGAGGCGGCGAAGAAGGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.gggcc                                                                                   >  2:3167502/1‑98 (MQ=255)
   aaaTACGGCCGCACCAAGGAGGCGGCGAAGAAGGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.gggcc                                                                                   <  2:353745/101‑4 (MQ=255)
   aaaTACGGCCGCACCAAGGAGGCGGCGAAGAAGGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.gggcc                                                                                   >  1:2111390/1‑98 (MQ=255)
   aaaTACGGCCGCACCAAGGAGGCGGCGAAGAAGGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.gggcc                                                                                   >  1:524709/1‑98 (MQ=255)
   aaaTACGGCCGCACCAAGGAGGCGGCGAAGAAGGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.gggcc                                                                                   <  1:516183/101‑4 (MQ=255)
    aaTACGGCCGCACCAAGGAGGCGGCGAAGAAGGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.gggccc                                                                                  <  1:1456305/101‑5 (MQ=255)
    aaTACGGCCGCACCAAGGAGGCGGCGAAGAAGGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.gggccc                                                                                  <  2:4664679/101‑5 (MQ=255)
     aTACGGCCGCACCAAGGAGGCGGCGAAGAAGGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.gggcccc                                                                                 <  2:252505/101‑6 (MQ=255)
     aTACGGCCGCACAAAGGAGGCGGCGAAGAAGGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.gggcccc                                                                                 >  1:2529940/1‑96 (MQ=255)
      aaCGGCCGCACCAAAGAGGCCGCGAAGAAGGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.gggccccc                                                                                <  2:1012307/100‑7 (MQ=255)
       aCGGCCGCACCAAGGAGGCGGCGAAGAAGGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.gggcccccg                                                                               <  2:2622116/101‑8 (MQ=255)
       aCGGCCGCACCAAGGAGGCGGCGAAGAAGGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.gggcccccg                                                                               <  2:2846657/101‑8 (MQ=255)
       aCGGCCGCACCAAGGAGGCGGCGAAGAAGGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.gggcccccg                                                                               <  2:532773/101‑8 (MQ=255)
       aCGGCCGCACCAAGGAGGCGGCGAAGAAGGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.gggcccccg                                                                               >  1:2410258/1‑94 (MQ=255)
       aCGGCCGCACCAAGGAGGCGGCGAAGAAGGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.gggcccccg                                                                               <  1:2210507/101‑8 (MQ=255)
       aCGGCCGCACCAAGGAGGCGGCGAAGAAGGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.gggcccccg                                                                               >  1:2204400/1‑94 (MQ=255)
         ggCCGCACCAAGGAGGCGGCGAAGAAGGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.gggcccccgct                                                                             <  2:3132219/101‑10 (MQ=255)
         ggCCGCACCAAGGAGGCGGCGAAGAAGGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.gggcccccgct                                                                             <  1:3736971/101‑10 (MQ=255)
         ggCCGCACCAAGGAGGCGGCGAAGAAGGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.gggcccccgct                                                                             >  1:4605002/1‑92 (MQ=255)
          gCCGCACCAAGGAGGCGGCGAAGAAGGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.gggcccccgctt                                                                            <  2:470742/101‑11 (MQ=255)
          gCCGCACCAAGGAGGCGGCGAAGAAGGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.gggcccccgctt                                                                            <  1:2529107/101‑11 (MQ=255)
          gCCGCACCAAGGAGGCGGCGAAGAAGGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.gggcccccgctt                                                                            >  2:1680391/1‑91 (MQ=255)
          gCCGCACCAAGGAGGCGGCGAAGAAGGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.gggcccccgctt                                                                            <  1:1651661/101‑11 (MQ=255)
                              aagaagGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAgggg                                                      >  1:3086497/1‑101 (MQ=255)
                              aagaagGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAgggg                                                      >  2:895453/1‑101 (MQ=255)
                               agaagGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAggggg                                                     >  1:2999180/1‑101 (MQ=255)
                                gaaggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAgggggg                                                    <  1:662126/101‑1 (MQ=255)
                                gaaggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGGGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAgggggg                                                    >  2:4591084/1‑101 (MQ=255)
                                 aaggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAggggggg                                                   >  1:2622892/1‑101 (MQ=255)
                                 aaggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAggggggg                                                   <  2:2126450/101‑1 (MQ=255)
                                 aaggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAggggggg                                                   >  2:1777990/1‑101 (MQ=255)
                                 aaggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAggggggg                                                   <  1:651559/101‑1 (MQ=255)
                                  aggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAgggggggg                                                  >  2:3268860/1‑100 (MQ=255)
                                  aggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGgc                                                  >  2:404222/1‑101 (MQ=255)
                                  aggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGgc                                                  <  2:4528672/101‑1 (MQ=255)
                                   ggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGgcg                                                 <  2:3371638/101‑1 (MQ=255)
                                   ggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGgcg                                                 <  2:4015066/101‑1 (MQ=255)
                                   ggaagTAGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGgcg                                                 <  1:1108747/101‑1 (MQ=255)
                                     aagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGgcgcc                                               >  1:1630644/1‑100 (MQ=255)
                                     aagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCt                                               <  1:889204/101‑1 (MQ=255)
                                     aagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCt                                               <  2:4235542/101‑1 (MQ=255)
                                     aagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCt                                               <  1:2741890/101‑1 (MQ=255)
                                      agTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtg                                              >  2:1633078/1‑101 (MQ=255)
                                      agTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtg                                              <  2:2345220/101‑1 (MQ=255)
                                         cgacgaCTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtgccc                                           <  2:410178/101‑4 (MQ=255)
                                         cgacgaCTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtgccc                                           <  2:3943798/101‑4 (MQ=255)
                                          gacgacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtgcccc                                          >  1:3314977/1‑97 (MQ=255)
                                          gacgacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtgcccc                                          >  1:2461293/1‑97 (MQ=255)
                                            cgacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtgcccccc                                        <  1:3720408/101‑7 (MQ=255)
                                            cgacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtgcccccc                                        <  2:2111240/101‑7 (MQ=255)
                                             gccTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGGGCCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtgccccccg                                       <  1:2188219/99‑8 (MQ=255)
                                             gacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCCCCGG.TT.TCGGTCCGCCCGGCGAGGGGGGGTGCtctcccccc                                       >  1:1409147/1‑93 (MQ=255)
                                             gacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtgccccccg                                       >  1:660020/1‑94 (MQ=255)
                                             gacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtgccccccg                                       >  1:2967542/1‑94 (MQ=255)
                                             gacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGGCCGCCCGAAGAGGGGGGGCGCtgccccccg                                       <  1:1791089/101‑8 (MQ=255)
                                              acTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCAGGTCC..CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtgccccccac                                      >  1:3766505/1‑93 (MQ=255)
                                              acTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtgccccccgt                                      >  1:914153/1‑93 (MQ=255)
                                              acTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtgccccccgt                                      <  1:1149137/101‑9 (MQ=255)
                                               cTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtgccccccgtc                                     <  2:3154124/101‑10 (MQ=255)
                                               cTTCCTGTCGAAACTCTGACCAAGGGCTGCGGGACCGGAACCGGGGGCTTTC.GGGCCC.CCGCTT.TAGGTCCGCCCGAAGAGGGGGGGCGCtgccccccgtc                                     <  2:4461020/101‑10 (MQ=255)
                                                ttCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtgccccccgtcc                                    >  1:1421850/1‑91 (MQ=255)
                                                ttCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtgccccccgtcc                                    <  1:2454384/101‑11 (MQ=255)
                                                ttCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtgccccccgtcc                                    <  1:3167746/101‑11 (MQ=255)
                                                                ggccTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCG.CCTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTcccccc                  <  1:2810349/99‑1 (MQ=2)
                                                                gACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTcccccc                  <  2:2817969/101‑1 (MQ=9)
                                                                 gccTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCg                 <  2:2998951/100‑1 (MQ=2)
                                                                 aCCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCg                 <  2:262693/101‑1 (MQ=9)
                                                                   cTAGGGGTGCCGGACCGGAACGGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGC..CCCTGCGGGACTCCCCCCggg               >  2:169957/1‑101 (MQ=2)
                                                                    tAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGa              >  2:3177762/1‑101 (MQ=2)
                                                                    tAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGa              >  2:3257410/1‑101 (MQ=2)
                                                                    tAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGa              >  1:2740666/1‑101 (MQ=2)
                                                                     aGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGat             <  2:3528564/101‑1 (MQ=2)
                                                                      ggggTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCgggcaa            >  2:4286807/1‑98 (MQ=2)
                                                                       gggTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGatat           <  2:2204250/101‑1 (MQ=2)
                                                                         gTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATAttt         <  2:4666616/101‑1 (MQ=0)
                                                                         gTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATAttt         <  1:2707785/101‑1 (MQ=0)
                                                                          tGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTc        >  2:409650/1‑101 (MQ=2)
                                                                            ccggaccggaACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCg      <  1:1997034/101‑1 (MQ=2)
                                                                             cggaccggaACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGCTATTTCCGc     >  1:1594806/1‑101 (MQ=2)
                                                                             cggaccggaACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGc     <  1:90527/101‑1 (MQ=2)
                                                                             cggaccggaACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGc     >  1:260502/1‑101 (MQ=2)
                                                                              ggaccggaACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGcc    <  1:2163655/101‑1 (MQ=2)
                                                                                accggaACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCaa  <  2:3490024/101‑1 (MQ=2)
                                                                                accggaACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCaa  <  1:2491559/101‑1 (MQ=2)
                                                                                accggaACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCaa  <  2:2536190/101‑1 (MQ=2)
                                                                                                     |                                                                                    
CAGAAATACGGCCGCACCAAGGAGGCGGCGAAGAAGGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GG.....CTG.TC.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..TGGTAC..CCCTGCGGGACTCCCCCCGGGATATTTCCGCCAA  >  NC_014034/1255765‑1255938

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG