New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 3344114 =29 (0.790)4 (0.110) 4/166 2.5 12.1% coding (2271/2337 nt) arcB aerobic respiration two‑component sensor histidine kinase ArcB
?NZ_CP009273 3576788 = NA (NA)coding (391/417 nt) BW25113_RS17825 YrhA family protein
Rejected: Frequency below/above cutoff threshold.

GCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAA  >  NZ_CP009273/3344028‑3344203
                                                                                      |                                                                                         
gcACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCAtc                                                                                        <  2:754470/90‑1 (MQ=255)
gcACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCAtc                                                                                        <  1:210700/90‑1 (MQ=255)
gcACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCAtc                                                                                        <  2:189302/90‑1 (MQ=255)
  aCCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCAtctc                                                                                      >  1:736428/1‑90 (MQ=255)
             gCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCAtt                                                                           >  2:992607/1‑90 (MQ=255)
             gCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCAtt                                                                           >  1:1022166/1‑90 (MQ=255)
                  ggccATTTTTTAGTGGCTTTTGCCCCCCCCGCTTTCAGCACTTCTGCGTCGTGACGCCACTCTTCTTTGATCTCTTCCATCCATTCACCg                                                                      >  1:129712/4‑90 (MQ=255)
                           ttAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCAtc                                                                                        >  2:551723/1‑63 (MQ=255)
                           ttAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCAtc                                                                                        <  1:551723/63‑1 (MQ=255)
                            cactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCATCTCTTCAATCCATTCACCGACGTTATCtt                                                            <  2:752675‑M1/32‑1 (MQ=255)
                                       gCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGa                                                 <  1:944464/90‑1 (MQ=255)
                                                  tttCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGgaga                                      >  2:697690/1‑90 (MQ=255)
                                                         aCTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAAtt                               >  1:969236/1‑90 (MQ=255)
                                                             aaatcagtaagttggcagcatcaccCATCTCTTCAATCCATTCACCGCCGTTATCTTCCCAGGCCGGAAGGTCAGGGGTCTGAATTTGCt                           <  2:707813‑M1/65‑1 (MQ=255)
                                                               aCGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGa                                                 <  1:1025226/66‑1 (MQ=255)
                                                               aCGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGa                                                 >  2:1025226/1‑66 (MQ=255)
                                                               aCGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGa                                                 <  1:196319/66‑1 (MQ=255)
                                                               aCGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGa                                                 >  2:196319/1‑66 (MQ=255)
                                                               aCGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGa                         >  2:185114/1‑90 (MQ=255)
                                                                   cgtGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCa                     >  1:1064930/1‑90 (MQ=255)
                                                                   cgtGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCa                     <  1:1077587/90‑1 (MQ=255)
                                                                    gtGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAg                    >  1:891622/1‑90 (MQ=255)
                                                                    gtGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAg                    >  2:696864/1‑90 (MQ=255)
                                                                    gtGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAg                    >  1:634867/1‑90 (MQ=255)
                                                                     tGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGc                   <  2:634867/90‑1 (MQ=255)
                                                                     tGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGc                   <  2:1064930/90‑1 (MQ=255)
                                                                      gACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCt                  >  2:335580/1‑90 (MQ=255)
                                                                              gcatcaccCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGa          <  1:437548‑M1/82‑1 (MQ=255)
                                                                               cttcttTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGAt         <  1:298101/90‑1 (MQ=255)
                                                                               cttcttTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGAt         <  1:335580/90‑1 (MQ=255)
                                                                               cttcttTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGAt         <  2:736428/90‑1 (MQ=255)
                                                                               cttcttTCATCTCTTCAATCCATCCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGAt         <  1:445464/90‑1 (MQ=255)
                                                                                 tcaccCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATgg       <  1:963429‑M1/85‑1 (MQ=255)
                                                                                      cATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTaa  <  2:891622/90‑1 (MQ=255)
                                                                                      |                                                                                         
GCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAA  >  NZ_CP009273/3344028‑3344203

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