New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? minE 2533382 =76 (1.470)3 (0.060) 3/96 NT 4.3% coding (6/1035 nt) rfaQ lipopolysaccharide core biosynthesis protein
?minE 2533402 = 63 (1.290)coding (26/1035 nt) rfaQ lipopolysaccharide core biosynthesis protein
Rejected: Frequency below/above cutoff threshold.

ATTTTTGCATCAGGATAATTCTGCTTGAGCGTACTGATGACAGGAGTAGTTAATAACATATCCCCATGATAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  minE/2533453‑2533382
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tatTAACTACTCCTGTCATCAGTACGCTCAAGCAGAATTATCCTGATGCAAAAATCGATATGCTGCTTTAT  >  minE/2533402‑2533469
                                                                                                                                            
ATTTTTGCATCAGGATAATTCTGCTTGAGCGTACTGATGACAGGAGTAGTTAATAACATATCCCCATGATA                                                                       <  1:1235935/71‑1
ATTTTTGCATCAGGATAATTCTGCTTGAGCGTACTGATGACAGGAGTAGTTAATAACATATCCCCATGATA                                                                       <  1:1427788/71‑1
 TTTTTGCATCAGGATAATTCTGCTTGAGCGTACTGATGACAGGAGTAGTTAATAACATATCCCCATGATAT                                                                      >  1:77219/1‑71
    TTGCATCAGGATAATTCTGCTTGAGCGTACTGATGACAGGAGTAGTTAATAACATATCCCCATGATAT                                                                      <  1:3000763/68‑1
                   TCTGCTTGAGCGTACTGATGACAGGAGTAGTTAATAACATATCCCCATGATAT                                                                      >  1:1387777/1‑53
                                                 TTAATAACATATCCCCATGATATTAACTACTCCTGTCATCAGTACGCTCAAGCAGAATTATCC                              <  1:2957863/63‑1
                                                 TTAATAACATATCCCCATGATATTAACTACTCCTGTCATCAGTACGCTCAAGCAGAATTA                                 <  1:3527541/60‑1
                                                          TATCCCCATGATATTAACTACTCCTGTCATCAGTACGCTCAAGCAGAATTATCCTGATGCAAAAATCGATA             >  1:2249011/1‑71
                                                                     TATTAACTACTCCTGTCATCAGTACGCTCAAGCAGAATTATC                               >  1:1127937/1‑42
                                                                     TATTAACTACTCCTGTCATCAGTACGCTCAAGCAGAATTATC                               >  1:1566592/1‑42
                                                                     TATTAACTACTCCTGTCATCAGTACGCTCAAGCAGAATTATC                               >  1:884507/1‑42
                                                                      ATTAACTACTCCTGTCATCAGTACGCTCAAGCAGAATTATCCTGATGCAAAAATCGATATGCTGCTTTAT  >  1:2091077/1‑70
                                                                      ATTAACTACTCCTGTCATCAGTACGCTCAAGCAGAATTATCCTGATGCAAAAATCGATATGCTGCTTTAT  <  1:3134449/70‑1
                                                                      ATTAACTACTCCTGTCATCAGTACGCTCAAGCAGAATTATCCTGATGCAAAAATCGATATGCTGCTTTAT  <  1:426508/70‑1
                                                                      ATTAACTACTCCTGTCATCAGTACGCTCAAGCAGAATTATCCTGAT                          >  1:9382/1‑46
                                                                       TTAACTACTCCTGTCATCAGTACGCTCAAGCAGAATTATCCTGATGCAAAAATCGATATGCTGCTTTAT  >  1:164161/1‑69
                                                                                                                                            
ATTTTTGCATCAGGATAATTCTGCTTGAGCGTACTGATGACAGGAGTAGTTAATAACATATCCCCATGATAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  minE/2533453‑2533382
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tatTAACTACTCCTGTCATCAGTACGCTCAAGCAGAATTATCCTGATGCAAAAATCGATATGCTGCTTTAT  >  minE/2533402‑2533469

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