New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3423374 =85 (1.280)8 (0.120) 6/280 8.5 14.7% intergenic (+182/+47) yhdZ/rrfF putative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
?NC_000913 = 3946821 9 (0.140)intergenic (+3/‑50) rrfC/aspT 5S ribosomal RNA of rrnC operon/tRNA‑Asp
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGG  >  NC_000913/3946682‑3946960
                                                                                                                                           |                                                                                                                                           
tCTGATGAAACAGACTTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAAtt                                                                                                                                          <  2:642322/143‑1 (MQ=17)
   gatgaGACAGAATTTGGGTGGCGGCCGGAGCGCGGTGGTGCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAg                                                                                                                                       <  1:862212/143‑1 (MQ=11)
   cacgaAACAGAATCTGCCAGGCGGCCGGAGCGCGGTGGTCCCACCCGACCCCATGCCGAACTGAGCTGTCTTACTGCGTAGCGCCGAGGGTAGTGTGTGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAg                                                                                                                                       <  1:409381/140‑1 (MQ=11)
      gaAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCAGGCCGAACTCAGAAGTGAAAGGCCGTAGCGCCGCTGGTAGTGTGGGGTCTCCCCATGCGAGAGGAAGGAACTGCCAGGCACCAGATTAAGCAg                                                                                                                                    >  1:537543/1‑143 (MQ=11)
                                                                                                                                 ccAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTaa         >  2:427270/1‑143 (MQ=33)
                                                                                                                                       aTCAAACTAAGCAGTAAGCCGGTCACAAGACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCGGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAggg   <  2:237663/143‑1 (MQ=255)
                                                                                                                                        tCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCGCCCTAATTAgggg  >  2:451314/1‑143 (MQ=38)
                                                                                                                                        tCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAgggg  >  2:673832/1‑143 (MQ=39)
                                                                                                                                        tCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTAGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGGGGGTTCGAGTGCCGTCCGTTCCGCACGCCTAAGTAgcgg  >  1:790097/1‑143 (MQ=1)
                                                                                                                                           |                                                                                                                                           
TCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGG  >  NC_000913/3946682‑3946960

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