New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3796153295 (1.120)3 (0.010) 3/182 30.3 1.0% coding (787/1074 nt) rfaK putative hexose transferase
?NC_000913 = 3796172 285 (1.100)coding (768/1074 nt) rfaK putative hexose transferase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GATATGAACAAGAGCTGTATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3796061‑3796153
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cAATACCAATACTTGAAATCGAAACACTACCATTTGATGAGTTTATTTATACCGTTGCGTTGACCAAGTATAGTGATTTTGTCATTTCTGTTG  <  NC_000913/3796172‑3796081
                                                                                                                                                                                         
GATATGAACAAGAGCTGTATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTC                                                                                              >  1:2888092/1‑93
GATATGAACAAGAGCTGTATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTC                                                                                              <  1:6491307/93‑1
GATATGAACAAGAGCTGTATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTC                                                                                              <  2:2710967/93‑1
GATATGAACAAGAGCTGTATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTC                                                                                              <  2:4503301/93‑1
                                                              TATAAATAAACTCATCAAATGGTAGTGTTTCAATACCAATACTTGAAATCGAAACACTACCATTTGATGAGTTTATTTATACCGTTGCGTTGA                                >  1:2756547/1‑93
                                                                TAAATAAACTCATCAAATGGTAGTGTTTCAATACCAATACTTGAAATCGAAACACTACCATTTGATGAGTTTATTTATACCGTTGCGTTGACC                              >  1:3103484/1‑93
                                                                                  GGTAGTGTTTCAATACCAATACTTGAAATCGAAACACTACCATTTGATGAGTTTATTTATACCGTTGCGTTGACCAAGTATAGTGATTTTGTC            <  2:3103484/93‑1
                                                                                            CAATACCAATACTTGAAATCGAAACACTACCATTTGATGAGTTTATTTATACCGTTGCGTTGACCAAGTATAGTGATTTTGTCATTTCTGTTG  >  2:1292822/1‑93
                                                                                            CAATACCAATACTTGAAATCGAAACACTACCATTTGATGAGTTTATTTATACCGTTGCGTTGACCAAGTATAGTGATTTTGTCATTTCTGTTG  >  2:482063/1‑93
                                                                                                                                                                                         
GATATGAACAAGAGCTGTATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3796061‑3796153
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cAATACCAATACTTGAAATCGAAACACTACCATTTGATGAGTTTATTTATACCGTTGCGTTGACCAAGTATAGTGATTTTGTCATTTCTGTTG  <  NC_000913/3796172‑3796081

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