New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? USA300TCH1516_ALE 895785 =26 (1.150)5 (0.250) 5/244 1.9 18.2% coding (806/822 nt) metQ_2 Methionine‑binding lipoprotein MetQ
?USA300TCH1516_ALE = 909740 22 (1.110)intergenic (‑663/‑371) USA300HOU_RS04500/csbD_1 hypothetical protein/Stress response protein CsbD
Rejected: Frequency below/above cutoff threshold.

AAAAGGGCATAATTTTGAAACAAAGGGCAAGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTATTACGGATAAAGTAAAATCCTATATGCATGCGTGCT  >  USA300TCH1516_ALE/909620‑909878
                                                                                                                        |                                                                                                                                          
aaaaGGGTATAATTTTGATATAAAGGGCAAGGGTATGAATAATAAATTATCCTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTaaataaat                                                                                                                         <  1:101458/140‑1 (MQ=255)
  aaGGGCATAATTTTGAAACAAAGGGCAAGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAAttattcctgctaaataatacc                                                                                                                       <  1:13679‑M2/140‑22 (MQ=255)
     ggCATAATTTTGAAACAAAGGGCAAGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCatat                                                                                                                    <  1:240880/140‑1 (MQ=255)
           aTTTTGAAACAAAGGTCAAGGGTATGAATAATAAATTATCTTTAAAATTTTATAAATGCATATAAAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATgg                                                                                                              <  1:2025/140‑1 (MQ=255)
             tttGAAACAAAGGGCAAGGGTATGAATAATAAATTATCTTGAAACTTTTATAAATTCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCt                                                                                                            <  2:226437/140‑1 (MQ=255)
                  aaCAAAGGGCAAGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTa                                                                                                       >  2:237968/1‑140 (MQ=255)
                           cAAGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGt                                                                                              <  1:43991/140‑1 (MQ=255)
                                       taataaATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCtttt                                                                                  <  1:262161/140‑1 (MQ=255)
                                            aaTTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTtttaatt                                                                             >  1:64041/1‑140 (MQ=255)
                                                               taAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTACTGGGCTTTTTAATTTAAACAGGTATTTATATGa                                                          >  2:150318/1‑140 (MQ=255)
                                                                 aaTGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGata                                                        <  2:212586/140‑1 (MQ=255)
                                                                    gCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATAttt                                                     <  1:237968/140‑1 (MQ=255)
                                                                              aaatGTTTTTATTGCCCTTGACCTTCCTTGATAAGCTGGCTAAttattactgctaaataataccaattaaataaaaataccctttaataacaacgattaaggggttatttttgtggtattgaaaataaaaaagggcacaa                                           >  2:34319‑M2/2‑43 (MQ=255)
                                                                               aaTGTTTTAAAGGACCTTGAACTTCCTTGATAATCTGGCTAAttattcctgctaaataataccaattaaataaaaataccctttaataacaacgattaaggggttatttttgtggtattgaaaataaaaaagggcacaaa                                          >  2:156315‑M2/1‑42 (MQ=255)
                                                                                aTGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTTTTTTGAAGAGAAAATAATATGATATTTAGGAATGAGATg                                         >  1:36355/1‑140 (MQ=255)
                                                                                aTGTGTTGATTGCCCTGGAACTTCCTTGATAATCTGGCTAAttattcctgctaaataagaccaattaaataaaaataccctttaataacaacgattaaggggttatttttgtggtattgaaaataaaaaagggcacaaaa                                         >  2:111001‑M2/1‑41 (MQ=255)
                                                                                 tGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGc                                        >  1:246573/1‑140 (MQ=255)
                                                                                  gTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAAttattcctgctaaataataccaattaaataaaaataccctttaataacaacgattaaggggttatttttgtggtattgaaaataaaaaagggcacaaaagg                                       <  2:73569‑M2/140‑102 (MQ=255)
                                                                                       gATTGCACTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTCAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAg                                                                      >  2:119543/1‑104 (MQ=255)
                                                                                       gATTGCACTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAg                                                                      <  1:119543/104‑1 (MQ=255)
                                                                                           gCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTATTACGGAt                              >  1:246015/1‑140 (MQ=255)
                                                                                              cgtGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTATTACGGATaaa                           <  1:6167/138‑1 (MQ=255)
                                                                                                 gAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGCTATTTAGGAATGAGATGCTATTACGGATAAAGTa                        <  1:68367/140‑1 (MQ=255)
                                                                                                   aCTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTATTACGGATAAAGTaaa                      >  1:243916/1‑140 (MQ=255)
                                                                                                          tGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTATTACGGATAAAGTAACATCCtat               >  2:9545/1‑140 (MQ=255)
                                                                                                          tGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGAAATTTAGGAATGAGATGCTATTACGGATAAAGTAAAATCCtat               >  1:95494/1‑140 (MQ=255)
                                                                                                                       aaTTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTATTACGGATAAAGTAAAATCCTATATGCATGCGTGCt  <  2:64041/140‑1 (MQ=255)
                                                                                                                        |                                                                                                                                          
AAAAGGGCATAATTTTGAAACAAAGGGCAAGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTATTACGGATAAAGTAAAATCCTATATGCATGCGTGCT  >  USA300TCH1516_ALE/909620‑909878

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