New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? USA300TCH1516_ALE 512620 =9 (0.650)10 (1.070)
+TTATCTTGTCCATTTTATTTTTTAACCAAAATTTGATTAAAAAAC
7/186 NT 55.3% intergenic (+688/‑6075) recR/speA_2 Recombination protein RecR/Arginine decarboxylase
?USA300TCH1516_ALE 2176447 = 15 (1.080)intergenic (+209/+5086) ilvA/USA300TCH1516_02069 L‑threonine dehydratase biosynthetic IlvA/tRNA‑Gly

ACATTTCTCGTTTCGTCAGATTCAAACGTTTTCACTTCGCCAAGCCATCTTTCTTTGTGTTTGCTTTTATTTTGACGTTTTAGACATAAAAAAAGAGACCTTGCGGTCTCAATGCGGCTCATCGCATCCACTTTTTGCCTGGCAACGTTCTACTCTAGCGGAACGTAAGTTCGACTACCATCGACGCTAAGGAGCTTAACTTCTGTGTTCGGCATGGGAACAGGTGTGACCTCCTTGCTATAGTCACCAGACATAT  >  USA300TCH1516_ALE/2176312‑2176567
                                                                                                                                       |                                                                                                                        
aCATTTCTCGTTTCGTCAGATTCAAACGTTTTCACTTCGCCAAGCCATCTTTCTTTGTGTTTGCTTTTATTTTGACGTTTTAGACATAAAAAAAGAGACCTTGCGGTCTCAATGCGGCTCATCGCATCCACTTTTTGCCTGGCAACGtt                                                                                                             <  1:159009/149‑1 (MQ=255)
  aTTTCTCGTTTCGTCAGATTCAAACGTTTTCACTTCGCCAAGCCATCTTTCTTTGTGTTTGCTTTTATTTTGACGTTTTAGACATAAAAAAAGAGACCTTGCGGTCTCAATGCGGCTCATCGCATCCACTTTTTGCCTGGCAACGTTCt                                                                                                           >  1:204184/1‑149 (MQ=40)
             cGTCAGATTCAAACGTTTTCACTTCGCCAAGCCATCTTTCTTTGTGTTTGCTTTTATTTTGACGTTTTAGACATAAAAAAAGAGACCTTGCGGTCTCAATGCGGCTCATCGCATCCACTTTTTGCCTGGCAACGTTCTACTCTAGCGGa                                                                                                <  1:188742/149‑1 (MQ=40)
               agttttcaatgttcatttttcttaccgacaagaattaattatacattttatcaacatttaagtcaataactttttttatcttgtccattttattttttaaccaaaatttgattaaaaaacTGCCTGGCAACGTTCTACTCTAGCGGAAc                                                                                              >  1:281285‑M2/121‑149 (MQ=255)
               agttttcaatgttcatttttcttaccgacaagaattaattatacattttatcaacatttaagtcaataactttttttatcttgtccattttattttttaaccaaaatttgattaaaaaacTGCCTGGCAACGTTCTACTCTAGCGGAAc                                                                                              >  1:100557‑M2/121‑149 (MQ=255)
                        aaCGTTTTCACTTCGCCAAGCCATCTTTCTTTGTGTTTGCTTTTATTTTGACGTTTTAGACATAAAAAAAGAGACCTTGCGGTCTCAATGCGGCTCATCGCATCCACTTTTTGCCTGGCAACGTTCTACTCTAGCGGAACGTAAGTTCg                                                                                     <  1:215924/149‑1 (MQ=39)
                              tttttcttaccgacaagaattaattatacattttatcaacatttaagtcaataactttttttatcttgtccattttattttttaaccaaaatttgattaaaaaacTGCCTGGCAACGTTCTACTCTAGCGGAACGTAAGTTCGACTAcc                                                                               >  1:147553‑M2/106‑149 (MQ=255)
                              tttttcttaccgacaagaattaattatacattttatcaacatttaagtcaataactttttttatcttgtccattttattttttaaccaaaatttgattaaaaaacTGCCTGGCAACGTTCTACTCTAGCGGAACGTAAGTTCGACTAcc                                                                               >  1:215432‑M2/106‑149 (MQ=255)
                                     taccgacaagaattaattatacattttatcaacatttaagtcaataactttttttatcttgtccattttattttttaaccaaaatttgattaaaaaacTGCCTGGCAACGTTCTACTCTAGCGGAACGTAAGTTCGACTACCATCGACg                                                                        <  1:29206‑M2/51‑1 (MQ=255)
                                      gCCAAGCCATCTTTCTTTGTGTTTGCTTTTATTTTGACGTTTTAGACATAAAAAAAGAGACCTTGCGGTCTCAATGCGGCTCATCGCATCCACTTTTTGCCTGGCAACGTTCTACTCTAGCGGAACGTAAGTTCGACTACCATCGACGc                                                                       <  1:293508/149‑1 (MQ=40)
                                        cgacaagaattaattatacattttatcaacatttaagtcaataactttttttatcttgtccattttattttttaaccaaaatttgattaaaaaacTGCCTGGCAACGTTCTACTCTAGCGGAACGTAAGTTCGACTACCATCGACGCTa                                                                     >  1:79065‑M2/96‑149 (MQ=255)
                                              aTCTTTCTTTGTGTTTGCTTTTATTTTGACGTTTTAGACATAAAAAAAGAGACCTTGCGGTCTCAATGCGGCTCATCGCATCCACTTTTTGCCTGGCAACGTTCTACTCTAGCGGAACGTAAGTTCGACTACCATCGACGCTAAGGAGc                                                               <  1:226126/149‑1 (MQ=39)
                                               tctttcttTGTGTTTGCTTTTATTTTGACGTTTTAGACATAAAAAAAGAGACCTTGCGGTCTCAATGCGGCTCATCGCATCCACTTTTTGCCTGGCAACGTTCTACTCTAGCGGAACGTAAGTTCGACTACCATCGACGCTAAGGAGCt                                                              <  1:18489/149‑1 (MQ=38)
                                                      ttGTGTTTGCTTTTATTTTGACGTTTTAGACATAAAAAAAGAGACCTTGCGGTCTCAATGCGGCTCATCGCATCCACTTTTTGCCTGGCAACGTTCTACTCTAGCGGAACGTAAGTTCGACTACCATCGACGCTAAGGAGCTTAACTTc                                                       <  1:197324/149‑1 (MQ=37)
                                                             tttatcaacatttaagtcaataactttttttatcttgtccattttattttttaaccaaaatttgattaaaaaacTGCCTGGCAACGTTCTACTCTAGCGGAACGTAAGTTCGACTAGCATCGACGCTAAGGAGCTTAACTTCTGTGTTc                                                >  1:171568‑M2/75‑149 (MQ=255)
                                                             tttatcaacatttaagtcaataactttttttatcttgtccattttattttttaaccaaaatttgattaaaaaacTGCCTGGCAACGTTCTACTCTAGCGGAACGTAAGTTCGACTACCATCGACGCTAAGGAGCTTAACTTCTGTGTTc                                                >  1:278941‑M2/75‑149 (MQ=255)
                                                                ttttATTTTGACGTTTTAGACATAAAAAAAGAGACCTTGCGGTCTCAATGCGGCTCATCGCATCCACTTTTTGCCTGGCAACGTTCTACTCTAGCGGAACGTAAGTTCGACTACCATCGACGCTAAGGAGCTTAACTTCTGTGTTCGGc                                             >  1:23921/1‑149 (MQ=32)
                                                                ttttATTTTGACGTTTTAGACATAAAAAAAGAGACCTTGCGGTCTCAATGCGGCTCATCGCATCCACTTTTTGCCTGGCAACGTTCTACTCTAGCGGAACGTAAGTTCGACTACCATCGACGCTAAGGAGCTTAACTTCTGTGTTCGGc                                             >  1:126520/1‑149 (MQ=32)
                                                                atcaacatttaagtcaataactttttttatcttgtccattttattttttaaccaaaatttgattaaaaaacTGCCTGGCAACGTTCTACTCTAGAAAGTCTTAAGTTCGACTACCATCGACGCTAAGGAGCTTAACTTCTGTGTTCGGc                                             <  1:296459‑M2/78‑1 (MQ=255)
                                                                        ttaagtcaataactttttttatcttgtccattttattttttaaccaaaatttgattaaaaaacTGCCTGGCAACGTTCTACTCTAGCGGAACGTAAGTTCGACTACCATCGACGCTAAGGAGCTTAACTTCTGTGTTCGGCATGGGAAc                                     >  1:277274‑M2/64‑149 (MQ=255)
                                                                          aCGTTTTAGACATAAAAAAAGAGACCTTGCGGTCTCAATGCGGCTCATCGCATCCACTTTTTGCCTGGCAACGTTCTACTCTAGCGGAACGTAAGTTCGACTACCATCGACGCTAAGGAGCTTAACTTCTGTGTTCGGCATGGGAACAg                                   <  1:229117/149‑1 (MQ=35)
                                                                                          ttatcttgtccattctattttttaaccaaaatttgattaaaaaacTGCCTGGCAACGTTCTACTCTAGCGGAACGTAAGTTCGACTACCATCGACGCTAAGGAGCTTAACTTCTGTGTTCGGCATGGGAACAGGTGTGACCTCCTTGCt                   <  1:259295‑M2/104‑1 (MQ=255)
                                                                                                        ggTCTCAATGCGGCTCATCGCATCCACTTTTTGCCTGGCAACGTTCTACTCTAGCGGAACGTAAGTTCGACTACCATCGACGCTAAGGAGCTTAACTTCTGTGTTCGGCATGGGAACAGGTGTGACCTCCTTGCTATAGTCACCAGACa     <  1:255633/149‑1 (MQ=32)
                                                                                                        ggTCTCAATGCGGCTCATCGCATCCACTTTTTGCCTGGCAACGTTCTACTCTAGCGGAACGTAAGTTCGACTACCATCGACGCTAAGGAGCTTAACTTCTGTGTTCGGCATGGGAACAGGTGTGACCTCCTTGCTATAGTCACCAGACa     <  1:255663/149‑1 (MQ=32)
                                                                                                         gTCTCAATGCGGCTCATCGCATCCACTTTTTGCCTGGCAACGTTCTACTCTAGCGGAACGTAAGTTCGACTACCATCGACGCTAAGGAGCTTAACTTCTGTGTACGGCATGGGAACAGGTGTGACCTCCTTGCTATAGTCACCAGACat    <  1:218993/149‑1 (MQ=17)
                                                                                                           ctcAATGCGGCTCATCGCATCCACTTTTTGCCTGGCAACGTTCTACTCTAGCGGAACGTAAGTTCGACTACCATCGACGCTAAGGAGCTTAACTTCTGTGTTCGGCATGGGAACAGGTGTGACCTCCTTGCTATAGTCACCAGACatat  <  1:164171/149‑1 (MQ=32)
                                                                                                                                       |                                                                                                                        
ACATTTCTCGTTTCGTCAGATTCAAACGTTTTCACTTCGCCAAGCCATCTTTCTTTGTGTTTGCTTTTATTTTGACGTTTTAGACATAAAAAAAGAGACCTTGCGGTCTCAATGCGGCTCATCGCATCCACTTTTTGCCTGGCAACGTTCTACTCTAGCGGAACGTAAGTTCGACTACCATCGACGCTAAGGAGCTTAACTTCTGTGTTCGGCATGGGAACAGGTGTGACCTCCTTGCTATAGTCACCAGACATAT  >  USA300TCH1516_ALE/2176312‑2176567

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