New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 4652552 =2 (0.070)4 (0.200) 3/216 2.6 51.7% intergenic (+10/+121) fusB/PP_4112 elongation factor G 2/sulfate transporter
?NC_002947 4652660 = 6 (0.300)intergenic (+118/+13) fusB/PP_4112 elongation factor G 2/sulfate transporter

GTGAAATGTTCGGCTACGCCACCTCGATGCGCTCGATGACCCAGGGCCGGGCCAGCTTCTCGATGGAGTTCACCCGCTATGCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCAAACATTGTGCGCCCCCCTGTAGGAGCGGGT  >  NC_002947/4652406‑4652635
                                                                                                                                                  |                                                                                   
tgtcgattattcggccgttgcggcgctgaaaaccctgcgcgagcgctatgccaaggcgggcaagcatttgcgcgtggtgcatttgtcggagcggtgcaagaagctgttgaagcgggcgggcgagcagcattgaaaactgttgcggcCtgt                                                                                  >  1:438341‑M1/147‑150 (MQ=255)
  tcgattattcggccgttgcggcgctgaaaaccctgcgcgagcgctatgccaaggcgggcaagcatttgcgcgtggtgcatttgtcggagcggtgcaagaagctgttgaagcgggcgggcgagcagcattgaaaactgttgcggtCTGTGc                                                                                <  1:311210‑M1/6‑1 (MQ=255)
            ggccgttgcggcgctgaaaaccctgcgcgagcgctatgccaaggcgggcaagcatttgcgcgtggtgcatttgtcggagcggtgcaagaagctgttgaagcgggcgggcgagcagcattgaaaactgttgcggtCTGTGCCGGCCTCTTc                                                                      <  1:465229‑M1/16‑1 (MQ=255)
              ccgttgcggcgctgaaaaccctgcgcgagcgctatgccaaggcgggcaagcatttgcgcgtggtgcatttgtcggagcggtgcaagaagctgttgaagcgggcgggcgagcagcattgaaaactgttgcggtCTGTGCCGGCCTCTTcgc                                                                    >  2:465229‑M1/133‑150 (MQ=255)
                      cTCGATGCGCTCGATGACCCAGGGCCGGGCCAGCTTCTCGATGGAGTTCACCCGCTATGCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAAC                                                            <  1:51848‑M1/150‑1 (MQ=255)
                                                    cAGCTTCTCGATGGTGTTCACCCGCTATGCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCaaa                              <  1:132996/150‑1 (MQ=255)
                                                    aaggcgggcaagcatttgcgcgtggtgcatttgtcggagcggtgcaagaagctgttgaagcgggcgggcgagcagcattgaaaactgttgcggtCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCaaa                              >  1:324413‑M1/95‑150 (MQ=255)
                                                     aGCTTCTCGATGGAGTTCACCCGCTATGCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGggggcgcacaatgtttgg                             <  2:330135‑M1/150‑19 (MQ=255)
                                                          ctcGATGGAGTTCACCCGCTATGCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCAAACATtgt                        <  1:292202/150‑1 (MQ=255)
                                                                       gcgtggtgcatttgtcggagcggtgcaagaagctgttgaagcgggcgggcgagcagcattgaaaactgttgcggtCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCAAACATTGTGCGCCCCCCTGTa           >  1:131132‑M1/76‑150 (MQ=255)
                                                                       gcgtggtgcatttgtcggagcggtgcaagaagctgttgaagcgggcgggcgagcagcattgaaaactgttgcggtCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCAAACATTGTGCGCCCCCCTGTa           >  1:131128‑M1/76‑150 (MQ=255)
                                                                       aCCCGCTATGCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGggggcgcacaatgtttgggcctgcggccatcctgta           <  2:439978‑M1/150‑37 (MQ=255)
                                                                          tggtgcatttgtcggagcggtgcaagaagctgttgaagcgggcgggcgagcagcattgaaaactgttgcggtCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCAAACATTGTGCGCCCCCCTGTAGGa        >  2:339147‑M1/73‑150 (MQ=255)
                                                                          cGCTATGCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGggggcgcacaatgtttgggcctgcggccatcctgtagga        <  2:141048‑M1/150‑40 (MQ=255)
                                                                          cGCTATGCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGggggcgcacaatgtttgggcctgcggccatcctgtagga        <  1:233697‑M1/150‑40 (MQ=255)
                                                                              aTGCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGggggcgcacaatgtttgggcctgcggccatcctgtaggagcgg    <  2:234979‑M1/150‑44 (MQ=255)
                                                                               tGCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGggggcgcacaatgtttgggcctgcggccatcctgtaggagcggg   <  2:54907‑M1/150‑45 (MQ=255)
                                                                                gCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCAAACATTGTGCGCCCCcctgtaggagcgggt  >  1:386541‑M2/1‑135 (MQ=255)
                                                                                gCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCAAACATTGTGCGCCCCcctgtaggagcgggt  >  1:386500‑M2/1‑135 (MQ=255)
                                                                                                                                                  |                                                                                   
GTGAAATGTTCGGCTACGCCACCTCGATGCGCTCGATGACCCAGGGCCGGGCCAGCTTCTCGATGGAGTTCACCCGCTATGCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCAAACATTGTGCGCCCCCCTGTAGGAGCGGGT  >  NC_002947/4652406‑4652635

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.