New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_012967 = 452354042 (0.150)9 (0.040) 8/126 16.5 11.0% coding (590/597 nt) fimE type 1 fimbria switch DNA invertase FimE
?NC_012967 = 4523778 109 (0.430)intergenic (+231/‑251) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_012967/4523463‑4523540
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑agagaagaATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTAC  <  NC_012967/4523778‑4523708
                                                                                                                                                     
TACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTNAAAAG                                                                               >  2:427100/1‑72
TACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAG                                                                               <  1:5514100/72‑1
 ACCGCAAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGA                                                                              <  1:5939434/72‑1
     CCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAG                                                                          >  2:1185851/1‑72
     CCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAG                                                                          >  2:1724322/1‑72
     CCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAAGAATCTCATAAACGAAAAATTAAAAAGAGAAG                                                                          >  1:313219/1‑72
      CAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGA                                                                         >  2:9423122/1‑72
                          CCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAAGTAAAAAGAGAAGAATGAGGAAGATGCGTCGAGC                                                     >  2:9734226/1‑72
                                    GGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACG                                           >  2:4000838/1‑72
                                                        ACGAAAAATTAAAAAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAAAGTTAGCTTTACATATAGCGGA                       >  2:2370709/1‑72
                                                          GAAAAATTAAAAAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCCTTACATATAGCGGAGG                     >  2:339458/1‑72
                                                           AAAAATTAAAAAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGT                    >  2:9839812/1‑72
                                                              AATTAAAAAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGNAACGTTAGCTTTACATATAGCGGAGGTGAT                 >  2:607864/1‑72
                                                                 TAAAAAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATCTAGCGGAGGTGATGTG              >  2:2524252/1‑72
                                                                 TAAAAAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTG              >  2:5393522/1‑72
                                                                  AAAAAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGA             >  2:2017287/1‑72
                                                                      AGAGAAGAATGAGGAAGATGGGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATT         <  2:8965665/72‑1
                                                                      AGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATT         <  1:2179723/72‑1
                                                                       cAGAAGAAGGAGGAAGAGGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTA        <  2:3423504/71‑1
                                                                         GNAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGANGTGAACTTAAT      <  1:967041/72‑1
                                                                           aaAATGAGGAAGGGGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTT    <  2:4675476/70‑1
                                                                            GAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTA   >  1:2308324/1‑72
                                                                            GAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTA   <  1:1748444/72‑1
                                                                            GAAGGAGGAGGAGGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTA   <  2:977495/72‑1
                                                                            GAAGGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTA   <  2:5127290/72‑1
                                                                             AATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTAC  <  1:1246823/72‑1
                                                                                                                                                     
TACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_012967/4523463‑4523540
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑agagaagaATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTAC  <  NC_012967/4523778‑4523708

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