New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 592057NA (NA)5 (0.100) 5/60 NT 10.9% noncoding (1343/1343 nt) IS186 repeat region
?REL606 1391825 = 41 (0.680)coding (230/1614 nt) mppA murein tripeptide (L‑ala‑gamma‑D‑glutamyl‑meso‑DAP) transporter subunit

GTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGAGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTA  >  REL606/1391790‑1391860
                                   |                                   
tcagatggttcaacccttaagttagcgcttatggga                                     >  1:7450619‑M2/36‑36 (MQ=255)
tcagatggttcaacccttaagttagcgcttatggga                                     <  1:1632272‑M2/1‑1 (MQ=255)
tcagatggttcaacccttaagttagcgcttatggga                                     >  1:5361957‑M2/36‑36 (MQ=255)
tcagatggttcaacccttaagttagcgcttatggga                                     >  1:200019‑M2/36‑36 (MQ=255)
tcagatggttcaacccttaagttagcgcttatggga                                     >  1:2500076‑M2/36‑36 (MQ=255)
tcagatggttcaacccttaagtaagcgcttatggga                                     <  1:793746‑M2/1‑1 (MQ=255)
gTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGga                                     >  1:5669571/1‑36 (MQ=255)
gTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGga                                     >  1:1918871/1‑36 (MQ=255)
gTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGga                                     <  1:4612823/36‑1 (MQ=255)
 cagatggttcaacccttaagttagcgcttatgggag                                    >  1:6204150‑M2/35‑36 (MQ=255)
 cagatggttcaacccttaagttagcgcttatgggag                                    <  1:5135958‑M2/2‑1 (MQ=255)
    gAAGGTCTGGTGAATCAGAACGAAAAAGGGGAGAtt                                 >  1:1607571/1‑36 (MQ=255)
    gAAGGTCTGATGAATCAGAACGAAAAAGGGGAGAtt                                 <  1:4988533/36‑1 (MQ=255)
     aaGGTCTGGTGAATCAGAACGAAAAAGGGGAGATTg                                >  1:831473/1‑36 (MQ=255)
      aGGTCTGGTGAATCAGAACGAAAAAGGGGAGATTGt                               >  1:2003296/1‑36 (MQ=255)
       ggTCTGGTGAATCAGAACGAAAAAGGGGAGATTGTc                              >  1:7424469/1‑36 (MQ=255)
       ggTCTGGTGAATCAGAACGAAAAAGGGGAGATTGTc                              <  1:4595596/36‑1 (MQ=255)
        gTCTGGTGAATCAGAACGAAAAAGGGGAGATTGTcc                             <  1:1624369/36‑1 (MQ=255)
        gTCTGGTGAATCAGAACGAAAAAGGGGAGATTGTcc                             >  1:5429642/1‑36 (MQ=255)
        gTCTGGTGAATCAGAACGAAAAAGGGGAGATTGTcc                             <  1:2163520/36‑1 (MQ=255)
          caacccttaagttagcgcttatgggAGATTGTcccc                           <  1:6483187‑M2/11‑1 (MQ=255)
          cTGGTGAATCAGAACGAAAAAGGGGAGATTGTcccc                           <  1:2253245/36‑1 (MQ=255)
          cTGGTGAATCAGAACGAAAAAGGGGAGATTGTcccc                           <  1:3351789/36‑1 (MQ=255)
          cTGGTGAATCAGAACGAAAAAGGGGAGATTGTcccc                           >  1:5301279/1‑36 (MQ=255)
           tGGTGAATCAGAACGAAAAAGGGGAGATTGTCCCCg                          >  1:7297699/1‑36 (MQ=255)
            ggTGAATCAGAACGAAAAAGGGGAGATTGTCCCCgg                         >  1:7276161/1‑36 (MQ=255)
            ggTGAATCAGAACGAAAAAGGGGAGATTGTCCCCgg                         <  1:1709984/36‑1 (MQ=255)
            ggTGAATCAGAACGAAAAAGGGGAGATTGTCCCCgg                         >  1:5097148/1‑36 (MQ=255)
              tGAATCAGAACGAAAAAGGGGAGATTGTCACCGGCg                       >  1:6123486/1‑36 (MQ=255)
              ccttaagttagcgcttatgggAGATTGTCCCCGGCg                       <  1:2515403‑M2/15‑1 (MQ=255)
               gAATCAGAACGAAAAAGGGTAGATTGTCCCCGGCGt                      <  1:5637328/36‑1 (MQ=255)
               gAATCAGAACGAAAAAGGGGAGATTGTCCCCGGCGt                      >  1:7419090/1‑36 (MQ=255)
                aaTCAGAACGAAAAAGGGGAGATTGTCCCCTGCGtt                     >  1:6545095/1‑36 (MQ=255)
                 ataaGAACGAAAAAGGGTAGATTGTCCCCGGCGTTg                    <  1:2213168/33‑1 (MQ=37)
                 aTCAGAACGAAAAAGGGGAGATTGTCCCCGGCGTTg                    <  1:4225131/36‑1 (MQ=255)
                    aGAACGAAAAAGGGGAGATTGTCCCCGGCGTCGCGa                 <  1:1020324/36‑1 (MQ=255)
                     gAACGAAAAAGGGGAGATTGTCCCCGGCGTTGCGAc                <  1:1219131/36‑1 (MQ=255)
                     gAACGAAAAAGGGGAGATTGTCCCCGGCGTTGCGAc                <  1:2238156/36‑1 (MQ=255)
                      aaCGAAAAAGGGGAGATTGCCCCTGGGGTTGCGACt               >  1:7026986/1‑36 (MQ=25)
                       aCGAAAAAGGGGAGATTGTCCCCGGCGTTGCGACTc              >  1:3469/1‑36 (MQ=255)
                        cGAAAAAGGGGAGATTGTCCCCGGCGTTGCGACTCa             >  1:5208951/1‑36 (MQ=255)
                        cGAAAAAGGGGAGATTGTCCCCGGCGTTGCGACTCa             >  1:3394575/1‑36 (MQ=255)
                         gAAAAAGGGGAGATTGTCCCCGGCGTTGCGACTCAg            >  1:5496398/1‑36 (MQ=255)
                          aaaaaGGGGAGATTGTCCCCGGCGTTGCGACTCAGt           <  1:827539/36‑1 (MQ=255)
                          aaaaaGGGGAGATTGTCCCCGGCGTTGCGACTCAGt           >  1:4741932/1‑36 (MQ=255)
                          aaaaaGGGGAGATTGTCCCCGGAGTTGCGACTCAgg           >  1:1281159/1‑35 (MQ=255)
                           aaaaGGGGAGATTGTCCCCGGCGTTGCGACTCAGTg          >  1:6629210/1‑36 (MQ=255)
                            aaaGGGGAGATTGTCCCCGGCGTTGCGACTCAGTgg         >  1:3180104/1‑36 (MQ=255)
                             tatgggAGATTGTCCCCGGCGTTGCGACTCAGTGGa        <  1:6746571‑M2/30‑1 (MQ=255)
                             aaGGGGAGATTGTCCCCGGCGTTGCGACTCAGTGGa        >  1:318550/1‑36 (MQ=255)
                             aaGGGGAGATTGTCCCCGGCGTTGCGACTCAGTGGa        >  1:5720529/1‑36 (MQ=255)
                              atgggAGATTGTCCCCGGCGTTGCGACTCAGTGGaa       <  1:53044‑M2/31‑1 (MQ=255)
                               tgggAGATTGTCCCCGGCGTTGCGACTCAGTGGaaa      <  1:5299704‑M2/32‑1 (MQ=255)
                               ggggAGATTGTCCCCGGCGTTGCGACTCAGTGGaaa      >  1:1694201/1‑36 (MQ=255)
                               ggggAGATTGTCCCCGGCGTTGCGACTCAGTGGaaa      >  1:7358787/1‑36 (MQ=255)
                               ggggAGATTGTCCCCGGCGTTGCGACTCAGTGGaaa      <  1:483447/36‑1 (MQ=255)
                               cgggAGATTGTCCCCTGCGTTGCGACTCATTGGaaa      <  1:2343842/35‑1 (MQ=255)
                                 ggAGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAg    >  1:6246133‑M2/3‑36 (MQ=255)
                                 ggAGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAg    >  1:6156147‑M2/3‑36 (MQ=255)
                                  gAGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGt   >  1:3673900‑M2/2‑36 (MQ=255)
                                   AGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTa  >  1:2865751‑M2/1‑36 (MQ=255)
                                   |                                   
GTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGAGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTA  >  REL606/1391790‑1391860

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