New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 30307421 (1.020)4 (0.210) 3/154 1.6 16.3% coding (1970/2526 nt) ecpC fimbrial usher EcpC
?NZ_CP009273 = 303087 22 (1.180)coding (1957/2526 nt) ecpC fimbrial usher EcpC
Rejected: Frequency below/above cutoff threshold.

CGTTGATTTTGGCGCTGATCTGACCGTCGTCCTCCAGCCCGGTGTTGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGT  >  NZ_CP009273/302988‑303162
                                                                                      |                                                                                        
cGTTGATTTTGGCGCTGATCTGACCGTCGTCCTCCAGCCCGGTGTTGAATATCCCCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCaa                                                                                       <  1:411267/90‑1 (MQ=255)
cGTTGATTTTGGCGCTGATCTGACCGTCGTCCTCCAGCCCGGTGTTGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCaa                                                                                       <  1:131635/90‑1 (MQ=255)
cGTTGATTTTGGCGCTGATCTGACCGTCGTCCTCCAGCCCGGTGTTGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCaa                                                                                       <  1:408026/90‑1 (MQ=255)
cGTTGATTTTGGCGCTGATCTGACCGTCGTCCTCCAGCCCGGTGTTGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCaa                                                                                       <  2:274437/90‑1 (MQ=255)
              cTGATCTGACCGTCGTCCTCCAGCCCGGTGTTGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCctgc                                                                         >  2:44005‑M2/1‑86 (MQ=255)
                   cTGACCGTCGTCCTCCAGCCCGTTGTTGAATATCCCCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGcc                                                                    <  1:434103/90‑1 (MQ=255)
                        cgtcgtCCTCCGGCCCGGTGTTGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGc                                                               <  2:334450/90‑1 (MQ=255)
                          tcgtcCTCCAGCCCGGTGTTGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTg                                                             >  1:303070/1‑90 (MQ=255)
                               ctccAGCCCGGTGTTGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCctgccagcgggcggactgatg                                                        <  2:505409‑M2/90‑22 (MQ=255)
                                    gCCCGGTGTTGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCgg                                                   >  1:480143/1‑90 (MQ=255)
                                       cGGTGTTGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGggtaaaaacattgctgccagcgggcggactgatggcaacgct                                                >  2:191610‑M1/1‑48 (MQ=255)
                                        ggTGTTGAATATCCCCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCaa                                               <  2:14816/90‑1 (MQ=255)
                                         gtgtTGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCctgccagcgggcggactgatggcaacgctgg                                              >  1:541028‑M2/1‑59 (MQ=255)
                                         gtgtTGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCctgccagcgggcggactgatggcaacgctgg                                              >  1:505409‑M2/1‑59 (MQ=255)
                                              gAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAggtggt                                         >  1:363270/1‑90 (MQ=255)
                                                   tCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGa                                    >  1:178589/1‑90 (MQ=255)
                                                      ccccAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGt                                 <  2:142962/90‑1 (MQ=255)
                                                        ccAGCGTTGCCATCAGTCCGCCCGCTGGCAGggtaaaaacattgctgccagcgggcggactgatggcaacgctggggtgatattcaacac                               >  1:541028‑M1/1‑31 (MQ=255)
                                                        ccAGCGTTGCCATCAGTCCGCCCGCTGGCAGggtaaaaacattgctgccagcgggcggactgatggcaacgctggggtgatattcaacac                               >  1:505409‑M1/1‑31 (MQ=255)
                                                          aGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCctgccagcgggcggactgatggcaacgctggggtgatattcaacaccg                             >  2:191610‑M2/1‑42 (MQ=255)
                                                          aGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGcc                             >  2:250730/1‑90 (MQ=255)
                                                              ttGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTTAAGTAGCCGGcc                         >  2:516040/1‑90 (MQ=255)
                                                                 ccATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTccgcc                      <  2:302043/90‑1 (MQ=255)
                                                                  cATCAGTCCGCCCGCTGGCAGggtaaaaacattgctgccagcgggcggactgatggcaacgctggggtgatattcaacaccgggctggag                     <  2:505409‑M1/90‑70 (MQ=255)
                                                                   aTCAGTCCGCCCGCTGGCAGCAATGTTTTTCCCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCcgccgc                    <  1:488753/90‑1 (MQ=255)
                                                                     cAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGAATTGACGTAGCCGTCCGCCgcgc                  >  2:217235/1‑90 (MQ=255)
                                                                             cggcTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATTGACg          >  2:527331/3‑90 (MQ=255)
                                                                                cTGGCAGCAATGTTTTTCCCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGACCGCCGCGCTATTGACGTTc       <  1:93502/90‑1 (MQ=255)
                                                                                   gcagcaATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGc    <  1:86431/90‑1 (MQ=255)
                                                                                   gCAGggtaaaaacattgctgccagcgggcggactgatggcaacgctggggtgatattcaacaccgggctggaggacgacggtcagatcag    >  2:44005‑M1/1‑4 (MQ=255)
                                                                                     agcaATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGt  <  2:143596/90‑1 (MQ=255)
                                                                                      |                                                                                        
CGTTGATTTTGGCGCTGATCTGACCGTCGTCCTCCAGCCCGGTGTTGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGT  >  NZ_CP009273/302988‑303162

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