New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 3454958 =44 (1.380)3 (0.100) 3/152 2.5 7.1% coding (132/510 nt) gspH type II secretion system minor pseudopilin GspH
?NZ_CP009273 3455001 = 39 (1.350)coding (175/510 nt) gspH type II secretion system minor pseudopilin GspH
Rejected: Frequency below/above cutoff threshold.

GTGCTGGCGCTGGTCGCTATCACGGCAAGCGTGGTGCTCTTCACATATGGCCGGGAGGATGTCGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAATTAGCCATCGACCGCGCCACGCTCAGTGGTCAGCCTGTCGGTATCCATTTCTCTGATTCAGCATG  >  NZ_CP009273/3454869‑3455047
                                                                                         |                                                                                         
cagcgccaggctgatggtgttttgcccggtaccataatgcgccatgctgaatcagagaaatggataccgacaggctgaccactgagcgtg                                                                                           <  2:94929‑M1/1‑1 (MQ=255)
cagcgccaggctgatggtgttttgcccggtaccataatgcgccatgctgaatcagagaaatggataccgacaggctgaccactgagcgtg                                                                                           <  2:147201‑M1/1‑1 (MQ=255)
cagcgccaggctgatggtgttttgcccggtaccataatgcgccatgctgaatcagagaaatggataccaacaggctgaccactgagcgtg                                                                                           <  2:1623‑M1/1‑1 (MQ=255)
 tgctggAGCTGGTCGCTATCACGGCAAGCGTGGTGCTCTTCACATATGGCCGGGAGGATGTCGCCAGCACTCGTGCGCGCGAAACAgcgg                                                                                          <  1:126506/90‑1 (MQ=255)
    tggCGCTGGTCGCTATCACGGCAAGCGTGGTGCTCTTCACATATGGCCGGGAGGATGTCGCCAGCACTCGTGCGCGCGAAACAGCGgcgc                                                                                       <  2:90440/90‑1 (MQ=255)
    tggCGCTGGTCGCTATCACGGCAAGCGTGGTGCTCTTCACATATGGCCGGGAGGATGTCGCCAGCACTCGTGCGCGCGAAACAGCGgcgc                                                                                       <  1:903334/90‑1 (MQ=255)
    gccaggctgatggtgttttgcccggtaccataatgcgccatgctgaatcagagaaatggataccgacaggctgaccactgagcgtGgcgc                                                                                       >  1:828522‑M1/86‑90 (MQ=255)
                     aCGGCAAGCGTGGTGCTCTTCACATATGGCCGGGAGGATATCGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACtt                                                                      >  1:527442/1‑90 (MQ=255)
                          aaGCGTGGTGCTCTTCACATATGGCCGGGAGGATGTCGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAAtt                                                                 >  1:799794/1‑90 (MQ=255)
                          aaGCGTGGTGCTCTTCACATATGGCCGGGAGGATGTCGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAAtt                                                                 >  1:604409/1‑90 (MQ=255)
                          aaGCGTGGTGCTCTTCACATATGGCCGGGAGGATGTCGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAAtt                                                                 >  2:580246/1‑90 (MQ=255)
                          aaGCGTGGTGCTCTTCACATATGGCCGGGAGGATGTCGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAAtt                                                                 >  1:474153/1‑90 (MQ=255)
                          aaGCGTGGTGCTCTTCACATATGGCCGGGAGGATGTCGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAAtt                                                                 <  2:895329/90‑1 (MQ=255)
                                ggtgCTCTTCACATATGGCCGGGAGGATGTCGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAATTAGCCAt                                                           >  1:715930/1‑90 (MQ=255)
                                    ctctTCACATATGGCCGGGAGGATGTCGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAATTAGCCATCGAc                                                       <  1:352222/90‑1 (MQ=255)
                                      gcgccatgctgaatcagagaaatggataccgacaggctgaccactgagcgtGGCGCGGTTTACCGCTGCACTTGAATTAGCCATCGACcg                                                     <  2:719626‑M1/39‑1 (MQ=255)
                                      ctTCACATATGGCCGGGAGGATGTCGCCAGCCCTAGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAATTAGCCATCGACcg                                                     <  2:137407/90‑1 (MQ=255)
                                          acaTATGGCCGGGAGGATGTCGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAATTAGCCATCGACCGCGcc                                                 <  1:439002/90‑1 (MQ=255)
                                           caTATGGCCGGGAGGATGTCGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAATTAGCCATCGACCGCGCCa                                                >  2:512251/1‑90 (MQ=255)
                                           atgctgaatcagagaaatggataccgacaggctgaccactgagcgtGGCGCGGTTTACCGCTGCACTTGAATTAGCCATCGACCGCGCCa                                                <  2:34802‑M1/44‑1 (MQ=255)
                                             tatGGCCGGGAGGATGTCGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGc                                                                             >  1:506505/1‑59 (MQ=255)
                                             tatGGCCGGGAGGATGTCGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGc                                                                             <  2:506505/59‑1 (MQ=255)
                                             tatGGCCGGGAGGATGTCGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAATTAGCCATCGACCGCGCCACg                                              <  1:817909/90‑1 (MQ=255)
                                             tatGGCCGGGAGGATGTCGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAATTAGCCATCGACCGCGCCACg                                              >  1:329633/1‑90 (MQ=255)
                                                ggCCGGGAGGATGTCGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAATTAGCCATCGACCGCGCCACGCTc                                           >  2:562977/1‑90 (MQ=255)
                                                ggCCGGGAGGATGTCGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAATTAGCCATCGACCGCGCCACGCTc                                           >  2:45693/1‑90 (MQ=255)
                                                  ccGGGAGGATGTCGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAATTAGCCATCGACCGCGCCACGCTCAg                                         >  1:147201/1‑90 (MQ=255)
                                                    gggAGGATGTCGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAATTAGCCATCGACCGCGCCACGCTCAGTg                                       >  2:835040/1‑90 (MQ=255)
                                                    gggAGGATGTCGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAATTAGCCATCGACCGCGCCACGCTCAGTg                                       >  1:835585/1‑90 (MQ=255)
                                                    gggAGGATGTCGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAATTAGCCATCGACCGCGCCACGCTCAGTg                                       >  1:894530/1‑90 (MQ=255)
                                                    gggAGGATGTCGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAATTAGCCATCGACCGCGCCACGCTCAGTg                                       >  2:24809/1‑90 (MQ=255)
                                                     ggaggaTGTCGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAAtt                                                                 <  1:150844/63‑1 (MQ=255)
                                                     ggaggaTGTCGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAAtt                                                                 >  2:150844/1‑63 (MQ=255)
                                                     ggaggaTGTCGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAATTAGCCATCGACCGCGCCACGCTCAGTgg                                      >  1:562587/1‑90 (MQ=255)
                                                     ggaggaTGTCGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAATTAGCCATCGACCGCGCCACGCTCAGTgg                                      >  2:557827/1‑90 (MQ=255)
                                                       aggaTGTCGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAATTAGCCATCGACCGCGCCACGCTCAGTGGTc                                    <  2:850176/90‑1 (MQ=255)
                                                       aggaTGTCGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAATTAGCCATCGACCGCGCCACGCTCAGTGGTc                                    <  2:716841/90‑1 (MQ=255)
                                                       aggaTGTCGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAATTAGCCATCGACCGCGCCACGCTCAGTGGTc                                    <  1:61262/90‑1 (MQ=255)
                                                       aggaTGTCGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAATTAGCCATCGACCGCGCCACGCTCAGTGGTc                                    <  1:72177/90‑1 (MQ=255)
                                                       aggaTGTCGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAATTAGCCATCGACCGCGCCACGCTCAGTGGTc                                    <  1:772326/90‑1 (MQ=255)
                                                           tGTCGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAATTAGCCATCGACCGCGCCACGCTCAGTGGTCAGcc                                >  1:429312/1‑90 (MQ=255)
                                                            gTCGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAATTAGCCATCGACCGCGCCACGCTCAGTGGTCAGCCt                               <  1:739022/90‑1 (MQ=255)
                                                              cGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAATTAGCCATCGACCGCGCCACGCTCAGTGGTCAGCCTGt                             <  1:110297/90‑1 (MQ=255)
                                                               gCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCCCTTGAATTAGCCATCGACCGCGCCACGCTCAGTGGTCAGCCTGTc                            <  2:36530/90‑1 (MQ=255)
                                                               gCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAATTAGCCATCGACCGCGCCACGCTCAGTGGTCAGCCTGTc                            <  2:799794/90‑1 (MQ=255)
                                                                   gCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAATTAGCCATCGACCGCGCCACGCTCAGTGGTCAGCCTGTCGGTa                        <  1:618985/90‑1 (MQ=255)
                                                                       aggctgaccactgagcgtGGCGCGGGTTACCGCCGCACTTGAATTAGCCACCGCCCGCCCCACGCGCAGTGGACACCCTGACGGTATCCa                    <  2:603158‑M1/72‑1 (MQ=255)
                                                                            cgcgcgAAACAGCGGCGCGGTTTACCGCTGCACTTGAATTAGCCATCGACCGCGCCACGCTCAGTGGTCAGCCTGTCGGTATCCATTtct               >  1:650022/1‑90 (MQ=255)
                                                                             gcgcgAAACAGCGGCGCGGTTTACCGCTGCACTTGAATTAGCCATCGACCGCGCCACGCTCAGTGGTCAGCCTGTCGGTATCCATTtctc              >  1:771135/1‑90 (MQ=255)
                                                                                   aaCAGCGGCGCGGTTTCCCGCTGCACTTGAATTAGCCATCGACCGCGCCACGCTCAGTGGTCAGCCTGTCGGTATCCATTTCTCTGATTc        <  1:309950/90‑1 (MQ=255)
                                                                                       gcggcgCGGTTTACCGCTGCACTTGAATTAGCCATCGACCGCGCCACGCTCAGTGGTCAGCCTGTCGGTATCCATTTCTCTGATTCAGCa    <  1:557827/90‑1 (MQ=255)
                                                                                         GGCGCGGTTTACCGCTGCACTTGAATTAGCCATCGACCGCGCCACGCTCAGTGGTCAGCCTGTCGGTATCCATTTCTCTGATTCAGCATg  >  2:506757‑M1/1‑90 (MQ=255)
                                                                                         |                                                                                         
GTGCTGGCGCTGGTCGCTATCACGGCAAGCGTGGTGCTCTTCACATATGGCCGGGAGGATGTCGCCAGCACTCGTGCGCGCGAAACAGCGGCGCGGTTTACCGCTGCACTTGAATTAGCCATCGACCGCGCCACGCTCAGTGGTCAGCCTGTCGGTATCCATTTCTCTGATTCAGCATG  >  NZ_CP009273/3454869‑3455047

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