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.

CAACACGATGTAAAAGGTGATGTTTATCAATCTTTATCACAATTACCTGAGTCAGTAAGAGAAATTATTGACGTAGATCATTCTAAAAACTTAGTAATTCAACATAATAATACGATTGCGTACACACAAGTTGATGATAATGCATCGAAAGATGGCGTTATCGTTGAAGGTTTAGCAGACGCTCTTATGAACCATAGTGATTTAGTACAAAAGTACTTTATGAAAGATGCAGTAACAGTAGATGAACATCGTATCACAGCGCTACACACGGCATTAG  >  USA300TCH1516_ALE/912753‑913029
                                                                                                                                           |                                                                                                                                          
caacaCGATGTAAAAGGTGATGTTTATCAATCTTTATCACAATTACCTGAGTCAGTAAGAGAAATTATTGACGTAGATCATTCTAAACACTTAGTAATTCAACATAATAATACGATTGCGTACACACAAGTTGATGATAA                                                                                                                                            >  1:83527‑M1/1‑140 (MQ=255)
  acacGATGTAAAAGGTGATGTTTATCAATCTTTATCACAATTACCTGAGTCAGTAAGAGAAATTATTGACGTAGATCATTCTAAAAACTTAGTAATTCAACATAATAATACGATTGCGTACACACAAGTTGATGATAATg                                                                                                                                          >  2:145747/1‑140 (MQ=255)
                    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)
                                    tCACAATTACCTGAGTCAGTAAGAGAAATTATTGACGTAGATAATTCTATATACTTAGTAATTCAACATAATAATACGATTGCGTACACACAAGTTGATGATAAcgccatctttcgatgcattatcatcaacttgtgtgt                                                                                                        <  1:86552‑M1/140‑37 (MQ=255)
                                               tGAGTCAGTAAGAGAAATTATTGACGTAGATCATTCTAAAAACTTAGTAATTCAACATAATAATACGATTGCGTACACACAAGTTGATGATAACGCCATCTTTCGATGCGttatcatcaacttgtgtgtacgcaatcgta                                                                                            <  2:86653‑M2/140‑31 (MQ=255)
                                               tGAGTCAGTAAGAGAAATTATTGACGTAGATCATTCTAAAAACTTAGTAATTCAACATAATAATACGATTGCGTACACACAAGCAGATGATAGctcatcggaagatggcgttatcatcaacttgtgtgtacgcaatcgta                                                                                             >  1:86653‑M1/1‑93 (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)
                                                                                                              taCGATTGCGTACACACAAGTTGATGATAAcgcatcgaaagatggcgttatcatcaacttgtgtgtacgcaatcgtattattatgttgaattactaagtttttagaatgatctacgtcaataatttctcttactgactca                              <  2:86653‑M1/140‑111 (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)
                                                                                                                                         tAAcgccatctttcgatgcattatcatcaacttgtgtgtacgcaatcgtattattatgttgaattactaagtttttagaatgatctacgtcaataatttctcttactgactcaggtaattgtgctcaagctagagaaaca   >  1:60506‑M1/1‑3 (MQ=255)
                                                                                                                                         tAATGCATCGAAAGATGGCGTTATCGTTGTTGTTTTTGCAGACACTCTTATGAACCATAGTGATTTAGTACAAAAGTACTTTATGAAAGATGCAGTAACAGTAGATGAACATCGTATCACAGCGCTACACACGGCATTAg  <  1:8620/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                          
CAACACGATGTAAAAGGTGATGTTTATCAATCTTTATCACAATTACCTGAGTCAGTAAGAGAAATTATTGACGTAGATCATTCTAAAAACTTAGTAATTCAACATAATAATACGATTGCGTACACACAAGTTGATGATAATGCATCGAAAGATGGCGTTATCGTTGAAGGTTTAGCAGACGCTCTTATGAACCATAGTGATTTAGTACAAAAGTACTTTATGAAAGATGCAGTAACAGTAGATGAACATCGTATCACAGCGCTACACACGGCATTAG  >  USA300TCH1516_ALE/912753‑913029

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.