New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3855060 =92 (1.050)3 (0.040) 3/244 NT 3.7% coding (1139/1185 nt) emrD multidrug efflux system protein
?NC_000913 3855100 = 70 (0.840)coding (1179/1185 nt) emrD multidrug efflux system protein
Rejected: Frequency below/above cutoff threshold.

TCAGGCGTTTTATTCTTACTCCGCACCGCATATGCGGGTGATGCGTGATCTGATAAAGCAGCATCGTTCTCCCATGGAGCTGATGACGATGCTGCGGTGACGTGCGCTTAAACGGGCTGCCCCTGATGCGACATCCGCGTCGCCAGCGGCAGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3855213‑3855060
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcagccCGTTTAAGCGCACGTCACCGCAGCATCGTCATCAGCTCCATGGGAGAACGATGCTGCTTTATCAGATCACGCATCACCCGCATATGCGGTGCGGAGTAAGAATAAAACGCCTGATAGCCCGCACAAAGCACGCTGTCTGGCGTCTCCT  >  NC_000913/3855100‑3855247
                                                                                                                                                                                                                                                                                                              
TCAGGCGTTTTATTCTTACTCCGCACCGCATATGCGGGTGATGCGTGATCTGATAAAGCAGCATCGTTCTCCCATGGAGCTGATGACGATGCTGCGGTGACGTGCGCTTAAACGGGCTGCCCCTGATGCGACATCCGCGTCGCCAGCG                                                                                                                                                            <  4:23907/148‑1
                                                                                                               ACGGGCTGCCCCTGATGCGACATCCGCGTCGCCAGCGGCAGCCCGTTTAAGCGCACGTCACCGCAGCATCGTCATCAGCTCCATGGGAGAACGATGCTGCTTTATCAGATCACGCATCACCCGCATATGCGGTGCGGAGTAAGAATAA                                             >  4:235758/1‑148
                                                                                                                   GCTGCCCCTGATGCGACATCCGCGTCGCCAGCGGCAGCCCGTTTAAGCGCACGTCACCGCAGCATCGTCATCAGCTCCATGGGAGAACGATGCTGCTTTATCAGATCACGCATCACCCGCATATGCGGTGCGGA                                                       >  1:413705/1‑134
                                                                                                                   GCTGCCCCTGATGCGACATCCGCGTCGCCAGCGGCAGCCCGTTTAAGCGCACGTCACCGCAGCATCGTCATCAGCTCCATGGGAGAACGATGCTGCTTTATCAGATCACGCATCACCCGCATATGCGGTGCGGA                                                       <  2:413705/134‑1
                                                                                                                                                     CAGCCCGTTTAAGCGCACGTCACCGCAGCATCGTCATCAGCTCCATGGGAGAACGATGCTGCTTTATCAGATCACGCATCACCCGCATATGCGGTGCGGAGTAAGAATAACACGCCTGATAGCCCGCACAAAGCACGCTGTCTGGCGTC      <  2:172080/149‑1
                                                                                                                                                       GCCCGTTTAAGCGCACGTCACCGCAGCATCGTCATCAGCTCCATGTGAGAACGATGCTGCTTTATCAGATCACGCATCACCCGCATATGCGGTGCGGAGTAAGAATAAAACGCC                                       <  6:16133/114‑1
                                                                                                                                                       GCCCGTTTAAGCGCACGTCACCGCAGCATCGTCATCAGCTCCATGGGAGAACGATGCTGCTTTATCAGATCACGCATCACCCGCATATGCGGTGCGGAGTAAGAATAAAACGCCTGATAGCCCGCACAAAGCACGCTGTCTGGCGTCTC    <  1:344245/149‑1
                                                                                                                                                       GCCCGTTTAAGCGCACGTCACCGCAGCATCGTCATCAGCTCCATGGGAGAACGATGCTGCTTTATCAGATCACGCATCACCCGCATATGCGGTGCGGAGTAAGAATAAAACGCC                                       >  5:16133/1‑114
                                                                                                                                                        CCCGTTTAAGCGCACGTCACCGCAGCATCGTCATCAGCTCCATGGGAGAACGATGCTGCTTTATCAGATCACGCATCACCCGCATATGCGGTGCGGAGTAAGAATAAAACGCCTGATAGCCCGCACAAAGCACGCTGTCTGGCGTCTCC   >  5:78464/1‑149
                                                                                                                                                         CCGTTTAAGCGCACGTCACCGCAGCATCGTCATCAGCTCCATGGGAGAACGATGCTGCTTTATTAGATCACGCATCACCCGCATATGCGGTGCGGAGTAAGAATAAAACGCCTGATAGCCCGCACAAAGCACGCTGTCTGGCGTCTCCT  >  3:101391/1‑149
                                                                                                                                                                                                                                                                                                              
TCAGGCGTTTTATTCTTACTCCGCACCGCATATGCGGGTGATGCGTGATCTGATAAAGCAGCATCGTTCTCCCATGGAGCTGATGACGATGCTGCGGTGACGTGCGCTTAAACGGGCTGCCCCTGATGCGACATCCGCGTCGCCAGCGGCAGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3855213‑3855060
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcagccCGTTTAAGCGCACGTCACCGCAGCATCGTCATCAGCTCCATGGGAGAACGATGCTGCTTTATCAGATCACGCATCACCCGCATATGCGGTGCGGAGTAAGAATAAAACGCCTGATAGCCCGCACAAAGCACGCTGTCTGGCGTCTCCT  >  NC_000913/3855100‑3855247

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