New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1955860 =25 (0.490)7 (0.140) 3/58 NT 15.9% coding (50/888 nt) motA flagellar motor protein MotA
?REL606 = 4141016 49 (0.790)coding (89/705 nt) yijC DNA‑binding transcriptional repressor

AGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCT  >  REL606/1955825‑1955895
                                   |                                   
agagTGCTCCAAGGCTTCCACCGGTCATCAAATAAc                                     <  1:6138057/36‑1 (MQ=255)
 gattGCTCCAAGGCTTCAACCGGTCATCAAATAAcc                                    <  1:857919/33‑1 (MQ=37)
 gagTGCTCCAAGGCTTCCACCGGTCATCAAATAAcc                                    <  1:7177741/36‑1 (MQ=255)
 gagTGCTCCAAGGCTTCCACCGGTCATCAAATAAcc                                    >  1:2198608/1‑36 (MQ=255)
 gagTGCTCCAAGGCTTCCACCGGTCATCAAATAAcc                                    <  1:4478368/36‑1 (MQ=255)
 gagTGCTCCAAGGCTTCCACCGGTCATAAAATAAcc                                    >  1:4168769/1‑36 (MQ=255)
   gTGCTCCAAGGCTTCCACCGGTCATCAAATAAccgc                                  <  1:5887129/36‑1 (MQ=255)
      cTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGa                               >  1:7397683/1‑36 (MQ=255)
      cTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGa                               >  1:5927709/1‑36 (MQ=255)
        ccAAGGCTTCCACCGGTCATCAAATAACCGCCGaaa                             >  1:805552/1‑36 (MQ=255)
        ccAAGGCTTCCACCGGTCATCAAATAACCGCCGaaa                             >  1:4553757/1‑36 (MQ=255)
        caaaGCCTCCCACAGGTAAAAAAAAAACCGCCGaaa                             <  1:222562/34‑1 (MQ=0)
         cAAGGCTTCCACCGGTCATCAAATAACCGCCGaaaa                            >  1:5323910/1‑36 (MQ=255)
             gCTTCCACCGGTCATCAAATAACCGCCTAAAACTGt                        >  1:8177049/1‑36 (MQ=255)
             gCTTCCACCGGTCATCAAATAACCGCCGAAAACTGt                        >  1:3417305/1‑36 (MQ=255)
              cTTCCACCGGTCATCAAATAACCGCCGAAAACTGTa                       >  1:7639050/1‑36 (MQ=255)
                 ccACCGGTCATCAAATAACCGCCGAAAACTGTACCg                    >  1:1293811/1‑36 (MQ=255)
                  cACCGGTCATCAAATAACCGCCGAAAACTGTACCga                   >  1:3347377/1‑36 (MQ=255)
                   aCCGGTCATCAAATAACCGCCGAAAACTGTACCgag                  >  1:4583578/1‑36 (MQ=255)
                    ccGGTCATCAAATAACCGCCGAAAACTGTACCgaga                 >  1:638923/1‑36 (MQ=255)
                       gaaaaaaaaaacCCGCCGAAAACTGTACCGAGaaca              >  1:1713578‑M1/13‑36 (MQ=255)
                       gTCATCAAATAACCGCCGAAAACTGTACCGAGaaca              >  1:6488897/1‑36 (MQ=255)
                       gTCATCAAAGAAAAGCCGAATACTGTACCGAGaaca              >  1:7603105/1‑36 (MQ=19)
                       aaaaaaaaaaacCCGCCGAAAACTGTACCGAGaaca              >  1:2798616‑M1/13‑36 (MQ=255)
                        tcatcaAATAACCGCCGAAAACTGTACCGAGaacaa             >  1:7449109/1‑36 (MQ=255)
                        tcatcaAATAACCGCCGAAAACTGTACCGAGaacaa             >  1:4992618/1‑36 (MQ=255)
                        aaaaaaaaaacCCGCCGAAAACTGTACCGAGaacaa             >  1:4701667‑M1/12‑36 (MQ=255)
                         catcaAATAACCGCCGAAAACTGTACCGAGaacaac            >  1:2849193/1‑36 (MQ=255)
                          aaaaaaAAACCGCCGAAAACTGTACCGAGAACAAcc           <  1:7586754/33‑1 (MQ=37)
                           tcaAATAACCGCCGAAAACTGTACAGAGAACAACCa          >  1:8473146/1‑36 (MQ=255)
                           aaaaaaacCCGCCGAAAACTGTACCGAGAACAACCa          >  1:6338915‑M1/9‑36 (MQ=255)
                            aaaaTAACCGCCGAAAACTGTACCGAGAACAACCAg         >  1:2280523/2‑36 (MQ=255)
                             aaaTAACCGCCGAAAACTGTACCGAGAACAACCAgg        >  1:6381845/1‑36 (MQ=255)
                              aaaaaCCGCCGAAAACTGTACCGAGAACAACCAGGt       >  1:4500260‑M1/6‑36 (MQ=255)
                              aaaaaCCGCCGAAAACTGTACCGAGAACAACCAGGt       >  1:6997068‑M1/6‑36 (MQ=255)
                              aaTAACCGCCGAAAACTGTACCGAGAACAACCAGGt       >  1:8002876/1‑36 (MQ=255)
                               aaacCCGCCGAAAACTGTACCGAGAACAACCAGGTa      >  1:5485810‑M1/5‑36 (MQ=255)
                               aTAACCGCCGAAAACTGTACCGAGAACAACCAGGTa      >  1:1922767/1‑36 (MQ=255)
                                 aaCCGCCGAAAACTGTACCGAGAACAACCAGGTAAc    >  1:6901818/1‑36 (MQ=255)
                                   CCGCCGAAAACTGTACCGAGAACAACCAGGTAACCt  <  1:2151918‑M1/36‑1 (MQ=255)
                                   CCGCCGAAAACTGTACCGAGAACAACCAGGTAACCt  <  1:7707359‑M1/36‑1 (MQ=255)
                                   CCGCCGAAAACTGTACCGAGAACAACCAGGTAACCt  <  1:7971511‑M1/36‑1 (MQ=255)
                                   |                                   
AGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCT  >  REL606/1955825‑1955895

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG