New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? BX571857 980997 =NA (NA)5 (0.320) 5/118 NT NA noncoding (1/27 nt) Repeat flanking prophage phiSa4(476) Repeat flanking prophage phiSa4(476)
?BX571857 = 1025750 NA (NA)intergenic (+519/‑57) SAS0954/SAS0955 amidase/putative protease

ACACGAGTGGAAAATATTTAAAGTTAAAGTGATAAAAACTTTGATATTACAGCGTTTTTAAACGTTCGTGTACGTTTGAGAATGGCCGTTTTACCATTAGGTGTTAATGTGTCATAATGATAAGTGATTTTATTGAATGGAGTGGACATTAGTGGATATTGGTAAAAAACATGTAATTCCTAAAAGTCAGT  >  BX571857/1025656‑1025846
                                                                                              |                                                                                                
acacGAGTGGAAAATATTTAAAGTTAAAGTGATAAAAACTTTGATATTACAGCGTTTTTAAACGTTCGTGTACGTTTGAGAATGGCCGTTTTACCat                                                                                                >  2:382346‑M2/1‑95 (MQ=255)
acacGAGTGGAAAATATTTAAAGTTAAAGTGATAAAAACTTTGATATTACAGCGTTTTTAAACGTTCGTGTACGTTTGAGAATGGCCGTTTTACCat                                                                                                >  1:306783‑M2/1‑95 (MQ=255)
 cacGAGTGGAAAATATTTAAAGTTAAAGTGATAAAAACTTTGATATTACAGCGTTTTTAAACGTTCGTGTACGTTTGAGAATGGCCGTTTTACCatt                                                                                               >  2:147504‑M2/1‑94 (MQ=255)
  acGAGTGGAAAATATTTAAAGTTAAAGTGATAAAAACTTTGATATTACAGCGTTTTTAAACGTTCGTGTACGTTTGAGAATGGCCGTTTTACCatta                                                                                              >  1:360743‑M2/1‑93 (MQ=255)
            aaTATTTAAAGTTAAAGTGATAAAAACTTTGATATTACAGCGTTTTTAAACGTTCGTGTACGTTTGAGAATGGCCGTTTTACCattaggtgttaatg                                                                                    <  2:348908‑M2/97‑15 (MQ=255)
             atatTTAAAGTTAAAGTGATAAAAACTTTGATATTACAGCGTTTTTAAACGTTCGTGTACGTTTGAGAATGGCCGTTTTACCattaggtgttaatgt                                                                                   >  2:200042‑M2/1‑82 (MQ=255)
                tttAAAGTTAAAGTGATAAAAACTTTGATATTACAGCGTTTTTAAACGTTCGTGTACGTTTGAGAATGGCCGTTTTACCattaggtgttaatgtgtc                                                                                >  2:250193‑M2/1‑79 (MQ=255)
                tttAAAGTTAAAGTGATAAAAACTTTGATATTACAGCGTTTTTAAACGTTCGTGTACGTTTGAGAATGGCCGTTTTACCattaggtgttaatgtgtc                                                                                <  2:74834‑M2/97‑19 (MQ=255)
                  taaagttaaagtGATAAAAACTTTGATATTACAGCGTTTTTAAACGTTCGTGTACGTTTGAGAATGGCCGTTTTACCattaggtgttaatgtgtcat                                                                              >  2:320617‑M2/1‑77 (MQ=255)
                         aaagtGATAAAAACTTTGATATTACAGCGTTTTTAAACGTTCGTGTACGTTTGAGAATGGCCGTTTTACCattaggtgttaatgtgtcataatgata                                                                       >  2:141938‑M2/1‑70 (MQ=255)
                             tGATAAAAACTTTGATATTACAGCGTTTTTAAACGTTCGTGTACGTTTGAGAATGGCCGTTTTACCATTAGGTGTTAATGTGTCATAATGATCAGTg                                                                   >  2:270005/1‑97 (MQ=255)
                                      cTTTGATATTACAGCGTTTTTAAACGTTCGTGTACGTTTGAGAATGGCCGTTTTACCATTAGGTGTTAATGTGTCATAATGATAAGTGATTTTATTg                                                          <  1:320021/97‑1 (MQ=255)
                                              ttACAGCGTTTTTAAACGTTCGTGTACGTTTGAGAATGGCCGTTTTACCATTAGGTGTTAATGTGTCATAATGATAAGTGATTTTATTGAATGGAGt                                                  >  2:308921/1‑97 (MQ=255)
                                              ttACAGCGTTTTTAAACGTTCGTGTACGTTTGAGAATGGCCGTTTTACCATTAGGTGTTAATGTGTCATAATGATAAGTGATTTTATTGAATGGAGt                                                  >  2:198214/1‑97 (MQ=255)
                                                 cAGCGTTTTTAAACGTTCGTGTACGTTTGAGAATGGCCGTTTTACCATTAGGTGTTAATGTGTCATAATGATAAGTGATTTTATTGAATGGAGTGGa                                               <  1:146918/97‑1 (MQ=255)
                                                         ttAAACGTTCGTGTACGTTTGAGAATGGCCGTTTTACCattaggtgttaatgtgtcataatgatacagataacgccttaattacagtgttttcaagg                                       >  2:225707‑M2/1‑38 (MQ=255)
                                                          tAAACGTTCGTGTACGTTTGAGAATGGCCGTTTTACCattaggtgttaatgtgtcataatgatacagataacgccttaattacagtgttttcaaggt                                      <  1:223502‑M2/97‑61 (MQ=255)
                                                           aaaCGTTCGTGTACGTTTGAGAATGGCCGTTTTACCattaggtgttaatgtgtcataatgatacagataacgccttaattacagtgttttcaaggtg                                     >  1:334231‑M2/1‑36 (MQ=255)
                                                                  cGTGTACGTTTGAGAATGGCCGTTTTACCATTAGGTGTTAATGTGTCATAATGATAAGTGATTTTATTGAATGGAGTGGACATTAGTGGATATTGGt                              <  1:244367/97‑1 (MQ=255)
                                                                  cGTGTACGTTTGAGAATGGCCGTTTTACCATTAGGTGTTAATGTGTCATAATGATAAGTGATTTTATTGAATGGAGTGGACATTAGTGGATATTGGt                              <  1:244366/97‑1 (MQ=255)
                                                                   gtgtACGTTTGAGAATGGCCGTTTTACCattaggtgttaatgtgtcataatgatacagataacgccttaattacagtgttttcaaggtgttcagaat                             <  1:228862‑M2/97‑70 (MQ=255)
                                                                   gtgtACGTTTGAGAATGGCCGTTTTACCATTAGGTGTTAATGTGTCATAATGATAAGTGATTTTATTGAATGGAGTGGACATTAGTGGATATTGGTa                             <  2:204341/97‑1 (MQ=255)
                                                                        cGTTTGAGAATGGCCGTTTTACCATTAGGTGTTAATGTGTCATAATGATAAGTGATTTTATTGAATGGAGTGGACATTAGTGGATATTGGTaaaaaa                        >  1:261713/1‑97 (MQ=255)
                                                                         gTTTGAGAATGGCCGTTTTACCATTAGGTGTTAATGTGTCATAATGATAAGTGATTTTATTGAATGGAGTGGACATTAGTGGATATTGGTAAAAAAc                       <  1:231153/97‑1 (MQ=255)
                                                                                aaTGGCCGTTTTACCattaggtgttaatgtgtcataatgatacagataacgccttaattacagtgttttcaaggtgttcagaatgtgttgcaaggta                >  2:34951‑M2/1‑15 (MQ=255)
                                                                                aaTGGCCGTTTTACCATTAGGTGTTAATGTGTCATAATGATAAGTGATTTTATTGAATGGAGTGGACATTAGTGGATATTGGTAAAAAACATGTAAt                <  1:51794/97‑1 (MQ=255)
                                                                                aaTGGCCGTTTTACCATTAGGTGTTAATGTGTCATAATGATAAGTGATTTTATTGAATGGAGTGGACATTAGTGGATATTGGTAAAAAACATGTAAt                >  2:204815/1‑97 (MQ=255)
                                                                                       gTTTTACCATTAGGTGTTAATGTGTCATAATGATAAGTGATTTTATTGAATGGAGTGGACATTAGTGGATATTGGTAAAAAACATGTAATTCCTaaa         >  1:150205/1‑97 (MQ=255)
                                                                                        ttttACCATTAGGTGTTAATGTGTCATAATGATAAGTGATTTTATTGAATGGAGTGGACATTAGTGGATATTGGTAAAAAACATGTAATTCCTaaaa        >  2:189693/1‑97 (MQ=255)
                                                                                              cattaggtgttaatgtgtcataatgatacagataacgccttaattacagtgttttcaaggtgttcagaatgtgttgcaaggtaaaaaagggcataaa  <  2:266598‑M2/97‑97 (MQ=255)
                                                                                              |                                                                                                
ACACGAGTGGAAAATATTTAAAGTTAAAGTGATAAAAACTTTGATATTACAGCGTTTTTAAACGTTCGTGTACGTTTGAGAATGGCCGTTTTACCATTAGGTGTTAATGTGTCATAATGATAAGTGATTTTATTGAATGGAGTGGACATTAGTGGATATTGGTAAAAAACATGTAATTCCTAAAAGTCAGT  >  BX571857/1025656‑1025846

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