New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 4649525 =12 (0.720)6 (0.490) 5/208 0.9 50.6% intergenic (+10/+121) fusB/PP_4112 elongation factor G 2/sulfate transporter
?NC_002947 4649633 = 3 (0.250)intergenic (+118/+13) fusB/PP_4112 elongation factor G 2/sulfate transporter

AACCCGCTCCTACAGGATGGCCGCAGGCCCAAACATTGTGCGCCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGACCGCAACAGTTTTCAATGCTGCTCGCCCGCCCGCTTCAACAGCTTCTTGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAA  >  NC_002947/4649548‑4649780
                                                                                     |                                                                                                                                                   
aacccgctcctacaggggggcgcacaatgtttgggcctgcggccatcctgtaggagcgggtttacccgcgaagaggccggcacagACCGCAACAGTTTTCAATGCTGCTCGCCCGCCCGCTTCAACAGCTTCTTGCACCGCTCCGACaaa                                                                                     >  1:366597‑M2/86‑150 (MQ=255)
   ccgctcctacaggggggcgcacaatgtttgggcctgcggccatcctgtaggagcgggtttacccgcgaagaggccggcacagACCGCAACAGtttt                                                                                                                                        <  2:325746‑M2/14‑1 (MQ=255)
   ccgctcctacaggggggcgcacaatgtttgggcctgcggccatcctgtaggagcgggtttacccgcgaagaggccggcacagACCGCAACAGtttt                                                                                                                                        >  1:325746‑M2/83‑96 (MQ=255)
         CTACAGGATGGCCGCAGGCCCAAACATTGTGCGCCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGACCGCAACAGTTTTCAATGCTGCTCGCCCGCCCGCTTCAACAGCTTCTTGCACCGCTCCGACAAATGCACCAcg                                                                            <  2:297511‑M1/150‑1 (MQ=255)
               gggggcgcacaatgtttgggcctgcggccatcctgtaggagcgggtttacccgcgaagaggccggcacagACCGCAACAGTTTTCAATGCTGCTCGCCCGCCCGCTTCAACAGCTTCTTGCACCGCTCCGACAAATGCACCACGCGCaaa                                                                      >  2:147932‑M2/71‑150 (MQ=255)
               gggggcgcacaatgtttgggcctgcggccatcctgtaggagcgggtttacccgcgaagaggccggcacagACCGCAACAGTTTTCAATGCTGCTCGCCCGCCCGCTTCAACAGCTTCTTGCACCGCTCCGACAAATGCACCACGCGCaaa                                                                      >  2:147942‑M2/71‑150 (MQ=255)
                ggggcgcacaatgtttgggcctgcggccatcctgtaggagcgggtttacccgcgaagaggccggcacagACCGCAACAGTTTTCAATGCTGCTCGCCCGCCCGCTTCAACAGCTTCTTGCACCGCTCCGACAAATGCACCACGCGCAAAt                                                                     <  2:366597‑M2/81‑1 (MQ=255)
                               aaCATTGTGCGCCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGACCGCAACAGTTTTCAATGCTGCTCGCCCGCCCGCTTCAACAGCTTCTTGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTgg                                                      >  2:175167/1‑150 (MQ=255)
                                   ttGTGCGCCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGACCGCAACAGTTTTCAATGCTGCTCGCCCGCCCGCTTCAACAGCTTCTTGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATa                                                  >  1:343919/1‑150 (MQ=255)
                                   ttGTGCGCCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGACCGCAACAGTTTTCAATGCTGCTCGCCCGCCCGCTTCAACAGCTTCTTGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATa                                                  >  1:75270/1‑150 (MQ=255)
                                                tgtaggagcgggtttacccgcgaagaggccggcacagACCGCAACAGTTTTCAATGCTGCTCGCCCGCCCGCTTCAACAGCTTCTTGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAgg                                     <  1:175167‑M2/113‑1 (MQ=255)
                                                       gcgggattccccgcgaagaggccggcacagACCGCAACAGTTTTCAATGCTGCTCGCCCGCCCGCTTCAACAGCTTCTTGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCa                              <  2:288311‑M2/120‑1 (MQ=255)
                                                          ggtttacccgcgaagaggccggcacagACCGCAACAGTTTTCAATGCTGCTCGCCCGCCCGCTTCAACAGCTTCTTGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCTCGCAGGGTTTTCAgcg                           >  2:157998‑M2/28‑150 (MQ=255)
                                                          ggtttacccgcgaagaggccggcacagACCGCAACAGTTTTCAATGCTGCTCGCCCGCCCGCTTCAACAGCTTCTTGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCTCGCAGGGTTTTCAgcg                           >  2:351032‑M2/28‑150 (MQ=255)
                                                           gtttacccgcgaagaggccggcacagACCGCAACAGTTTTCAATGCTGCTCGCCCGCCCGCTTCAACAGCTTCTTGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAgcgc                          >  1:262120‑M2/27‑150 (MQ=255)
                                                           gtttacccgcgaagaggccggcacagACCGCAACAGTTTTCAATGCTGCTCGCCCGCCCGCTTCAACAGCTTCTTGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGCCATAGCGCTCGCGCAGGGTTTTCAgcgc                          >  2:68835‑M2/27‑150 (MQ=255)
                                                                           gccggcacagACCGCAACAGTTTTCAATGCTGCTCGCCCGCCCGCTTCAACAGCTTCTTGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGaataa          >  1:105468‑M2/11‑150 (MQ=255)
                                                                                   agACCGCAACAGTTTTCAATGCTGCTCGCCCGCCCGCTTCAACAGCTTCTTGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACaaa  >  1:89245‑M2/3‑150 (MQ=255)
                                                                                     |                                                                                                                                                   
AACCCGCTCCTACAGGATGGCCGCAGGCCCAAACATTGTGCGCCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGACCGCAACAGTTTTCAATGCTGCTCGCCCGCCCGCTTCAACAGCTTCTTGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAA  >  NC_002947/4649548‑4649780

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