New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? minE = 141782120 (0.330)21 (0.370)
+TCAG
11/90 NT 38.9% intergenic (+102/‑231) proL/narP tRNA‑Pro/DNA‑binding response regulator in two‑component regulatory system with NarQ or NarX
?minE 2713590 = 52 (0.850)coding (771/2544 nt) nirB nitrite reductase, large subunit, NAD(P)H‑binding

TTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  minE/1417756‑1417821
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑AACGTCCAGACCACACTG  >  minE/2713590‑2713607
                                                                  ||||                  
TTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAG                    <  1:4676669/70‑1
TTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAG                    >  1:2496168/1‑70
TTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAG                    <  1:4142999/70‑1
      TATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTC              <  1:5638077/70‑1
       ATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTC              <  1:4456886/69‑1
          TGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGAC         >  1:2516611/1‑71
                 AATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACc   >  1:2040293/1‑69
                 AATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACc   >  1:2147127/1‑69
                  ATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACcc  <  1:3720370/70‑3
                       ATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCAC      <  1:2748099/61‑1
                          AACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAG                    <  1:3588710/44‑1
                                 CGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGAC         <  1:2816888/48‑1
                                 CGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGAC         <  1:5693106/48‑1
                                 CGTGACACATCACAGCCTGTTTATTTTCTGTTATCAG                    <  1:5686206/37‑1
                                   TGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCA       >  1:1944238/1‑48
                                       ACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCAC      >  1:4335409/1‑45
                                         ATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACc   >  1:177685/1‑45
                                         ATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACc   >  1:4603690/1‑45
                                         ATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGA          <  1:876423/39‑1
                                          TCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACA     <  1:3358197/43‑1
                                          TCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACA     <  1:3159619/43‑1
                                           CACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACAC    >  1:4676462/1‑43
                                            ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACc   >  1:5361505/1‑42
                                            ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACA     >  1:141749/1‑41
                                            ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGA          >  1:4684939/1‑36
                                            ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGA          >  1:5485718/1‑36
                                                                  ||||                  
TTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  minE/1417756‑1417821
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑AACGTCCAGACCACACTG  >  minE/2713590‑2713607

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.