New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 = 1748592139 (0.660)45 (0.210) 3/186 33.7 0.39 intergenic (‑296/‑102) RCAP_rcc01616/gppA hypothetical protein/guanosine‑5'‑triphosphate,3'‑diphosphate pyrophosphatase
?NC_014034 = 2600942 1 (0.000)intergenic (+4/‑137) RCAP_rcc02424/surE short‑chain dehydrogenase/reductase family oxidoreductase/5'‑nucleotidase SurE
Rejected: Position hash score below cutoff.

GCGGAGCCTCCGGCGGAGGTATTTAGGGCAAGATGAAACGGTGCCCCTTCATCTTGCCTCAAATATCCCGGGGGGAGGCCGCAGGCCCGGGGGCAGCGCCCCCCG....................................................................................................  >  NC_014034/1748488‑1748592
..................................................................................................CCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGAAGCGCGCCCGCACGGCGGTGGC  <  NC_014034/2600942‑2600843
                                                                                                                                                                                                             
                                                                               gcgaaGCCGGGGGGCAGAGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGA                           <  2:1821377/96‑1 (MQ=1)
                                                                                cgaaGCCGGGGGGCAGAGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAA                          >  2:1439911/5‑101 (MQ=1)
                                                                                cgaaGCCGGGGGGCAGAGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAA                          <  2:2962821/97‑1 (MQ=1)
                                                                                 gaaGCCGGGGGGCAGAGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAG                         >  1:1007707/4‑101 (MQ=1)
                                                                                 gaaGCCGGGGGGCAGAGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAG                         >  2:106601/4‑101 (MQ=1)
                                                                                 gaaGCCGGGGGGCAGAGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAG                         >  1:1632469/4‑101 (MQ=1)
                                                                                 gaaGCCGGGGGGCAGAGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAG                         >  1:1865277/4‑101 (MQ=1)
                                                                                 gaaGCCGGGGGGCAGAGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAG                         >  2:482636/4‑101 (MQ=1)
                                                                                 gaaGCCGGGGGGCAGAGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAG                         >  2:3846544/4‑101 (MQ=1)
                                                                                  aaGCCGGGGGGCAGAGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGA                        >  2:3737302/3‑101 (MQ=1)
                                                                                  aaGCCGGGGGGCAGAGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGA                        <  1:3352598/99‑1 (MQ=1)
                                                                                  aaGCCGGGGGGCAGAGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGA                        <  1:2768641/99‑1 (MQ=1)
                                                                                  aaGCCGGGGGGCAGAGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGA                        >  1:2736017/3‑101 (MQ=1)
                                                                                  aaGCCGGGGGGCAGAGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGA                        <  1:1545095/99‑1 (MQ=1)
                                                                                  aaGCCGGGGGGCAGAGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGA                        >  1:2339738/3‑101 (MQ=1)
                                                                                   aGCCGGGGGGCAGAGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGAA                       >  2:1781560/2‑101 (MQ=1)
                                                                                   aGCCGGGGGGCAGAGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGAA                       <  2:1042565/100‑1 (MQ=1)
                                                                                   aGCCGGGGGGCAGAGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGAA                       >  1:2653115/2‑101 (MQ=1)
                                                                                   aGCCGGGGGGCAGACCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGAA                       <  2:2587881/100‑1 (MQ=1)
                                                                                    GCCGGGGGGCAGAGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGAAG                      <  2:1154956/101‑1 (MQ=1)
                                                                                    GCCGGGGGGCAGAGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGAAG                      >  1:3452705/1‑101 (MQ=1)
                                                                                     ccgGGGGGCAGAGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGAAGC                     <  1:981544/98‑1 (MQ=1)
                                                                                      ngGGGGGCAGAGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGANNANGCn                    >  1:1831462/3‑100 (MQ=1)
                                                                                      ggGGGGGCAGCCCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGAAGCG                    <  1:1551736/99‑1 (MQ=1)
                                                                                      cgGGGGGCAGAGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGAAGCG                    <  2:1277049/99‑1 (MQ=1)
                                                                                      cgGGGGGCAGAGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGAAGCG                    >  2:1963757/3‑101 (MQ=1)
                                                                                      cgGGGGGCAGAGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGCCCTCGGACGGGAAGGAGCG                    >  1:1249898/3‑101 (MQ=0)
                                                                                       gGGGGGCAGAGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGAAGCGC                   >  1:3456140/2‑101 (MQ=1)
                                                                                       gGGGGGCAGAGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGAAGCGC                   >  2:1221863/2‑101 (MQ=1)
                                                                                        GGGGGCAGAGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGAAGCGCG                  >  1:1298877/1‑101 (MQ=1)
                                                                                        GGGGGCAGAGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGAAGCGCG                  <  2:2177369/101‑1 (MQ=1)
                                                                                        GGGGGCAGAGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGAAGCGCG                  <  2:2793324/101‑1 (MQ=1)
                                                                                         GGGGCAGAGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGAAGCGCGC                 >  1:3139658/1‑101 (MQ=1)
                                                                                          GGGCAGAGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGAAGCGCGCC                <  1:2437751/101‑1 (MQ=1)
                                                                                           GGCAGCCCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGAAGCGCGCCC               <  1:2975436/101‑1 (MQ=1)
                                                                                            ggccCGCCCCCGGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGAAGCGCGCCCG              <  2:876355/97‑1 (MQ=0)
                                                                                              agaGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGAAGCGCGCCCGCA            >  2:397666/4‑101 (MQ=1)
                                                                                              agaGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGAAGCGCGCCCGCA            <  2:4003418/98‑1 (MQ=1)
                                                                                               gaGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGAAGCGCGCCCGCAC           <  1:2478824/99‑1 (MQ=1)
                                                                                               gaGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGAAGCGCGCCCGCAC           <  2:4286471/99‑1 (MQ=1)
                                                                                               gaGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGAAGCGCGCCCGCAC           >  1:1921924/3‑101 (MQ=1)
                                                                                               gaGCCCCCCGTCCGTCAGATCCCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGAAGCGCGCCCGCAC           <  1:570262/99‑1 (MQ=1)
                                                                                                aGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGAAGCGCGCCCGCACG          >  2:2250945/2‑101 (MQ=1)
                                                                                                aGCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGAAGCGCGCCCGCACG          >  1:3380892/2‑101 (MQ=1)
                                                                                                 GCCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGAAGCGCGCCCGCACGG         >  1:203381/1‑101 (MQ=1)
                                                                                                                                                                                                             
GCGGAGCCTCCGGCGGAGGTATTTAGGGCAAGATGAAACGGTGCCCCTTCATCTTGCCTCAAATATCCCGGGGGGAGGCCGCAGGCCCGGGGGCAGCGCCCCCCG....................................................................................................  >  NC_014034/1748488‑1748592
..................................................................................................CCCCCCGTCCGTCAGATCTCGGCAAGAAGGCGTTTCGCTTCGTCGCGGCAATCGATCAGCGGCGTGCGGTCCTCGGACGGGAAGAAGCGCGCCCGCACGGCGGTGGC  <  NC_014034/2600942‑2600843

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