New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? minE = 254058451 (1.040)5 (0.110) 3/90 NT 8.1% coding (75/699 nt) rfaY lipopolysaccharide core biosynthesis protein
?minE = 2540601 68 (1.550)coding (92/699 nt) rfaY lipopolysaccharide core biosynthesis protein

GATTCAGAAGAGCAAGATCAAAGACTTGGTTGTTTTTACCGATGAAAACAATTCAAAGTACCTCAATGTATTA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  minE/2540512‑2540584
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tattaTAAGACAAGAAGTCATTTAATACATTGAGGTACTTTGAATTGTTTTCATCGGTAAAAACAACCAAGTC  <  minE/2540601‑2540534
                                                                                                                                             
GATTCAGAAGAGCAAGATCAAAGACTTGGTTGTTTTTACCGATGAAAACAATTCAAAGTACCTCAATGTAT                                                                        >  1:1203539/1‑71
  TTCAGAAGAGCAAGATCAAAGACTTGGTTGTTTTTACCGATGAAAACAATTCAAAGTACCTCAATGTAcat                                                                      >  1:2969095/1‑68
  TTCAGAAGAGCAAGATCAAAGACTTGGTTGTTTTTACCGATGAAAACAATTCAAAGTACCTCAATGTATTA                                                                      <  1:2963187/71‑1
                ATCAAAGACTTGGTTGTTTTTACCGATGAAAACAATTCAAAGTACCTCAATGTATTA                                                                      <  1:1738618/57‑1
                ATCAAAGACTTGGTTGTTTTTACCGATGAAAACAATTCAAAGTACCTCAATGTATTA                                                                      <  1:322928/57‑1
                      GACTTGGTTGTTTTTACCGATGAAAACAATTCAAAGTACCTCAATGTATTA                                                                      <  1:1630641/51‑1
                            GTTGTTTTTACCGATGAAAACAATTCAAAGTACCTCAATGTATTATAAGACAAGAAGTCATTT                                                    <  1:945876/63‑1
                                       CGATGAAAACAATTCAAAGTACCTCAATGTATTATAAGACAAGAAGTCATTTAA                                                  <  1:379125/54‑1
                                       CGATGAAAACAATTCAAAGTACCTCAATGTATTATAAGACAAGAAGTCATTTAA                                                  <  1:549064/54‑1
                                         ATGAAAACAATTCAAAGTACCTCAATGTATTATAAGACAAGAAGTCATTTAATA                                                >  1:2880981/1‑54
                                         ATGAAAACAATTCAAAGTACCTCAATGTATTATAAGACAAGAAGTCATTTAATA                                                >  1:862247/1‑54
                                                                    TATTATAAGACAAGAAGTCATTTAATACATTGAGGTACTTTGAATTGTTTTCATCGGTAAAAACAACCA      <  1:1384372/69‑1
                                                                    TATTATAAGACAAGAAGTCATTTAATACATTGAGGTACTTTGAATTGTTTTCATCGGTAAAA             <  1:1872247/62‑1
                                                                      TTATAAGACAAGAAGTCATTTAATACATTGAGGTACTTTGAATTGTTTTCATCGGTAAAAACAACCAAGTC  >  1:1553824/1‑71
                                                                        ATAAGACAAGAAGTCATTTAATACGTTGAGGTACTTTGAATT                             >  1:2775889/1‑42
                                                                        ATAAGACAAGAAGTCATTTAATACATTGAGGTACTTTGAATT                             >  1:647725/1‑42
                                                                                                                                             
GATTCAGAAGAGCAAGATCAAAGACTTGGTTGTTTTTACCGATGAAAACAATTCAAAGTACCTCAATGTATTA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  minE/2540512‑2540584
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tattaTAAGACAAGAAGTCATTTAATACATTGAGGTACTTTGAATTGTTTTCATCGGTAAAAACAACCAAGTC  <  minE/2540601‑2540534

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