New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? USA300TCH1516_ALE 36148 =NA (NA)4 (0.040) 4/268 12.7 5.0% coding (404/675 nt) USA300HOU_RS00140 hypothetical protein
?USA300TCH1516_ALE = 76632 80 (0.750)intergenic (+172/‑444) USA300HOU_RS00140/USA300HOU_RS00350 hypothetical protein/Monoacylglycerol lipase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATC  >  USA300TCH1516_ALE/36010‑36272
                                                                                                                                          |                                                                                                                            
tcgtACCTTGTCTTTCTTACTTTAATATGACGGTGATCGTGCTCAATGAGGTTATTGAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGcc                                                                                                                            <  2:700159/138‑1 (MQ=255)
tGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTCAGCCATTGCTACCTTCGTTGAAGGTGcc                                                                                                                            <  2:913898/141‑1 (MQ=25)
   tcccTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCACTGCGGTTGTTAGGATGTGGCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGa                                                                                                                         <  2:546461/139‑1 (MQ=11)
   tACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGa                                                                                                                         <  2:5322/141‑1 (MQ=35)
      cTTGTCTTTTTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCt                                                                                                                      >  1:883533/1‑141 (MQ=25)
      cTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCt                                                                                                                      >  1:757325/1‑141 (MQ=34)
      cTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCt                                                                                                                      <  1:1042344/141‑1 (MQ=35)
            tttcttACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAAtt                                                                                                                >  2:789982/1‑141 (MQ=35)
            tttcttACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAAtt                                                                                                                <  1:1104869/141‑1 (MQ=35)
             ttcttACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTa                                                                                                               <  1:151300/141‑1 (MQ=35)
             ttcttACTTTAATATGACGGTGATCTTGCTCAATGAAGTTATACAAATATTACGATGGACAATGACAGTCAGGGTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTCGAAGGTGCCGGTTCTGttgttt                                                                                                               >  1:113969/1‑136 (MQ=11)
                ttACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCt                                                                                                            <  2:334109/141‑1 (MQ=35)
                   cTTTAATATGGCGGGTATCTTGCGCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCtttt                                                                                                         <  2:1058077/141‑1 (MQ=255)
                   cTTTAATATGACGGAGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCtttt                                                                                                         >  1:656053/1‑141 (MQ=21)
                    tttACTAGGACGGTGATCTTGCTCAATGGGGTTATTCAGATATTTCGATGGACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTg                                                                                                        <  2:788820/141‑1 (MQ=255)
                    tttAATTTGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTg                                                                                                        <  2:599647/141‑1 (MQ=21)
                     ttAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCCTCGTTGAAGGTGCCTGATCTGTAATTACCTGTTGa                                                                                                       <  2:141199/141‑1 (MQ=21)
                           tGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTa                                                                                                 <  2:679001/141‑1 (MQ=35)
                             aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAcc                                                                                               >  2:1187243/1‑141 (MQ=35)
                             aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAcc                                                                                               >  1:481779/1‑141 (MQ=34)
                             aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAcc                                                                                               >  2:1060815/1‑141 (MQ=35)
                             aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAcc                                                                                               >  1:1093058/1‑141 (MQ=35)
                                   ttCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTg                                                                                         <  2:52710/140‑1 (MQ=17)
                                   aTCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTg                                                                                         >  2:1185006/1‑141 (MQ=32)
                                    tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGt                                                                                        <  1:87997/141‑1 (MQ=32)
                                    tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGt                                                                                        >  1:586224/1‑141 (MQ=32)
                                      ttGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGttt                                                                                      >  1:146473/1‑141 (MQ=32)
                                      ttGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGttt                                                                                      >  1:502834/1‑141 (MQ=32)
                                      ttGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCGGATCTGTAATTACCTTTTGCGGTTTACAAAATTGctt                                                                                      >  1:691313/1‑141 (MQ=255)
                                        gCTCAATGGGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTaa                                                                                    <  2:1135009/141‑1 (MQ=17)
                                           cAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTTAAGGTGCCTGTTCTGTAATTACCTTTTTATGTTTACCAAATTGTTTAATga                                                                                 >  1:1141944/1‑141 (MQ=255)
                                                 ggTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGtt                                                                           <  1:704236/141‑1 (MQ=32)
                                                   ttATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTCTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTg                                                                         <  2:656053/141‑1 (MQ=17)
                                                   ttATTCAGATATTTCGATGTACAATGACAGGCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTg                                                                         <  2:601095/141‑1 (MQ=17)
                                                   ttATCCAGATAGTTCGATGTACAATGCCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATTAGACGTTTg                                                                         <  1:441737/141‑1 (MQ=255)
                                                       tCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTGTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATaa                                                                     <  2:904970/141‑1 (MQ=17)
                                                        cAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTGGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTCCCAAATTGTTTAATGAAACGTTGGATaaa                                                                    >  2:1072243/1‑141 (MQ=255)
                                                        cAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTCATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATaaa                                                                    >  1:1025290/1‑141 (MQ=17)
                                                           atatTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGc                                                                 >  2:304642/1‑141 (MQ=32)
                                                           atatTTCCTTGTACAATGACTGTGCTTTTTAAGTCTAATAGCTTTCTTTACTGTGGCCGTTGCTACCTTCGTTGAAGGTGTCTGATCTGTAATCACCTTTTGAGGTTTGCCAAATTGTTTAGTGAGGCGTTTGATTAACGc                                                                 <  2:1092287/141‑1 (MQ=2)
                                                              tttCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCata                                                              >  1:1118832/1‑141 (MQ=32)
                                                                tCGATTTACAATGACAGTCAGGGTTTAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATg                                                            <  2:934704/141‑1 (MQ=12)
                                                                tCGATGTACAATGACAGTGAGGTTGAAGTTTAAAAGCTTTAATCACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATg                                                            <  1:1204720/141‑1 (MQ=12)
                                                                tCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATg                                                            >  2:1000384/1‑141 (MQ=32)
                                                                tCGAGGTACAATGCCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATg                                                            <  1:3169/141‑1 (MQ=17)
                                                                   aTGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTg                                                         >  1:182639/1‑141 (MQ=32)
                                                                    tGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGa                                                        >  1:210675/1‑141 (MQ=32)
                                                                       aCAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATg                                                     >  1:673110/1‑141 (MQ=32)
                                                                         aaTGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGAt                                                   <  1:789982/141‑1 (MQ=32)
                                                                           tGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTa                                                 >  2:983029/1‑141 (MQ=31)
                                                                            gACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACACATATGCTGAATGATTAt                                                >  1:1114260/1‑141 (MQ=11)
                                                                            gACAGTAAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGACAATGATGCCTTCGTTGAAGGCGCCTGATCTGTACTTACCTTTTGTGGTTTAACAAATTGTTTACTGAGACGTTTGATAAATGCATATGCTGAATGATTAt                                                >  2:621575/1‑141 (MQ=11)
                                                                                                                tAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTc            >  1:1147431/1‑141 (MQ=32)
                                                                                                                  gCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTccc          >  1:406776/1‑141 (MQ=11)
                                                                                                                         cttccTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGCGGTTTACCCAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCAt   <  1:543761/138‑1 (MQ=0)
                                                                                                                          taacTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACTGTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCAtt  <  1:380604/138‑2 (MQ=0)
                                                                                                                                          |                                                                                                                            
TGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATC  >  USA300TCH1516_ALE/36010‑36272

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