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.

CTGATCTGACCGTCGTCCTCCAGCCCGGTGTTGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGTTCCGCTGGCGT  >  NZ_CP009273/303002‑303173
                                                                                     |                                                                                      
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)
                          ggTGTTGAATATCCCCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCaa                                                          <  2:14816/90‑1 (MQ=255)
                           gtgtTGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCctgccagcgggcggactgatggcaacgctgg                                                         >  1:505409‑M2/1‑59 (MQ=255)
                           gtgtTGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCctgccagcgggcggactgatggcaacgctgg                                                         >  1:541028‑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)
                                            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)
                                                     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)
                                                                       agcaATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGt             <  2:143596/90‑1 (MQ=255)
                                                                            tGTTTTTCCCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGTTCCGc        <  1:165772/90‑1 (MQ=255)
                                                                            tGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGTTCCGc        <  1:80465/90‑1 (MQ=255)
                                                                            tGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGTTCCGc        <  2:557230/90‑1 (MQ=255)
                                                                                 ttACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGTTCCGCTGGCg   <  2:95209/90‑1 (MQ=255)
                                                                                  tACCCTGCCAGCCGACGCTGCCATGGGCGGTCAAGTTGGGATTGACGTAGCCGTCCGCCGCGCTATTTACGTTCAGCGTCCCGCTGGCGt  >  1:282588/1‑90 (MQ=255)
                                                                                     |                                                                                      
CTGATCTGACCGTCGTCCTCCAGCCCGGTGTTGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGTTCCGCTGGCGT  >  NZ_CP009273/303002‑303173

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