New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? W3110S.gb = 315707NA (NA)31 (0.880) 7/108 NT 100% coding (654/687 nt) ykgA predicted DNA‑binding transcriptional regulator
?W3110S.gb 576911 = 0 (0.000)coding (76/498 nt) ybcS predicted lysozyme

CTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACA  >  W3110S.gb/315639‑315757
                                                                    |                                                  
cTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATa                                                  <  1:1154947/71‑1 (MQ=32)
cTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATa                                                  <  1:2674640/71‑1 (MQ=32)
   aGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAg                                               <  1:54905/71‑1 (MQ=38)
                  gAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTg                                 <  1:869810/70‑1 (MQ=255)
                  gAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTg                                 <  1:1698010/70‑1 (MQ=255)
                  gAAAACAAGAACCTCCCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTg                                 <  1:905864/70‑1 (MQ=255)
                      aCAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGAtaataa                                                 >  1:937251/1‑50 (MQ=255)
                      aCAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATaa                                                 >  1:2513264/1‑50 (MQ=38)
                      aCAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATaa                                                 >  1:737745/1‑50 (MQ=38)
                      aCAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATaa                                                 >  1:2962332/1‑50 (MQ=38)
                      aCAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATaa                                                 >  1:2570387/1‑50 (MQ=38)
                      aCAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATaa                                                 >  1:2418632/1‑50 (MQ=38)
                      aCAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATaa                                                 >  1:2335003/1‑50 (MQ=38)
                      aCAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATaa                                                 >  1:1253601/1‑50 (MQ=38)
                      aCAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATaa                                                 >  1:1216642/1‑50 (MQ=38)
                      aCAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATaa                                                 >  1:1155231/1‑50 (MQ=38)
                                          tgtgTCCATATTACGTGGGTAGGATCATAAAGTCCTCTGGTCGTTg                                 <  1:3072184/46‑1 (MQ=255)
                                                  tatTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAg            <  1:2374739/59‑1 (MQ=255)
                                                  tatTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAg            <  1:2014506/59‑1 (MQ=255)
                                                  tatTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAg            <  1:2647870/59‑1 (MQ=255)
                                                  tatTACGTGGGTAGGATCATAAAGTCCTCGGGGCGTTGTCCACTCTTTTGGCATGTCAg            <  1:887216/59‑1 (MQ=255)
                                                    ttACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGt                                >  1:1337632/1‑37 (MQ=255)
                                                    ttACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCa          >  1:1106708/1‑59 (MQ=255)
                                                    ttACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCa          >  1:1314533/1‑59 (MQ=255)
                                                    ttACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCa          >  1:1206867/1‑59 (MQ=255)
                                                    ttACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCa          >  1:1146386/1‑59 (MQ=255)
                                                               ggATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTcaca  >  1:837534/1‑56 (MQ=255)
                                                               ggATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTcaca  >  1:892220/1‑56 (MQ=255)
                                                                    |                                                  
CTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACA  >  W3110S.gb/315639‑315757

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