New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 2039283 =29 (0.530)10 (0.190) 7/142 5.2 23.1% coding (412/747 nt) wbbA putative glycosyltransferase WbbA
?REL606 2039547 = 39 (0.760)coding (148/747 nt) wbbA putative glycosyltransferase WbbA
Rejected: Coverage evenness skew score above cutoff.

TTCATATCGATATTTTCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTTATAAGATTTATAAAACAATTTT  >  REL606/2039474‑2039615
                                                                         |                                                                    
ttCATATCGATATTTTCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAg                                                                    <  1:969232/76‑1 (MQ=255)
 tCATATCGATATTTTCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGt                                                                   <  1:104299/76‑1 (MQ=255)
     atCGATATTTTCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAAt                                                               >  1:469318/1‑76 (MQ=255)
     atCGATATTTTCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAAt                                                               >  1:469301/1‑76 (MQ=255)
      tCGATATTTTCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATg                                                              <  1:1267132/76‑1 (MQ=255)
       cGATATTTTCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGc                                                             <  1:1464840/76‑1 (MQ=255)
         atatTTTCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTa                                                           <  1:1628628/76‑1 (MQ=255)
               tCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGt                                                     <  1:1433845/76‑1 (MQ=255)
                 aaaGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCt                                                   <  1:1092886/76‑1 (MQ=255)
                   aGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCttt                                                 >  2:1600758/1‑76 (MQ=255)
                    gtgtAGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTta                                                >  1:875184/1‑76 (MQ=255)
                      gtTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTtata                                              <  1:692992/76‑1 (MQ=255)
                      gtTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTtata                                              <  1:783653/76‑1 (MQ=255)
                       ttGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTatat                                             <  2:882840/76‑1 (MQ=255)
                       ttGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTatat                                             <  1:1239409/76‑1 (MQ=255)
                       ttGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTatat                                             <  1:576575/76‑1 (MQ=255)
                       ttGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTatat                                             <  1:808283/76‑1 (MQ=255)
                        tGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATAtt                                            <  1:1189280/76‑1 (MQ=255)
                        tGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATAtt                                            <  1:620276/76‑1 (MQ=255)
                            ataACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGa                                        <  1:1730677/76‑1 (MQ=255)
                               aCGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTc                                     >  2:1035432/1‑76 (MQ=255)
                                        aGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGt                            <  1:1355971/76‑1 (MQ=255)
                                             gCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTta                       <  2:1295291/76‑1 (MQ=255)
                                              cAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTtat                      <  2:1642593/76‑1 (MQ=255)
                                               aaCGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTtata                     >  1:1472953/1‑76 (MQ=255)
                                               aaCGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTtata                     <  1:1199815/76‑1 (MQ=255)
                                               aaCGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTtata                     <  1:497096/76‑1 (MQ=255)
                                               aaCGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTtata                     <  1:497079/76‑1 (MQ=255)
                                               aaCGTTTAATGAGCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTtata                     <  1:877811/76‑1 (MQ=255)
                                               aaCGTTTAATGAACCAACAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTtata                     <  1:877829/76‑1 (MQ=255)
                                                   tttAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTTATAAGAt                 >  2:1655992/1‑76 (MQ=255)
                                                   tttAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTTATAAGAt                 >  2:1655993/1‑76 (MQ=255)
                                                   tttAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTTATAAGAt                 <  1:619240/76‑1 (MQ=255)
                                                   tttAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTTATAAGAt                 <  1:619238/76‑1 (MQ=255)
                                                     tAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTTATAAGAttt               <  2:997639/76‑1 (MQ=255)
                                                      aaTGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTTATAAGATTta              <  1:1624640/76‑1 (MQ=255)
                                                          aTCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTTATAAGATTTATaaa          >  1:832272/1‑76 (MQ=255)
                                                               tCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTTATAAGATTTATAAAACAAt     <  1:915278/76‑1 (MQ=255)
                                                                  aaaaTAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTTATAAGATTTATAAAACAAtttt  >  2:503429/1‑76 (MQ=255)
                                                                         |                                                                    
TTCATATCGATATTTTCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTTATAAGATTTATAAAACAATTTT  >  REL606/2039474‑2039615

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.