New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? USA300TCH1516_ALE = 91289221 (0.920)3 (0.140) 3/268 2.7 12.4% coding (317/1308 nt) sufB_1 FeS cluster assembly protein SufB
?USA300TCH1516_ALE = 912909 22 (1.000)coding (334/1308 nt) sufB_1 FeS cluster assembly protein SufB
Rejected: Frequency below/above cutoff threshold.

TGTTTATCAATCTTTATCACAATTACCTGAGTCAGTAAGAGAAATTATTGACGTAGATCATTCTAAAAACTTAGTAATTCAACATAATAATACGATTGCGTACACACAAGTTGATGATAATGCATCGAAAGATGGCGTTATCGTTGAAGGTTTAGCAGACGCTCTTATGAACCATAGTGATTTAGTACAAAAGTACTTTATGAAAGATGCAGTAACAGTAGATGAACATCGTATCACAGCGCTACACACGGCATTAGTTAATGGTGGCG  >  USA300TCH1516_ALE/912773‑913041
                                                                                                                                         |                                                                                                                                    
tGTTTCTCTAGCTTGAGCACAATTACCTGAGTCAGTAAGAGAAATTATTGACGTAGATCATTCTAAAAACTTAGTAATTCAACATAATAATACGATTGCGTACACACAAGTTGATGATAATGCATCGAAAGATGGCGtta                                                                                                                                   >  1:60506‑M2/1‑137 (MQ=255)
        aaTCTTTATCACAATTACCTGAGTCAGTAAGAGAAATTATTGACGTAGATCATTCTAAAAACTTAGTAATTCAACATAATAATACGATTGCGTAAACACAAGTTGATGATAATGCATCGAAAGATGGCGTTATCGTTGaa                                                                                                                           >  1:644/1‑140 (MQ=255)
          tCTGTAACACAATTACCTGAGGCAGTAAGAGAAATTATTGACGTAGATCATTCTAAAAACTTAGTAATTCAACATAATAATACGATTGCGTACACACAAGTTGATGATAATGCATCGAAAGATGGCGTTATCGTTGAAgg                                                                                                                         <  1:109018/140‑1 (MQ=255)
             ttATCACAATTACCTGAGTCAGTAAGAGAAATTATTGACGTAGATCATTCTAAAAACTTAGTAATTCAACATAATAATACGATTGCGTACACACAAGTTGATGATAATGCATCGAAAGATGGCGTTATCGTTGAAGGttt                                                                                                                      <  1:161103/140‑1 (MQ=255)
               attaaaaTTACCTGAGTCAGTAAGAGAAATTATTGATGTAGATTATTCTAAAAACTTAGTAATTCAACATAATAATACGATTGCGTACACACAAGTTGATGATAATGCATCGAAAGATGGCGTTATCGTTGAAGGTTTAg                                                                                                                    <  2:224763/135‑1 (MQ=255)
                           tGAGTCAGTAAGAGAAATTATTGACGTAGATCATTCTAAAAACTTAGTAATTCAACATAATAATACGATTGCGTACACACAAGTTGATGATAACGCCATCTTTCGATGCGttatcatcaacttgtgtgtacgcaatcgta                                                                                                        <  2:86653‑M2/140‑31 (MQ=255)
                             aGTCAGTAAGAGAAATTATTGACGTAGATCATTCTAAAAACTTAGTAATTCAACATAATAATACGATTGCGTACACACAAGTTGATGATAATGCATCGAAAGATGGCGTTATCGTTGAAGGTTTAGCAGACGCTCTTATg                                                                                                      <  1:82324/140‑1 (MQ=255)
                                           attaatGACGTAGATCATTCTAAAAACTTAGTAATTCAACATAATAATACGATTGCGTACACACAAGTTGATGATAATGCGTCGAAAGATGGCGTTATCGTTGAAGGTTTAGCAGACGCTCTTATGAACCATAGTGAttt                                                                                        <  2:195110/140‑1 (MQ=255)
                                               ttGACGTAGATCATTCTAAAAACTTAGTAATTCAACATAATAATACGATTTCGTACACACAAGTTGATGATAATGCATCGAAAGATAGCGTTGTCGTTGAAGGTTTAGCAGACGCTCTTATGAACCATAGTGATTTAGTa                                                                                    >  2:23319/1‑140 (MQ=255)
                                                                      ttAGTAATTCAACATAATAATACGATTGCGTACACACAAGTTGATGATAATACATAGAAAGATGGCGTTATCGTTGAAGGTTTAGCCCACGCACGTATGAACAATAGTCATTTAGTACAAAAGTACTTTATGAAAGAttc                                                             >  2:159885/1‑138 (MQ=255)
                                                                        aGTAATTCAACATAATAATACGATTGCGTACACACAAGTTGATGATAATGCATCGAAAGCTGGCGGTATCGTTGAAGGTTTAGCAGACACTCTTATGAACCATAGTGATTTAGTACAAAAGTACTTTATGAAAGATGCAg                                                           >  2:16239/1‑140 (MQ=255)
                                                                        aGTAATACAACATAATAATACGATTGCGTACACACAAGTTGATGATAATGCATCGAAAGATGGCGTTATCGTTGAAGGTTTAGCAGACGCTCTTATGAACCATAGTGATTTAGTACAAAAGTACTTTATGAAAGATGCAg                                                           <  1:85575/140‑1 (MQ=255)
                                                                                  cATAATAATACGATTGCGTACACACAAGTTGATGATAATGCATCGAAAGATGGCGTTATCGTTGAAGGTTTAGCAGACGCTCTTATGACCCATAGTGATTTAGTACACAAGTACTTTCTGAAAGATGCAGTCACAGTAGa                                                 >  1:234367/1‑140 (MQ=255)
                                                                                  cATAATAATACGATTGCGTACACACAAGTTGATGATAATGCATCGAAAGATGGCGTTATCGTTGAAGGTTTAGCAGACGCTCTTATGAACCATAGAGATTTAGTACAAAAGTACTTTATGAACGATGCAGTAACAGTAGa                                                 >  2:207637/1‑140 (MQ=255)
                                                                                      ataataCGATTGCGTACACACAAGTTGATGATAATGCATCGAAAGATGGCGTTATCGTTGAAGGTTTAGCAGACGCTCTTATGAACCATAGTGATTTAGTACAAAAGTACTTTATGAAAGATGCAGTAACAGTAGATGaa                                             >  2:254233/1‑140 (MQ=255)
                                                                                      ataataCGATTGCGTACACACAAGTTGATGATAATGCATCGAAAGATGGCGTTATCGTTGAAGGTTTAGCAGACGCTCTTATGAACCATAGTGATTTAGTACAAAAGAAATAAATGGAAGATGCAGTAACAGTAGATGaa                                             >  2:1424/1‑140 (MQ=255)
                                                                                          taCGATTGCGTACACACAAGTTGATGATAACGCCATCTTCCGATGAGctatcatctgcttgtgtgtacgcaatcgtattattatgttgaattactaagtttttagaatgatctacgtcaataatttctcttactgactca                                         >  1:86653‑M2/1‑47 (MQ=255)
                                                                                                  cGTACACACAAGTTGATGATAATGCATCGAAAGATGGCGTTATCGTTGAAGGTTTAGCAGACGCTCTTATGAACCATAGTGATTTAGTACAAAAGTACTTTATGAAAGATGCAGTAACAGTAGATGAACATCGTATcaca                                 >  2:32855/1‑140 (MQ=255)
                                                                                                     acacacAAGTTGATGATAATGCATCGAAAGATGGCGttatcatcaacttgtgtgtacgcaatcgtattattatgttgaattactaagtatatagaattatctacgtcaataatttctcttactgactcaggtaattgtga                              <  1:86552‑M2/140‑105 (MQ=255)
                                                                                                         gcaAGTTGATGAGAATGCATCGAAAGATGGCGTTATCGTTGAAGGTTTAGCAGACGCTCTTATGAACCATAGTGATTTAGTACAAAAGTACTTTATGAAATATGCAGTAACAGTAGATGAACATCGTATCACAGCGCTac                          <  2:140516/139‑1 (MQ=255)
                                                                                                             gTTGATGATAATTCATTGTTAGATGGTGTTATCGTTGAAGGTTTAGCAGACGCTCGTATGAACCATAGTGATTTAGTACAAAAGTACTTTATGAAAGATGCAGTAACAGTAGATGAACATCGTATCACAGCGCTacacac                      <  1:206052/140‑1 (MQ=255)
                                                                                                                gatgatAATGCATCGAAAGATGGCGTTATCGTTGAAGGTTTAGCAGACGCTCTTATGAACCATAGTGATTTAGTACAAAAGTACTTTATGAAAGATGCAGTAACAGTAGATGAACATCGTATCACAGCGCTACACACGGc                   >  1:124606/1‑140 (MQ=255)
                                                                                                                    atAATGCATCGAAAGATGGCGTTATCGTTGAAGGTTTAGCAGACGCTCTTATGAACCATAGTGATTTAGTACAAAAGTACTTTATGAAAGATGCAGTAACAGTAGATGAACATCGTATCACAGCGCTACACACGGCATTa               >  1:27575/1‑140 (MQ=255)
                                                                                                                     tAATGCATCGAAAGATGGCGTTATCGTTGTTGTTTTTGCAGACACTCTTATGAACCATAGTGATTTAGTACAAAAGTACTTTATGAAAGATGCAGTAACAGTAGATGAACATCGTATCACAGCGCTACACACGGCATTAg              <  1:8620/140‑1 (MQ=255)
                                                                                                                                aaaGATGGCGTTATCGTTGAAGGTTTAGCAGACGCTCTTATGAACCATAGTGATTTAGTACAAAAGTACTTTATGAAAGATGCAGTAACAGTAGATGAACATCGTATCACAGCGCTACACACGGCATTAGTTAAtggtgg    <  1:145747/140‑1 (MQ=255)
                                                                                                                                  aGATGTCTTTATCGTTGAAGGTTTAGCAGACGCTTTTATGAACCATAGTGATTTAGTACAAAAGTACTTTATGAAAGATGCAGTAACAGTAGATGAACATCGTATCACATCGCTACACACGGCATTAGTTAATGGTGGCg  <  2:70472/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                    
TGTTTATCAATCTTTATCACAATTACCTGAGTCAGTAAGAGAAATTATTGACGTAGATCATTCTAAAAACTTAGTAATTCAACATAATAATACGATTGCGTACACACAAGTTGATGATAATGCATCGAAAGATGGCGTTATCGTTGAAGGTTTAGCAGACGCTCTTATGAACCATAGTGATTTAGTACAAAAGTACTTTATGAAAGATGCAGTAACAGTAGATGAACATCGTATCACAGCGCTACACACGGCATTAGTTAATGGTGGCG  >  USA300TCH1516_ALE/912773‑913041

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