New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 204106122 (0.340)5 (0.080) 3/140 7.2 22.0% coding (323/1113 nt) vioA VioA, involved in dTDP‑N‑acetylviosamine synthesis
?REL606 = 2041218 15 (0.250)coding (166/1113 nt) vioA VioA, involved in dTDP‑N‑acetylviosamine synthesis
Rejected: Coverage evenness skew score above cutoff.

TAGCAATGGACAGGTAAAATTGCTGAAGTTTTTGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGTCCATTCCATAAAATAGCATGTGAA  >  REL606/2040986‑2041120
                                                                           |                                                           
tAGCAATGGACAGGTAAAATTGCTGAAGTTTTTGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATc                                                             >  2:1139069/1‑76 (MQ=255)
                 aaTTGCTGAAGTTTTTGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATca                                            <  1:1286691/76‑1 (MQ=255)
                 aaTTGCTGAAGTTTTTGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATca                                            <  2:1028941/76‑1 (MQ=255)
                      cTGAAGTTTTTGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACaaa                                       <  1:1902131/76‑1 (MQ=255)
                        gAAGTTTTTGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGa                                     >  1:1334449/1‑76 (MQ=255)
                        gAAGTTTTTGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGa                                     >  1:1928966/1‑76 (MQ=255)
                        gAAGTTTTTGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGa                                     >  1:919497/1‑76 (MQ=255)
                         aaGTTTTTGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGAc                                    <  1:810086/76‑1 (MQ=255)
                              tttGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCg                               >  1:1655111/1‑76 (MQ=255)
                                 ggTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTaa                            <  2:1603840/76‑1 (MQ=255)
                                     tAATCGCCCGTTCCCTTTTTCTGTACTCAATATTATATCCCTCATTTTCACTATCACCAAAGACCGGCGTAAGtcc                        >  1:467125/1‑76 (MQ=255)
                                     tAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGtcc                        >  1:467104/1‑76 (MQ=255)
                                          gCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGtccattcc                   <  1:1892367/76‑1 (MQ=255)
                                             tGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGTCCATTCCATa                >  2:367287/1‑76 (MQ=255)
                                               ttCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGTCCATTCCATaaa              <  1:1972501/76‑1 (MQ=255)
                                               ttCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGTCCATTCCATaaa              <  1:1972478/76‑1 (MQ=255)
                                                tCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGTCCATTCCATaaaa             >  1:87737/1‑76 (MQ=255)
                                                 cTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGTCCATTCCATAAAAt            >  1:1007532/1‑76 (MQ=255)
                                                  tATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGTCCATTCCATAAAATa           >  1:161386/1‑76 (MQ=255)
                                                     ttttCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGTCCATTCCATAAAATAGCa        <  1:1446905/76‑1 (MQ=255)
                                                      tttCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGTCCATTCCATAAAATAGCAt       >  1:281776/1‑76 (MQ=255)
                                                          tGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGTCCATTCCATAAAATAGCATGTGa   <  1:709010/76‑1 (MQ=255)
                                                           gTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGTCCATTCCATAAAATAGCATGTGaa  >  1:1638960/1‑76 (MQ=255)
                                                                           |                                                           
TAGCAATGGACAGGTAAAATTGCTGAAGTTTTTGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGTCCATTCCATAAAATAGCATGTGAA  >  REL606/2040986‑2041120

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