New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 2507756 =NA (NA)5 (0.300) 5/162 0.9 14.7% intergenic (+33/‑54) nupC/BW25113_RS12490 nucleoside permease NupC/IS4‑like element IS421 family transposase
?NZ_CP009273 = 4533433 30 (1.760)coding (502/549 nt) fimA type 1 fimbrial major subunit FimA
Rejected: Frequency below/above cutoff threshold.

ACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCC  >  NZ_CP009273/4533349‑4533521
                                                                                    |                                                                                        
aCATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGt                                                                                     <  1:247652/90‑1 (MQ=255)
     tAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCataagcgcta                                                                                <  1:400952‑M2/90‑11 (MQ=255)
           aGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGc                                                                          >  1:352662/1‑90 (MQ=255)
            gAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCataagcgctaacttaag                                                                         >  1:229761‑M2/1‑73 (MQ=255)
            gAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCataagcgctaacttaag                                                                         <  2:229761‑M2/90‑18 (MQ=255)
            gAAACAACCCTGAATAACGGAACCAATACCATTCCGTCCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCg                                                                         >  1:291749/1‑90 (MQ=255)
                 aacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGc                                                                    >  2:53223/1‑90 (MQ=255)
                 aacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGc                                                                    >  2:382738/1‑90 (MQ=255)
                   cccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGa                                                                  <  1:327376/90‑1 (MQ=255)
                          tAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCaa                                                           <  1:191479/90‑1 (MQ=255)
                          tAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCaa                                                           <  1:97051/90‑1 (MQ=255)
                               gAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTc                                                      <  2:45334/90‑1 (MQ=255)
                                 aCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCataagcgctaacttaagggttgtggtattacgcctgat                                                    <  2:438356‑M2/90‑39 (MQ=255)
                                 aCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAg                                                    <  1:363435/90‑1 (MQ=255)
                                 aCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAg                                                    >  2:86756/1‑90 (MQ=255)
                                 aCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAg                                                    >  1:90836/1‑90 (MQ=255)
                                 aCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAg                                                    >  2:383326/1‑90 (MQ=255)
                                 aCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAg                                                    <  2:352662/90‑1 (MQ=255)
                                  ccAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCataagcgctaacttaagggttgtggtattacgcctgata                                                   >  2:8286‑M2/1‑51 (MQ=255)
                                  ccAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGt                                                   <  1:272163/90‑1 (MQ=255)
                                   cAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTa                                                  <  2:338643/90‑1 (MQ=255)
                                           ttCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAc                                          >  1:439263/1‑90 (MQ=255)
                                               gTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAc                                      >  2:445173/1‑90 (MQ=255)
                                               gTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAc                                      >  2:90037/1‑90 (MQ=255)
                                                 tCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAccc                                    <  1:10803/90‑1 (MQ=255)
                                                 tCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAccc                                    <  1:86756/90‑1 (MQ=255)
                                                  ccAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCa                                   <  2:231939/90‑1 (MQ=255)
                                                  ccAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCa                                   <  1:2897/90‑1 (MQ=255)
                                                           ttATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAggg                          <  2:90836/90‑1 (MQ=255)
                                                               tttGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGt                      <  2:291749/90‑1 (MQ=255)
                                                               tttGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGt                      <  2:196588/90‑1 (MQ=255)
                                                                              gCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGa       <  2:439263/90‑1 (MQ=255)
                                                                                  ccccGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGcc   <  2:241818/90‑1 (MQ=255)
                                                                                   cccGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGccc  <  1:283226/90‑1 (MQ=255)
                                                                                   cccGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGccc  <  1:456645/90‑1 (MQ=255)
                                                                                   ccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataactggtcagcaattc  >  2:483098‑M2/1‑2 (MQ=255)
                                                                                    |                                                                                        
ACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCC  >  NZ_CP009273/4533349‑4533521

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