New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 450564718 (0.600)3 (0.100) 3/162 2.4 14.5% coding (848/2325 nt) fecA Fe(3+) dicitrate transport protein FecA
?NZ_CP009273 4505651 = 18 (0.620)coding (844/2325 nt) fecA Fe(3+) dicitrate transport protein FecA
Rejected: Frequency below/above cutoff threshold.

GTGGCGCTGTGCTCGCGCCAGTCACTGCCGCGCGTGCCGGAGTAGAGCAGCGCGGTGCCAAAACCGTTGTCCGCTGTGCCGCCCACCATCAGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGC  >  NZ_CP009273/4505563‑4505736
                                                                                    |                                                                                         
gTGGCGCTGTGCTCGCGCCAGTCACTGCCGCGCGTGCCGGAGTAGAGCAGCGCGGTGCCAAAACCGTTGTCCGCTGTGCCGCCCACCATc                                                                                      >  1:743586/1‑90 (MQ=255)
 tGGCGCTGTGCTCGCGCCAGTCACTGCCGCGCGTGCCGGAGTAGAGCAGCGCGGTGCCAAAACCGTTGTCCGCTGTGCCGCCCACCATCa                                                                                     >  1:520573/1‑90 (MQ=255)
                  cAGTCACTGCCGCGCGTGCCGGAGTAGAGCAGCGCGGTGCCAAAACCGTTGTCCGCTGTGCCGCCCACCATCAGGTTGTGCGTCTCTTTc                                                                    >  1:728000/1‑90 (MQ=255)
                   aGTCACTGCCGCGCGTGCCGGAGTAGAGCAGCGCGGTGCCAAAACCGTTGTCCGCTGTGCCGCCCACCATCAGGTTGTGCGTCTCTTTCg                                                                   <  2:673316/90‑1 (MQ=255)
                        cTGCCGCGCGTGCCGGAGTAGAGCAGCGCGGTGCCAAAACCGTTGTCCGCTGTGCCGCCCACCATCAGGTTGTGCGTCTCTTTCGGgttg                                                              <  1:291454/90‑1 (MQ=255)
                                 gTGCCGGAGTAGAGCAGCGCGGTGCCAAAACCGTTGTCCGCTGTGCCGCCCACCATCAGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGaa                                                     <  2:271026/90‑1 (MQ=255)
                                      ggAGTAGAGCAGCGCGGTGCCAAAACCGTTGTCCGCTGTGCCGCCCACCATCAGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGt                                                <  1:867484/90‑1 (MQ=255)
                                      ggAGTAGAGCAGCGCGGTGCCAAAACCGTTGTCCGCTGTGCCGCCCACCATCAGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGt                                                <  2:633546/90‑1 (MQ=255)
                                           agagCAGCGCGGTGCCAAAACCGTTGTCCGCTGTGCCGCCCACCATCAGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGgc                                           <  1:61653/90‑1 (MQ=255)
                                           agagCAGCGCGGTGCCAAAACCGTTGTCCGCTGTGCCGCCCACCATCAGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGgc                                           <  2:856787/90‑1 (MQ=255)
                                                 gcgcGGTGCCAAAACCGTTGTCCGCTGTGCCGCCCACCATCAGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGgctgagct                                     >  2:820347/1‑90 (MQ=255)
                                                    cgGTGCCAAAACCGTTGTCCGCTGTGCCGCCCCCCATCAGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGAc                                  <  2:603469/90‑1 (MQ=255)
                                                    cgGTGCCAAAACCGTTGTCCGCTGTGCCGCCCACCATCAGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGAc                                  <  2:600693/90‑1 (MQ=255)
                                                           aaaaCCGTTGTCCGCTGTGCCGCCCACCATCAGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCAc                           >  2:334401/1‑90 (MQ=255)
                                                               ccGTTGTCCGCTGTGCCGCCCACCATCAGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCAcgccc                       <  2:156453/90‑1 (MQ=255)
                                                                cGTTGTCCGCTGTGCCGCCCAtcaggttgtgcgtctctttcgggttgttttgtgaagaggttgggctgagctgac                                     >  1:835655‑M1/1‑21 (MQ=255)
                                                                cGTTGTCCGCTGTGCCGCCCAtcaggttgtgcgtctctttcgggttgttttgtgaagaggttgggctgagctgac                                     <  2:835655‑M1/75‑55 (MQ=255)
                                                                       cGCTGTGCCGCCCACCATCAGGTTGTGCGTCTTTTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGAt               <  1:359291/90‑1 (MQ=255)
                                                                                 cccAtcaggttgtgcgtctctttcgggttgttttgtgaagaggttgggctgagctgaccttccacgcccgcctcgataccaaagtcctgc     <  1:135954‑M1/90‑87 (MQ=255)
                                                                                  ccAtcaggttgtgcgtctctttcgggttgttttgtgaagaggttgggctgagctgaccttccacgcccgcctcgataccaaagtcctgcg    >  2:616246‑M1/1‑3 (MQ=255)
                                                                                   caccaTCAGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTg   <  1:613242/90‑1 (MQ=255)
                                                                                    accaTCAGGTTGTGCGTCTCTTTCGGTTTTTTTTGTGAAGGGGTTGGGTTGAGCTGACCTCCCCCGCCCGCCTCGATACCAAAGTCCTGc  <  2:288536/90‑1 (MQ=255)
                                                                                    |                                                                                         
GTGGCGCTGTGCTCGCGCCAGTCACTGCCGCGCGTGCCGGAGTAGAGCAGCGCGGTGCCAAAACCGTTGTCCGCTGTGCCGCCCACCATCAGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGC  >  NZ_CP009273/4505563‑4505736

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