New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4511910 =83 (0.950)4 (0.050) 4/254 NT 4.9% coding (505/957 nt) fecD ferric citrate ABC transporter permease
?NC_000913 4511933 = 77 (0.920)coding (482/957 nt) fecD ferric citrate ABC transporter permease
Rejected: Frequency below/above cutoff threshold.

GCGGCATGGCGGGGTTGATATTGCTGAAGATGCTGGCAAAGACCCACCAGCCGATGAAGCTGGCGCTCACCGGCGTGGCGCTTTCTGCATGCTGGGCCAGCCTGACGGATTATCTGATGCTCTCGCGCCCGCAGGATGTGAACAACGCCCTGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4512062‑4511910
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cctgcGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGC  >  NC_000913/4511933‑4512079
                                                                                                                                                                                                                                                                                                            
GCGGCATGGCGGGGTTGATATTGCTGAAGATGCTGGCAAAGACCCACCAGCCGATGAAGCTGGCGCTCACCGGCGTGGCGCTTTCTGCATGCTGGGCCAGCCTGACGGATTATCTGATGCTCTCGCGCCCGCAGGATGTGAACAACGCC                                                                                                                                                         <  7:397744/149‑1
                                                                                                                            GCGCCCGCAGGATGTGAACAACGCCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTG                                                                                      <  1:396800/92‑1
                                                                                                                            GCGCCCGCAGGATGTGAACAACGCCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTG                                                                                      >  2:396800/1‑92
                                                                                                                              GCCCGCAGGATGTGAACAACGCCCTGCGGGCGCGAGCGCATCAGATAATCCGTCAGGCTGGTCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTGGCCAGCATCGTCAGCAATATCAACCCCGC                           >  4:7861/1‑149
                                                                                                                               CCCGCAGGATGTGAACAACGCCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCAAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGCTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGC                           >  2:243558/1‑148
                                                                                                                                                     CTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCA    >  8:145062/1‑149
                                                                                                                                                     CTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAG      >  8:131985/1‑147
                                                                                                                                                      TGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCAAGCAG   >  1:339515/1‑149
                                                                                                                                                       GCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGC  >  2:36060/1‑149
                                                                                                                                                                                                                                                                                                            
GCGGCATGGCGGGGTTGATATTGCTGAAGATGCTGGCAAAGACCCACCAGCCGATGAAGCTGGCGCTCACCGGCGTGGCGCTTTCTGCATGCTGGGCCAGCCTGACGGATTATCTGATGCTCTCGCGCCCGCAGGATGTGAACAACGCCCTGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4512062‑4511910
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cctgcGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGC  >  NC_000913/4511933‑4512079

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