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.

TTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATAAAGCGTTGCCATATCCCATCTTATATTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAA  >  REL606/2039208‑2039355
                                                                           |                                                                        
tttttAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATAAAGCGTTGCCATATCCCATCt                                                                          >  1:1469636/1‑76 (MQ=255)
 ttttAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATAAAGCGTTGCCATATCCCATCtt                                                                         <  2:1367894/76‑1 (MQ=255)
  tttAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATAAAGCGTTGCCATATCCCATCTta                                                                        <  2:346054/76‑1 (MQ=255)
   ttAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATAAAGCGTTGCCATATCCCATCTtat                                                                       <  2:897859/76‑1 (MQ=255)
       tGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATAAAGCGTTGCCATATCCCATCTTATattt                                                                   <  2:1569179/76‑1 (MQ=255)
               gATGTTCACAATATTCACTATATTTATCTGGATGATCATAAAGCGTTGCCATATCCCATCTTATATTTGATTTGAc                                                           >  2:692992/1‑76 (MQ=255)
               gATGTTCACAATATTCACTATATTTATCTGGATGATCATAAAGCGTTGCCATATCCCATCTTATATTTGATTTGAc                                                           >  2:1467176/1‑76 (MQ=255)
                    tCACAATATTCACTATATTTATCTGGATGATCATAAAGCGTTGCCATATCCCATCTTATATTTGATTTGACCAAGt                                                      <  1:1614573/76‑1 (MQ=255)
                      acaATATTCACTATATTTATCTGGATGATCATAAAGCGTTGCCATATCCCATCTTATATTTGATTTGACCAAGTTa                                                    <  2:831418/76‑1 (MQ=255)
                          tatTCACTATATTTATCTGGATGATCATAAAGCGTTGCCATATCCCATCTTATATTTGATTTGACCAAGTTAAata                                                <  1:625376/76‑1 (MQ=255)
                            ttCACTATATTTATCTGGATGATCATAAAGCGTTGCCATATCCCATCTTATATTTGATTTGACCAAGTTAAATAtc                                              >  1:437885/1‑76 (MQ=255)
                             tCACTATATTTATCTGGATGATCATAAAGCGTTGCCATATCCCATCTTATATTTGATTTGACCAAGTTAAATAtct                                             >  1:909169/1‑76 (MQ=255)
                                cTATATTTATCTGGATGATCATAAAGCGTTGCCATATCCCATCTTATATTTGATTTGACCAAGTTAAATATCtcat                                          <  1:623311/76‑1 (MQ=255)
                                cTATATTTATCTAGATGATCATAAAGCGTTGCCATATCCCATCTTATATTTGATTTGACCAAGTTAAATATCtcat                                          <  1:1163113/76‑1 (MQ=255)
                                  atatTTATCTGGATGATCATAAAGCGTTGCCATATCCCATCTTATATTTGATTTGACCAAGTTAAATATCccatca                                        <  2:1006883/76‑1 (MQ=255)
                                    atTTATCTGGATGATCATAAAGCGTTGCCATATCCCATCTTATATTTGATTTGACCAAGTTAAATATCTCATCAAt                                      >  1:881644/1‑76 (MQ=255)
                                     tttATCTGGATGATCATAAAGCGTTGCCATATCCCATCTTATATTTGATTTGACCAAGTTAAATATCTCATCAATc                                     >  1:143597/1‑76 (MQ=255)
                                           tGGATGATCATAAAGCGTTGCCATATCCCATCTTATATTTGATTTGACCAAGTTAAATATCTCATCAATCAAttgg                               <  2:270850/76‑3 (MQ=255)
                                            ggATGATCATAAAGCGTTGCCATATCCCATCTTATATTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAg                              <  1:1348386/76‑1 (MQ=255)
                                               tgatCATAAAGCGTTGCCATATCCCATCTTATATTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACt                           <  2:254788/76‑1 (MQ=255)
                                                   cATAAAGCGTTGCCATATCCCATCTTATATTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACTCATg                       <  1:1600331/76‑1 (MQ=255)
                                                      aaaGCGTTGCCATATCCCATCTTATATTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCaa                    >  1:1489125/1‑76 (MQ=255)
                                                           gTTGCCATATCCCATCTTATATTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAaataat               <  1:1467176/76‑1 (MQ=255)
                                                                  tatCCCATCTTATATTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTc        >  1:852488/1‑76 (MQ=255)
                                                                  tatCCCATCTTATATTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTc        <  1:333750/76‑1 (MQ=255)
                                                                    tCCCATCTTATATTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTa      <  1:938449/76‑1 (MQ=255)
                                                                     cccATCTTATATTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTaa     >  1:495321/1‑76 (MQ=255)
                                                                      ccATCTTATATTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTaaa    <  2:820201/76‑1 (MQ=255)
                                                                        aTCTTATATTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTaaaaa  <  1:928837/76‑1 (MQ=255)
                                                                           |                                                                        
TTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATAAAGCGTTGCCATATCCCATCTTATATTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAA  >  REL606/2039208‑2039355

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