New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AE015451 4960421 =72 (1.050)3 (0.050) 3/326 NT 4.1% coding (50/1779 nt) fliF Flagellar M‑ring protein
?AE015451 4960435 = 72 (1.090)coding (36/1779 nt) fliF Flagellar M‑ring protein
Rejected: Frequency below/above cutoff threshold.

AGCAAGCCTCGACCCAACTGGCCAACGCGTTCGAAATCGGCAAGAGCGGCGTCGACCTGACTGACGTGATGATCGCTTCGCAAAAGGCCAGCGTGTCGATGCAAGCCATGACCCAGGTACGCAACAAGCTGGTCCAGGCGTACCAGGACATCATGCAGATGCCGGTTTGAGGACGGATGTAACCCATGGCTGAAGCAGTCGTCGAAAACGCCCCCGCCAAAGCTGGCTCGCCAGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/4960655‑4960421
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gccagcTTTGGCGGGGGCGTTTTCGACGACTGCTTCAGCCATGGGTTACATCCGTCCTCAAACCGGCATCTGCATGAT  >  AE015451/4960435‑4960506
                                                                                                                                                                                                                                                                                                                   
AGCAAGCCTCGACCCAACTGGCCAACGCGTTCGAAATCGGCAAGAGCGGCGTCGACCTGACTGACGTGATGATCGCTTCGCAAAAGGCCAGCGTGTCGATGCAAGCCATGACCCAGGTACGCAACAAGCTGGTCCAGGCGTACCAGGACATCATGCAGATGCCGGTTTGAGGACGGATGTAACCCATGGCTGAAGCAGTCGTCGAAAACGCCCCCGCCAAAGCTGGCTCGCCAGC                                                                          >  2:611595/1‑235
                                                                                                                                       AGGCGTACCAGGACATCATGCAGATGCCGGTTTGAGGACGGATGTAACCCATGGCTGAAGCAGTCGTCGAAAACGCCCCCGCCAAAGCTGGCTCGCCAGCTTTGGCGGGGGCGTTTTCGACGACTGCTTCAGCCATGGGTT                                 >  1:719014/1‑141
                                                                                                                                       AGGCGTACCAGGACATCATGCAGATGCCGGTTTGAGGACGGATGTAACCCATGGCTGAAGCAGTCGTCGAAAACGCCCCCGCCAAAGCTGGCTCGCCAGCTTTGGCGGGGGCGTTTTCGACGACTGCTTCAGCCATGGGTT                                 <  2:719011/141‑1
                                                                                                                                                                                                    AGTCGTAGAAAACGCCCCCGCCAAAGCTGGCTCGCCAGCTTTGGCGGGGGCGTTTTCGACGAATGCTTCAGCCATGGGGTACATCCGTCCTCAAACCCGCATCTGCGTGgt  >  2:541112/1‑109
                                                                                                                                                                                                                                                                                                                   
AGCAAGCCTCGACCCAACTGGCCAACGCGTTCGAAATCGGCAAGAGCGGCGTCGACCTGACTGACGTGATGATCGCTTCGCAAAAGGCCAGCGTGTCGATGCAAGCCATGACCCAGGTACGCAACAAGCTGGTCCAGGCGTACCAGGACATCATGCAGATGCCGGTTTGAGGACGGATGTAACCCATGGCTGAAGCAGTCGTCGAAAACGCCCCCGCCAAAGCTGGCTCGCCAGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/4960655‑4960421
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gccagcTTTGGCGGGGGCGTTTTCGACGACTGCTTCAGCCATGGGTTACATCCGTCCTCAAACCGGCATCTGCATGAT  >  AE015451/4960435‑4960506

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