New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? minE = 270261771 (1.370)3 (0.060) 3/92 NT 4.3% coding (8/1089 nt) aroB 3‑dehydroquinate synthase
?minE = 2702628 69 (1.480)coding (19/1089 nt) aroB 3‑dehydroquinate synthase
Rejected: Frequency below/above cutoff threshold.

AAGCAACTAATTCTGGCTTTATATACACTCGTCTGCGGGTACAGTAATTAAGGTGGATGTCGCGTTATGGAGAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  minE/2702544‑2702617
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gagagTAACGACAATCCTCTCCATAACGCGACATCCACCTTAATTACTGTACCCGCAGACGAGTGTATATA  <  minE/2702628‑2702563
                                                                                                                                            
AAGCAACTAATTCTGGCTTTATATACACTCGTCTGCGGGTACAGTAATTAAGGTGGATGTCGCGTTATGGA                                                                       <  1:539482/71‑1
AAGCAACTAATTCTGGCTTTATATACACTCGTCTGCGGGTACAGTAATTAAGGTGGATGTCGCGTTATGGA                                                                       <  1:605543/71‑1
   CAACTAATTCTGGCTTTATATACACTCGTCTGCGGGTACAGTAATTAAGGTGGATGTCGCGTTATGGAG                                                                      <  1:3347341/69‑1
     ACTAATTCTGGCTTTATATACACTCGTCTGCGGGTACAGTAATTAAGGTGGATGTCGCGTTATGGAGAG                                                                    >  1:3034233/1‑69
              GGCTTTATATACACTCGTCTGCGGGTACAGTAATTAAGGTGGATGTCGCGTTATGGAGAG                                                                    <  1:2923986/60‑1
                       TACACTCGTCTGCGGGTACAGTAATTAAGGTGGATGTCGCGTTATGGAGAGTAACGACAATCCTCTCCATA                                                <  1:2061132/71‑1
                         CACTCGTCTGCGGGTACAGTAATTAAGGTGGATGTCGCGTTATGGAGAG                                                                    <  1:474062/49‑1
                              GTCTGCGGGTACAGTAATTAAGGTGGATGTCGCGTTATGGAGAGTAACGACAATCCTC                                                      <  1:1060086/58‑1
                                                     TGGATGTCGCGTTATGGAGAGTAACGACAATCCTCT                                                     >  1:1384327/1‑36
                                                                      AGAGTAACGACAATCCTCTCCATAACGCGACATCCACCTTAATTACTGTACCCGCAGACGAGTGTATATA  >  1:1301002/1‑70
                                                                      AGAGTAACGACAATCCTCTCCATAACGCGACATCCACCTTAATTACTGTACCCGCAGACGAGTGTATATA  >  1:489/1‑70
                                                                      AGAGTAACGACAATCCTCTCCATAACGCGACATCCACCTTAATT                            <  1:3556273/44‑1
                                                                         GTAACGACAATCCTCTCCATAACGCGACATCCACCTT                                <  1:2308183/37‑1
                                                                                                                                            
AAGCAACTAATTCTGGCTTTATATACACTCGTCTGCGGGTACAGTAATTAAGGTGGATGTCGCGTTATGGAGAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  minE/2702544‑2702617
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gagagTAACGACAATCCTCTCCATAACGCGACATCCACCTTAATTACTGTACCCGCAGACGAGTGTATATA  <  minE/2702628‑2702563

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ 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.