New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1286705 =26 (0.540)8 (0.170) 4/56 NT 27.1% intergenic (+29/‑154) narI/ECB_01206 nitrate reductase 1, gamma (cytochrome b(NR)) subunit/hypothetical protein
?REL606 3092577 = 17 (0.280)coding (756/927 nt) yqhG hypothetical protein

GTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAA  >  REL606/3092542‑3092612
                                   |                                   
gTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGAt                                     >  1:5316007/1‑36 (MQ=255)
 ttGATGAAATTCGAGTCCTGGAGAAAACCGCCGAtt                                    >  1:5512097/1‑36 (MQ=255)
 ttGATGAAATTCGAGTCCTGGAGAAAACCGCCGAtt                                    >  1:7919234/1‑36 (MQ=255)
 ttGATAAAATTCGATTCCTGGAAAAAACCGCCGAtt                                    <  1:1791/36‑1 (MQ=21)
  tgatgaAATTCGAGTCCTGGAGAAAACCGCCGATTc                                   >  1:5075734/1‑36 (MQ=255)
  tgatgaAATTCGAGTCCTGGAGAAAACCGCCGATTc                                   >  1:7623664/1‑36 (MQ=255)
   gatgaAATTCGAGTCCTGGAGAAAACCGCAGATTcc                                  >  1:837819/1‑36 (MQ=255)
    atgaAATTCGAGTCCTGGAGAAAACCGCCGATTccg                                 >  1:6949889/1‑36 (MQ=255)
     tgaAATTCGAGTCCTGGAGAAAACCGCCGATTccgc                                >  1:6420046/1‑36 (MQ=255)
      gaAATTCGAGTCCTGGAGAAAACCGCCGATTccgcc                               >  1:6805454/1‑36 (MQ=255)
      gaAATTCGAGTCCTGGAGAAAACCGCCGATTccgcc                               >  1:7423489/1‑36 (MQ=255)
      gaAATTCGAGTCCTGGAGAAAACCGCCGATTccgcc                               >  1:1998049/1‑36 (MQ=255)
       aaaTTCGAGTCCTGGAGAAAACCGCCGATTccgccc                              >  1:4284667/1‑36 (MQ=255)
         aTTCGAGTCCTGGAGAAAACCGCCGATTccgcccgc                            >  1:5804611/1‑36 (MQ=255)
          ttCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCa                           >  1:721938/1‑36 (MQ=255)
           tCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCAt                          >  1:3906038/1‑36 (MQ=255)
           tCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCAt                          >  1:3674788/1‑36 (MQ=255)
           tCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCAt                          >  1:7461987/1‑36 (MQ=255)
           tCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCAt                          >  1:2908688/1‑36 (MQ=255)
               gTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCAt                      >  1:3565403/1‑36 (MQ=255)
                cacaaaaaaaccccgccgaTTCCGCCCGCATTCATg                     <  1:3935195‑M2/17‑1 (MQ=255)
                 ccTGGAGAAAACCGCCGATTCCGCCCGCATTCATGt                    >  1:2836072/1‑36 (MQ=255)
                   tGGAGAAAACCGCCGCTTCCGCCCTCCTTCATGTTc                  >  1:2976535/1‑36 (MQ=25)
                   tGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTc                  >  1:2472738/1‑36 (MQ=255)
                   gaaaaaaccccgccgaTTCCGCCCGCATTCATGTTc                  <  1:7754492‑M2/20‑1 (MQ=255)
                   aaaaaaaccccgccgaTTCCGCCCGCATTCATGTTc                  <  1:2993200‑M2/20‑1 (MQ=255)
                     aaaaaccccgccgaTTCCGCCCGCATTCATGTTCGt                <  1:1354095‑M2/22‑1 (MQ=255)
                     aaaaaacccgccgaTTCCGCCCGCATTCATGTTCGt                <  1:6360370‑M2/22‑1 (MQ=255)
                       gaAAACCGCCGATTCCGCCCGCATTCATGTTCGttt              >  1:7594650/1‑36 (MQ=255)
                       gaAAACCGCCGATTCCGCCCGCATTCATGTTCGttt              >  1:1456442/1‑36 (MQ=255)
                       aaaccccgccgaTTCCGCCCGCATTCATGTTCGttt              <  1:1020722‑M2/24‑1 (MQ=255)
                       aaaacccgccgaTTCCGCCCGCATTCATGTTCGttt              <  1:153809‑M2/24‑1 (MQ=255)
                        aaaaCCGCCGATTCCGCCCGCATGCATGTTCGtttg             >  1:807768/1‑35 (MQ=255)
                          ccccgccgaTTCCGCCCGCATTCATGTTCGTTTTAc           <  1:5302233‑M2/27‑1 (MQ=255)
                           aCCGCCGATTCCGCCCGCATTCATGTTCGTTTTACg          >  1:3388899/1‑36 (MQ=255)
                           aCCGCCGATTCCGCCCGCATTCATGTTCGTTTTACg          >  1:7018354/1‑36 (MQ=255)
                           aCCGCCGATTCCGCCCGCATTCATGTTCGTTTTACg          >  1:5875654/1‑36 (MQ=255)
                            ccgccgaTTCCGCCCGCATTCATGTTCGTTTTACGc         <  1:6402519‑M2/29‑1 (MQ=255)
                             agccgaTTCCGCCCGCATTCATGTTCGTTTTACGCt        >  1:5504185‑M2/7‑36 (MQ=255)
                              gccgaTTCCGCCCGCATTCATGTTCGTTTTACGCTg       >  1:3492753‑M2/6‑36 (MQ=255)
                              gccgaTTCCGCCCGCATTCATGTTCGTTTTACGCTg       >  1:2562114‑M2/6‑36 (MQ=255)
                               ccgaTTCCGCCCGCATTCATGTTCGTTTTACGCTGa      <  1:3447958‑M2/32‑1 (MQ=255)
                               ccgaTTCCGCCCGCATTCATGTTCGTTTTACGCTGa      <  1:3125037‑M2/32‑1 (MQ=255)
                               ccgaTTCCGCCCGCATTCATGTTCGTTTTACGCTGa      >  1:7946372‑M2/5‑36 (MQ=255)
                                cgaTTCCGCCCGCATTCATGTTCGTTTTACGCTGAc     >  1:2044552‑M2/4‑36 (MQ=255)
                                 gaTTCCGCCCGCATTCATGTTCGTTTTACGCTGAcc    >  1:768829‑M2/3‑36 (MQ=255)
                                 gaTTCCGCCCGCATTCATGTTCGTTTTACGCTGAcc    <  1:5892520‑M2/34‑1 (MQ=255)
                                  aTTCCGCCCGCATTCATGTTCGTTTTACGCTGACCa   <  1:7079218‑M2/35‑1 (MQ=255)
                                  aTTCCGCCCGCATTCATGTTCGTTTTACGCTGACCa   >  1:6285585‑M2/2‑36 (MQ=255)
                                   ttCCGCCCGTATTCATGTTCGTTTTACGCTGACCaa  <  1:2650416/36‑1 (MQ=255)
                                   ttCCGCCCGCATTCATGTTCGTTTTACGCTGACCaa  >  1:2280167/1‑36 (MQ=255)
                                   |                                   
GTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAA  >  REL606/3092542‑3092612

Base quality scores:  ATCG  < 3 ≤  ATCG  < 13 ≤  ATCG  < 25 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG