New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? USA300TCH1516_ALE 42499 =235 (1.150)19 (0.120)
+ATAAATTTTTTGATTATTTT
9/120 NT 17.8% coding (1467/1524 nt) USA300TCH1516_00035 hypothetical protein
?USA300TCH1516_ALE = 1912410 0 (0.000)coding (532/609 nt) USA300HOU_RS15290 hypothetical protein

ACGTCAGTTCGAGGAATAAATAAAGTTAAACGAGAGCTAGGTTTTGTATTAATGGCACTTAATATAAGGAAAATAGCAGCTCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTA  >  USA300TCH1516_ALE/1912320‑1912477
                                                                                           |                                                                   
aCGTCAGTTCGAGGAATAAATAAAGTTAAACGAGAGCTAGGTTTTGTATTAATGGCACTTAATATAAGGAAAATAGCAGCTCAACGAGCTg                                                                     <  1:1611225/91‑1 (MQ=32)
 cGTCAGTTCGAGGAATAAATAAAGTTAAACGAGAGCTAGGTTTTGTATTAATGGCACTTAATATAAGGAAAATAGCAGCTCAACGAGCTGt                                                                    <  2:2996323/91‑1 (MQ=31)
                                                  aaTGGCACTTAATATAAGAAAAGTAGTAGCTCAACGAGCTGaaaataatcaaaaaatttataaaaaagacaatttctattgaaataatata                   >  2:329138‑M2/1‑41 (MQ=255)
                                                  aaTGGCACTTAATATAAGAAAAGTAGTAGCTCAACGAGCTGaaaataatcaaaaaatttataaaaaagacaatttctattgaaataatata                   >  2:3704203‑M2/1‑41 (MQ=255)
                                                  aaTAGCACTTAATATAAGAAAAGTAGTAGCTCAACGAGCTGaaaataatcaaaaaatttataaaaaagacaatttctattgaaataatata                   >  2:3831091‑M2/1‑41 (MQ=255)
                                                   aTGGCACTTAATATAAGAAAAGTAGTAGCTCAACGAGCTGaaaataatcaaaaaatttataaaaaagacaatttctattgaaataatatag                  >  1:919242‑M2/1‑40 (MQ=255)
                                                      gCACTTAATATAAGAAAAGTAGTAGCTCAACGAGCTGaaaataatcaaaaaatttataaaaaagacaatttctattgaaataatatagaaa               >  2:1984925‑M2/1‑37 (MQ=255)
                                                      gCACTTAATATAAGAAAAGTAGTAGCTCAACGAGCTGaaaataatcaaaaaatttataaaaaagacaatttctattgaaataatatagaaa               >  2:2921893‑M2/1‑37 (MQ=255)
                                                      gCACTTAATATAAGAAAAGTAGTAGCTCAACGAGCTGaaaataatcaaaaaatttataaaaaagacaatttctattgaaataatatagaaa               >  2:2922004‑M2/1‑37 (MQ=255)
                                                       cACTTAATATAAGAAAAGTAGTAGCTCAACGAGCTGaaaataatcaaaaaatttataaaaaagacaatttctattgaaataatatagaaat              >  2:3021596‑M2/1‑36 (MQ=255)
                                                       aacTTAATATAAGAAAAGTAGTAGCTCAACGAGCTGaaaataatcaaaaaatttataaaaaagacaatttctattgaaataatatagaaat              >  2:1663003‑M2/2‑36 (MQ=255)
                                                        aCTTAATATAAGAAAAGTAGTAGCTCAACGAGCTGaaaataatcaaaaaatttataaaaaagacaatttctattgaaataatatagaaatt             >  2:73521‑M2/1‑35 (MQ=255)
                                                        aCTTAATATAAGAAAAGTAGTAGCTCAACGAGCTGaaaataatcaaaaaatttataaaaaagacaatttctattgaaataatatagaaatt             >  2:73485‑M2/1‑35 (MQ=255)
                                                        aCTTAATATAAGAAAAGTAGTAGCTCAACGAGCTGaaaataatcaaaaaatttataaaaaagacaatttctattgaaataatatagaaatt             >  2:506469‑M2/1‑35 (MQ=255)
                                                         cTTAATATAAGAAAAGTAGTAGCTCAACGAGCTGaaaataatcaaaaaatttataaaaaagacaatttctatagaaataatatagaaattg            >  2:2041147‑M2/1‑34 (MQ=255)
                                                          ttAATATAAGAAAAGTAGTAGCTCAACGAGCTGaaaataatcaaaaaatttataaaaaagacaatttctattgaaataatatagaaattgt           >  2:766059‑M2/1‑33 (MQ=255)
                                                          ttAATATAAGAAAAGTAGTAGCTCAACGAGCTGaaaataatcaaaaaatttataaaaaagacaatttctattgaaataatatagaaattgt           >  2:1545376‑M2/1‑33 (MQ=255)
                                                                taAGAAAAGTAGTAGCTCAACGAGCTGaaaataatcaaaaaatttataaaaaagacaatttctattgaaataatatagaaattgtcttttt     >  2:2145829‑M2/1‑27 (MQ=255)
                                                                 aaGAAAAGTAGTAGCTTAACGAGCTGaaaataatcaaaaaatttataaaaaagacaatttctattgaaataatatagaaattgtctttttt    >  2:2131496‑M2/1‑26 (MQ=255)
                                                                 aaGAAAAGTAGTAGCTCAACGAGCTGaaaataatcaaaaaatttataaaaaagacaatttctattgaaataatatagaaattgtctttttt    >  2:2980946‑M2/1‑26 (MQ=255)
                                                                    gAAAAGTAGTAGCTCAGCGAGCTGaaaataatcaaaaaatttataaaaaagacaatttctattgaaataatatagaaattgtcttttttat  >  2:3779251‑M2/1‑24 (MQ=255)
                                                                                           |                                                                   
ACGTCAGTTCGAGGAATAAATAAAGTTAAACGAGAGCTAGGTTTTGTATTAATGGCACTTAATATAAGGAAAATAGCAGCTCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTA  >  USA300TCH1516_ALE/1912320‑1912477

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