New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? USA300TCH1516_ALE = 5127670 (0.000)9 (0.100) 6/150 NT 100% intergenic (+835/‑5928) recR/speA_2 Recombination protein RecR/Arginine decarboxylase
?USA300TCH1516_ALE = 512806 NA (NA)intergenic (+874/‑5889) recR/speA_2 Recombination protein RecR/Arginine decarboxylase

AAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATAATTTTTATGGAGAGTTTGATCCTGGCTCAGGATGAACGCTGGCGGCGTGCCTAATACATGCAAGTCGAGCGAACG  >  USA300TCH1516_ALE/512677‑512857
                                                                                          |                                                                                          
aaaTGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGC                                                                                            <  1:2805837‑M1/91‑1 (MQ=255)
 aaTGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCt                                                                                           <  1:874396/91‑1 (MQ=31)
 aaTGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCt                                                                                           <  1:1650470/91‑1 (MQ=31)
                             tGTCAACGTTAATTCCAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCtaatcaaactctccataaaaattatgatgt                                                              <  1:1740696‑M1/91‑31 (MQ=255)
                               tCAACGTTAATTCCAAAAACGTAACTATTAGTTACAAACATTATTTAGTATTTATGAGCtaatcaaactctccataaaaattatgatgttt                                                            <  1:4079583‑M1/91‑33 (MQ=255)
                                cAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCtaatcaaactctccataaaaattatgatgttt                                                            <  1:2365646‑M1/91‑33 (MQ=255)
                                cAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCtaatcaaactctccataaaaattatgatgttt                                                            <  1:2365657‑M1/91‑33 (MQ=255)
                                                    gTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCcaggatcaaactctccataaaaattatgatgtttgattagctcataaatact                                        <  1:3822666‑M1/91‑53 (MQ=255)
                                                     tAACTATAAGTTACAAACATTATTTAGTATTTATGAGCcaggatcaaactctccataaaaattatgatgtttgattagttcataaatacta                                       <  1:3201991‑M1/91‑54 (MQ=255)
                                                      aaCTATAAGTTACAAACATTATTTAGTATTTATGAGCcaggatcaaactctccataaaaattatgatgtttgattagctcataaatactaa                                      <  1:510350‑M1/91‑55 (MQ=255)
                                                       aCTATAAGTTACAAACATTATTTAGTATTTATGAGCcaggatcaaactctccataaaaattatgatgtttgattagctcataaatactaaa                                     <  1:3269232‑M1/91‑56 (MQ=255)
                                                           taAGTTACAAACATTATTTAGTATTTATGAGCcaggatcaaactctccataaaaattatgatgtttgattagctcataaatactaaataat                                 <  1:881150‑M1/91‑60 (MQ=255)
                                                                                      tGAGCcaggatcaaactctccataaaaattatgatgtttgattagctcataaatactaaataatgtttgtaacttatagttacgttttttg      <  1:1011771‑M1/91‑87 (MQ=255)
                                                                                      tGAGCcaggatcaaactctccataaaaattatgatgtttgattagctcataaatactaaataatgtttgtaacttatagttacgttttttg      <  1:1011785‑M1/91‑87 (MQ=255)
                                                                                      tGAGCcaggatcaaactctccataaaaattatgatgtttgattagctcataaatactaaataatgtttgtaacttatagttacgttttttg      <  1:3937188‑M1/91‑87 (MQ=255)
                                                                                       gAGCcaggatcaaactctccataaaaattatgatgtttgattagctcataaatactaaataatgtttgtaacttatagttacgttttttgg     >  1:711879‑M1/1‑4 (MQ=255)
                                                                                       gAGCcaggatcaaactctccataaaaattatgatgtttgattagctcataaatactaaataatgtttgtaacttatagttacgttttttgg     >  1:711757‑M1/1‑4 (MQ=255)
                                                                                       gAGCcaggatcaaactctccataaaaattatgatgtttgattagctcataaatactaaataatgtttgtaacttatagttacgttttttgg     <  1:3136847‑M1/91‑88 (MQ=255)
                                                                                       gAGCcaggatcaaactctccataaaaattatgatgtttgattagctcataaatactaaataatgtttgtaacttatagttacgttttttgg     >  1:1153904‑M1/1‑4 (MQ=255)
                                                                                       gAGCcaggatcaaactctccataaaaattatgatgtttgattagctcataaatactaaataatgtttgtaacttatagttacgttttttgg     <  1:1428574‑M1/91‑88 (MQ=255)
                                                                                       gAGCcaggatcaaactctccataaaaattatgatgtttgattagctcataaatactaaataatgtttgtaacttatagttacgttttttgg     <  1:3365897‑M1/91‑88 (MQ=255)
                                                                                        aGCcaggatcaaactctccataaaaattatgatgtttgattagctcataaatactaaataatgtttgtaacttatagttacgttttttgga    >  1:1317392‑M1/1‑3 (MQ=255)
                                                                                        aGCcaggatcaaactctccataaaaattatgatgtttgattagctcataaatactaaataatgtttgtaacttatagttacgttttttgga    <  1:1311584‑M1/91‑89 (MQ=255)
                                                                                         gCcaggatcaaactctccataaaaattatgatgtttgattagctcataaatactaaataatgtttgtaacttatagttacgttttttggaa   <  1:2753972‑M1/91‑90 (MQ=255)
                                                                                         gCcaggatcaaactctccataaaaattatgatgtttgattagctcataaatactaaataatgtttgtaacttatagttacgttttttggaa   >  1:50502‑M1/1‑2 (MQ=255)
                                                                                         gCcaggatcaaactctccataaaaattatgatgtttgattagctcataaatactaaataatgtttgtaacttatagttacgttttttggaa   <  1:190571‑M1/91‑90 (MQ=255)
                                                                                         gCcaggatcaaactctccataaaaattatgatgtttgattagctcataaatactaaataatgtttgtaacttatagttacgttttttggaa   <  1:2753995‑M1/91‑90 (MQ=255)
                                                                                         gCcaggatcaaacgctccataaaaattatgatgtttgattagctcataaatactaaataatgtttgtaacttatagttacgttttttggaa   <  1:2753643‑M1/91‑90 (MQ=255)
                                                                                          ccaggatcaaactctccataaaaattatgatgtttgattagctcataaatactaaataatgtttgtaacttatagttacgttttttggaat  <  1:3565725‑M1/91‑91 (MQ=255)
                                                                                          ccaggatcaaactctccataaaaattatgatgtttgattagctcataaatactaaataatgtttgtaacttatagttacgttttttggaat  <  1:3022194‑M1/91‑91 (MQ=255)
                                                                                          ccaggatcaaactctccataaaaattatgatgtttgattagctcataaatactaaataatgtttgtaacttatagttacgttttttggaat  >  1:2286971‑M1/1‑1 (MQ=255)
                                                                                          |                                                                                          
AAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATAATTTTTATGGAGAGTTTGATCCTGGCTCAGGATGAACGCTGGCGGCGTGCCTAATACATGCAAGTCGAGCGAACG  >  USA300TCH1516_ALE/512677‑512857

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