New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 394682411 (0.190)8 (0.140) 7/274 6.4 18.9% intergenic (+5/‑48) rrfC/aspT 5S ribosomal RNA of rrnC operon/tRNA‑Asp
?NC_000913 = 4040681 58 (1.020)intergenic (+45/+225) rrfA/mobB 5S ribosomal RNA of rrnA operon/molybdopterin‑guanine dinucleotide biosynthesis protein B
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAAT  >  NC_000913/3946679‑3946972
                                                                                                                                                 |                                                                                                                                                    
gCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTaa                                                                                                                                                   >  2:792927/1‑149 (MQ=32)
  ggTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAGc                                                                                                                                                 <  1:283778/149‑1 (MQ=32)
          gaAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAAGTCAGAAGTGAACCGCCGTGGCGCCGATGGTAGTGTGGGGTCTCCCAATGCGAGAGTAGGGAACTGCCAGGCATCACATTAAGCAGTAAGcc                                                                                                                                         >  1:437017/1‑149 (MQ=11)
          gaAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCAGGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAGCAGTAAGcc                                                                                                                                         >  2:894938/1‑149 (MQ=11)
                                                                                                                              ggAACTGCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTa                     >  1:714964/1‑149 (MQ=32)
                                                                                                                                   tGCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAg                <  2:948321/149‑1 (MQ=37)
                                                                                                                                           aTCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCGGCCGCGCTAATGTGGGGCGTAg        >  1:345927/1‑149 (MQ=11)
                                                                                                                                           aTCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAg        >  1:382310/1‑149 (MQ=255)
                                                                                                                                            tCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGt       >  1:339148/1‑149 (MQ=255)
                                                                                                                                              aaaTTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTCGGGGCGTAGTTc     >  1:294832/1‑149 (MQ=255)
                                                                                                                                                 ttAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAAt  <  1:792927/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAAT  >  NC_000913/3946679‑3946972

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