New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 = 39246833 (0.160)46 (0.220) 27/174 18.1 0.73 intergenic (‑116/+27) RCAP_rcc00363/trmD hypothetical protein/tRNA (guanine‑N1)‑methyltransferase
?NC_014034 2880583 = 1 (0.000)coding (1002/1053 nt) pdxA 4‑hydroxythreonine‑4‑phosphate dehydrogenase
Rejected: Position hash score below cutoff.

GCCATTTCATCTTGCCTCAAATATCCCGGGGGTCCGGGGGCGGAGCCCCCGGCGG.................................................................................................................  >  NC_014034/392414‑392468
....................................................................CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGGC  >  NC_014034/2880583‑2880682
                                                       |||||||||||||                                                                                                    
tcctTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGC                                                                     <  2:2276628/97‑1 (MQ=37)
  ctTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCG                                                                   >  1:2986332/3‑101 (MQ=25)
  ctTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCG                                                                   >  1:3807905/3‑101 (MQ=25)
    TTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCG                                                                 >  2:2947242/1‑101 (MQ=17)
     TTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGA                                                                >  1:4091452/1‑101 (MQ=17)
      TCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGAT                                                               >  2:2107955/1‑101 (MQ=17)
        ATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCA                                                             >  2:4155814/1‑101 (MQ=17)
           TTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCCCGATCAGAC                                                          >  2:1138059/1‑101 (MQ=12)
            TGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACT                                                         >  1:1559229/1‑101 (MQ=17)
               tTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTG                                                      >  1:549121/2‑101 (MQ=11)
                TCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGC                                                     >  1:291401/1‑101 (MQ=1)
                 CAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCG                                                    >  1:886195/1‑101 (MQ=1)
                  AAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGG                                                   >  2:2213265/1‑101 (MQ=1)
                   AATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCCTCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGG                                                  >  1:428392/1‑101 (MQ=0)
                     TATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCT                                                >  2:1194213/1‑101 (MQ=1)
                       TCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGG                                              >  2:3637480/1‑101 (MQ=1)
                          CGGGGGTCGGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCT                                           <  2:1381539/101‑1 (MQ=0)
                          CGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCT                                           >  1:2826244/1‑101 (MQ=1)
                           GGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTT                                          >  1:485742/1‑101 (MQ=11)
                           GGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTT                                          >  2:213902/1‑101 (MQ=11)
                             GGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGACAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGC                                        >  2:786280/1‑101 (MQ=0)
                              GGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCC                                       >  1:193136/1‑101 (MQ=1)
                                TCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGG                                     >  1:642726/1‑101 (MQ=1)
                                 CCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGT                                    >  1:1924737/1‑101 (MQ=1)
                                  CGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTG                                   >  1:2121218/1‑101 (MQ=1)
                                   GGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTTGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGC                                  >  1:3792152/1‑101 (MQ=0)
                                   GGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGC                                  >  2:2797885/1‑101 (MQ=1)
                                     gggggcgcCCCCGGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCC                                <  1:1563982/93‑1 (MQ=255)
                                     GGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCC                                >  1:4331547/1‑101 (MQ=1)
                                        caGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCA                             >  1:2218332/3‑101 (MQ=14)
                                         aGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAA                            >  2:1123010/2‑101 (MQ=14)
                                            cCCCCGGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGT                         <  1:3841074/100‑1 (MQ=17)
                                             cccggGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTC                        <  1:2456568/96‑1 (MQ=17)
                                              CCCGGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCA                       <  2:4406201/101‑1 (MQ=17)
                                               CCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCAT                      >  1:3743941/1‑101 (MQ=32)
                                               CCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCAT                      >  1:3678501/1‑101 (MQ=32)
                                                CCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATA                     <  2:549073/101‑1 (MQ=32)
                                                CCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATA                     >  1:1283564/1‑101 (MQ=32)
                                                 CGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAG                    <  2:359677/101‑1 (MQ=32)
                                                  GGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGG                   >  2:3976101/1‑101 (MQ=32)
                                                  GGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGG                   >  2:1897858/1‑101 (MQ=32)
                                                   GCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGC                  >  2:3857720/1‑101 (MQ=32)
                                                   GCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGC                  <  2:2686490/101‑1 (MQ=32)
                                                    CGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCC                 <  2:642665/101‑1 (MQ=32)
                                                    CGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCC                 >  2:801628/1‑101 (MQ=32)
                                                    CGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCC                 <  1:14945/101‑1 (MQ=32)
                                                      gccCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGT               <  2:1910149/98‑1 (MQ=32)
                                                       CGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTG              >  1:1359278/1‑101 (MQ=31)
                                                       CGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTG              <  2:1559124/101‑1 (MQ=31)
                                                          GTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCA           <  2:4331255/101‑1 (MQ=1)
                                                          GTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCA           >  2:3926327/1‑101 (MQ=1)
                                                          GTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCA           >  2:1822218/1‑101 (MQ=1)
                                                           TCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCAT          <  1:3968534/101‑1 (MQ=1)
                                                             ACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGA        >  2:1821586/1‑101 (MQ=1)
                                                               AGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATC      <  1:4469748/101‑1 (MQ=1)
                                                                GCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCG     <  2:367633/101‑1 (MQ=1)
                                                                 CCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGG    >  1:2921637/1‑101 (MQ=1)
                                                                  CGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGG   >  1:3047942/1‑101 (MQ=1)
                                                                  CGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGG   >  1:1855608/1‑101 (MQ=1)
                                                                  CGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGG   >  2:683388/1‑101 (MQ=1)
                                                                  CGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGG   >  1:1126798/1‑101 (MQ=1)
                                                                  CGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCAATGATCGGG   >  1:784799/1‑101 (MQ=1)
                                                                   GCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGGC  <  2:1623910/101‑1 (MQ=1)
                                                                   GCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGGC  <  1:2798101/101‑1 (MQ=1)
                                                       |||||||||||||                                                                                                    
GCCATTTCATCTTGCCTCAAATATCCCGGGGGTCCGGGGGCGGAGCCCCCGGCGG.................................................................................................................  >  NC_014034/392414‑392468
....................................................................CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGGC  >  NC_014034/2880583‑2880682

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