New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 5099995 =22 (0.760)4 (0.170) 4/214 2.9 19.6% intergenic (+33/+120) PP_4488/phhR hypothetical protein/sigma‑54 dependent transcriptional regulator
?NC_002947 5100077 = 15 (0.640)intergenic (+115/+38) PP_4488/phhR hypothetical protein/sigma‑54 dependent transcriptional regulator
Rejected: Frequency below/above cutoff threshold.

ACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAA  >  NC_002947/5099955‑5100226
                                                                                                                          |                                                                                                                                                     
aCGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAAAGCACAAACATTTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCtt                                                                                                                            >  2:206486/1‑150 (MQ=255)
        tACCTTGGAGTCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACcaac                                                                                                                    <  1:510424/150‑1 (MQ=255)
                        gCCCTGCGGGCCCCTATCGCGACACAAGGCCGGTCCCACAGATAATGCGCTCAGCCAGAGCACCGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACt                                                                                                    >  1:213557/1‑150 (MQ=255)
                         cccTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGTACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACtt                                                                                                   <  2:293842/150‑1 (MQ=255)
                                          GCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATccc                                                                                  >  1:584277‑M1/1‑150 (MQ=255)
                                              CACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACtt                                                                                                   <  2:9852‑M1/129‑1 (MQ=255)
                                              CACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACtt                                                                                                   >  1:9852‑M1/1‑129 (MQ=255)
                                                      gatctaCACAGACAACGCGCTCAGCCATATCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTgg                                                                      >  2:226240‑M1/7‑150 (MQ=255)
                                                        ctcccacaggtaatacgctgtgctctggctgagcgcattatctgtgggagcggccttgtgtcgcgaAAGGACTGCAACTCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTc                                                                                <  2:99041‑M2/72‑1 (MQ=255)
                                                        ctcccacaggtaatacgctgtgctctggctgagcgcattatctgtgggagcggccttgtgtcgcgaAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTc                                                                                >  1:99041‑M2/67‑138 (MQ=255)
                                                           ccacaggtaatacgctgtgctctggctgagcgcattatctgtgggagcggccttgtgtcgcgaAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAAc                                                                             <  1:270686‑M2/75‑1 (MQ=255)
                                                           ccacaggtaatacgctgtgctctggctgagcgcattatctgtgggagcggccttgtgtcgcgaAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAAc                                                                             >  2:270686‑M2/64‑138 (MQ=255)
                                                           ccACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTATTTGACTACCACCTACACCGCCGCCATCACAACCTATTAATCCCGCAACTTGTTGGCAAtt                                                                 >  2:7129/1‑150 (MQ=255)
                                                           CCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAAtt                                                                 <  2:285473‑M1/150‑1 (MQ=255)
                                                            CACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTtggcagatc                                                                <  1:226240‑M1/150‑5 (MQ=255)
                                                             ACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGCCCTTGTGTCGCCAAAGGACTGCAACGCAGTCACAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCaatttt                                                               <  1:111752‑M1/150‑3 (MQ=255)
                                                              CAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTc                                                              <  2:161192‑M1/150‑1 (MQ=255)
                                                                   aaTGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGTAACGCAGTCCCACGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGtgtg                                                         <  1:271097/150‑1 (MQ=255)
                                                                        gcTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCAAGGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGaaa                                                    <  1:7129/150‑1 (MQ=255)
                                                                          tCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAaca                                                  >  1:324798/1‑150 (MQ=255)
                                                                            aGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACAcc                                                <  1:285473/150‑1 (MQ=255)
                                                                              ccAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGt                                                                                                                                       <  1:346687/61‑1 (MQ=255)
                                                                              ccAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGt                                                                                                                                       >  2:346687/1‑61 (MQ=255)
                                                                              ccAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCa                                              <  1:161192/150‑1 (MQ=255)
                                                                                       aGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGc                                     >  1:127207/1‑150 (MQ=255)
                                                                                        gCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGcc                                    >  1:174006/1‑150 (MQ=255)
                                                                                                ccTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGccc                            <  2:213557/150‑1 (MQ=255)
                                                                                                   gtgggagcggccttgtgtcgcgaAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGac                         >  2:379383‑M2/24‑150 (MQ=255)
                                                                                                        agcggccttgtgtcgcgaAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTg                    >  2:677546‑M2/19‑150 (MQ=255)
                                                                                                          cggccttgtgtcgcgaAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTggg                  <  1:379383‑M2/134‑1 (MQ=255)
                                                                                                            gccttgtgtcgcgaAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGaa                <  2:127207‑M2/136‑1 (MQ=255)
                                                                                                                         aAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTc                                                              >  2:463936‑M2/2‑91 (MQ=255)
                                                                                                                         aAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTc                                                              <  1:463936‑M2/90‑1 (MQ=255)
                                                                                                                          aaGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCaaa  <  2:324798/150‑1 (MQ=255)
                                                                                                                          |                                                                                                                                                     
ACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAA  >  NC_002947/5099955‑5100226

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