New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 2039335 =20 (0.900)3 (0.150) 3/138 2.6 13.2% coding (360/747 nt) wbbA putative glycosyltransferase WbbA
?REL606 2039563 = 21 (1.030)coding (132/747 nt) wbbA putative glycosyltransferase WbbA
Rejected: Frequency below/above cutoff threshold.

ATAAAGCGTTGCCATATCCCATCTTATATTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTAAATCATTTTTATCAATTTTTCCAACTCGACGATCTTCATCTACTAG  >  REL606/2039260‑2039409
                                                                           |                                                                          
cggcctgagtggtttagtcatgaaaattgttttataaatcttataaatacaataaaagaatctaaatataaagacc                                                                            >  2:539064‑M1/76‑76 (MQ=255)
  gccagagtggtttagtcatcaaaaatgttttataaatcttataaatacaataaaagaatctaaatataaagaccac                                                                          >  2:655857‑M1/74‑75 (MQ=255)
   cctgagtggtttagtcatggaaattgttttataaatcttataaatacaataaaagaatctaaatataaagacCaaa                                                                         <  1:435767‑M1/4‑1 (MQ=255)
   cctgagtggtttagtcatgaaaattgttttataaatcttataaatacaataaaagaatctaaatataaagacCaaa                                                                         <  1:435769‑M1/4‑1 (MQ=255)
   aaGCGTTGCCATATCCCATCTTATATTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCaaa                                                                         <  2:84462/76‑1 (MQ=255)
       gTTGCCATATCCCATCTTATATTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAaataat                                                                     >  1:180837/1‑76 (MQ=255)
          gCCATATCCCATCTTATATTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCAt                                                                  <  1:218596/76‑1 (MQ=255)
          gCCATATCCCATCTTATATTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCAt                                                                  <  1:218605/76‑1 (MQ=255)
              tatCCCATCTTATATTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTc                                                              >  1:697298/1‑76 (MQ=255)
                  ccATCTTATATTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTaaa                                                          <  1:25174/76‑1 (MQ=255)
                     tCTTATATTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTaaaaaa                                                       <  1:218786/76‑1 (MQ=255)
                      cTTATATTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAat                                                      >  1:187542/1‑76 (MQ=255)
                      cTTATATTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAat                                                      >  1:230757/1‑76 (MQ=255)
                      cTTATATTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAat                                                      >  1:230747/1‑76 (MQ=255)
                             ttgatttgaCCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTaa                                               <  1:23771/76‑1 (MQ=255)
                               gatttgaCCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTAAAt                                             <  1:331793/76‑1 (MQ=255)
                                   tgaCCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTAAATCAtt                                         <  1:694275/76‑1 (MQ=255)
                                   tgaCCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTAAATCAtt                                         >  2:372637/1‑76 (MQ=255)
                                   aaatcttataaatacaataaaagaatctaaatataaagacCAAATAATCATTTTCTAAAAAATATATTAAATCAtt                                         <  2:460494‑M1/36‑1 (MQ=255)
                                     aCCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTAAATCAtttt                                       <  1:361375/76‑1 (MQ=255)
                                        aaGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTAAATCATTTTTAt                                    <  1:543262/76‑1 (MQ=255)
                                        aaGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTAAATCATTTTTAt                                    <  2:537607/76‑1 (MQ=255)
                                               tacaataaaagaatctaaatataaagacCAAATAATCATTTTCTAAAAAATATATTAAATCATTTTTATCAAtttt                             <  2:355284‑M1/48‑1 (MQ=255)
                                                  aataaaagaatctaaatataaagacCAAATAATCATTTTCTAAAAAATATATTAAATCATTTTTATCAATTTTTcc                          <  2:356797‑M1/51‑1 (MQ=255)
                                                   ctcATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTAAATCATTTTTATCAATTTTTCCa                         >  1:443304/1‑76 (MQ=255)
                                                   ctcATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTAAATCATTTTTATCAATTTTTCCa                         >  1:443308/1‑76 (MQ=255)
                                                                          gCAAATAATCATTTTCTAAAAAATATATTAAATCATTTTTATCAATTTTTCCAACTCGACGATCTTCATCTACTAg  >  1:629204/1‑76 (MQ=255)
                                                                           |                                                                          
ATAAAGCGTTGCCATATCCCATCTTATATTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTAAATCATTTTTATCAATTTTTCCAACTCGACGATCTTCATCTACTAG  >  REL606/2039260‑2039409

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