New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 2039188 =35 (0.590)11 (0.200) 9/140 4.8 14.1% coding (507/747 nt) wbbA putative glycosyltransferase WbbA
?REL606 2039561 = 102 (1.840)coding (134/747 nt) wbbA putative glycosyltransferase WbbA
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

ACCCTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATAA  >  REL606/2039114‑2039263
                                                                          |                                                                           
aCCCTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACaa                                                                            <  1:43796/76‑1 (MQ=255)
aCCCTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACaa                                                                            <  1:1609730/76‑1 (MQ=255)
aCCCTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACaa                                                                            <  1:1770230/76‑1 (MQ=255)
aCCCTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACaa                                                                            <  1:1770251/76‑1 (MQ=255)
   cTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACAATAg                                                                         <  1:871389/76‑1 (MQ=255)
    tGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACAATAGt                                                                        >  1:1571766/1‑76 (MQ=255)
      gCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACAATAGttt                                                                      <  1:1451480/76‑1 (MQ=255)
      gCTTGCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCGACTCCCAGAATAAAATACAATAGttt                                                                      <  1:563668/76‑1 (MQ=255)
               atatatGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACAATAGTTTTTTTGtttt                                                             >  1:1619002/1‑76 (MQ=255)
                    tGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTAttt                                                        <  1:1825619/76‑1 (MQ=255)
                       tGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTAtttttt                                                     <  1:1296935/76‑1 (MQ=255)
                       tGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTAtttttt                                                     >  2:491656/1‑76 (MQ=255)
                        gCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTa                                                    <  1:621277/76‑1 (MQ=255)
                            aTGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTAATGa                                                >  1:431731/1‑76 (MQ=255)
                                 cTTGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTTATGAATTaa                                           >  1:1090391/1‑76 (MQ=255)
                                      atttATTTTCCAATGCCTACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGt                                      >  2:1438066/4‑76 (MQ=255)
                                          aTTTTCCAATGCCTACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTcac                                  >  1:538664/1‑76 (MQ=255)
                                                gaaTGCCTACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATAtt                            >  2:1535523/2‑76 (MQ=255)
                                                    gCCTACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACt                        >  1:812519/1‑76 (MQ=255)
                                                       tACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACtata                     <  1:340765/76‑1 (MQ=255)
                                                       tACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACtata                     <  1:1258397/76‑1 (MQ=255)
                                                         cTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATAtt                   >  1:289080/1‑76 (MQ=255)
                                                          tCCCAGAATAAAATACAATAGTTTTTTTGTTTTTTTTTTTTAATAAATCAAGATGTTCACAATATTCACTATAttt                  >  1:1255821/1‑76 (MQ=255)
                                                            ccAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTAt                <  1:1333378/76‑1 (MQ=255)
                                                               gAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTg             <  2:567327/76‑1 (MQ=255)
                                                                   aaaaTACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGgatg         <  1:456583/76‑1 (MQ=255)
                                                                   aaaaTACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGgatg         >  2:1146992/1‑76 (MQ=255)
                                                                   aaaaTACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGgatg         >  2:1146999/1‑76 (MQ=255)
                                                                   aaaaTACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGgatg         >  2:1603793/1‑76 (MQ=255)
                                                                     aaTACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGgatgat       <  2:1811570/76‑1 (MQ=255)
                                                                      aTACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATc      >  1:598810/1‑76 (MQ=255)
                                                                      aTACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATc      >  1:598796/1‑76 (MQ=255)
                                                                       tACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCa     >  1:222304/1‑76 (MQ=255)
                                                                       tACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCa     <  1:1780349/76‑1 (MQ=255)
                                                                        caaaTAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCAt    >  2:361676/3‑76 (MQ=255)
                                                                        aCAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCAt    <  1:439761/76‑1 (MQ=255)
                                                                          aaTAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATaa  >  2:354923/1‑76 (MQ=255)
                                                                          aaTAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATaa  <  1:1479991/76‑1 (MQ=255)
                                                                          |                                                                           
ACCCTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATAA  >  REL606/2039114‑2039263

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