New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4511992159 (1.230)7 (0.060) 6/170 NT 4.3% coding (423/957 nt) fecD ferric citrate ABC transporter permease
?NC_000913 = 4512009 163 (1.320)coding (406/957 nt) fecD ferric citrate ABC transporter permease
Rejected: Frequency below cutoff threshold.

AGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCAC‑‑GCCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4511892‑4511992
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gccgATGAAGCTGGCGCTCACCGGCGTGGCGCTTTCTGCATGCTGGGCCAGCCTGACGGATTATCTGATGCTCTCGCGCCCGCAGGATGTGAACAACGCCCTG  <  NC_000913/4512009‑4511911
                                                                                                                                                                                                          
AGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCAC‑‑GCC                                                                                                      >  1:1633216/1‑100
 GCTGCCGGTCAGCGACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCAC‑‑GCC                                                                                                      >  2:372137/1‑99
                       GCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCAC‑‑GCCGATGAAGCTGGCGCTCACCGGCG                                                                               <  2:3082417/100‑1
                       GCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCAC‑‑GCCGATGAAGCTGGCGCTCACCGGCG                                                                               <  1:826866/100‑1
                           TGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCAC‑‑GCCGATGAAGCTGGCGCTCACCGGCGTGGC                                                                           <  1:142029/100‑1
                            GTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCAC‑‑GCCGATGAAGCTGGCGCTCACCGGCGTGGCG                                                                          <  1:796914/100‑1
                             TTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCAC‑‑GCCGATGAAGCTGGCGCTCACCGGCGTGGCGC                                                                         <  1:1564939/100‑1
                                          GGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGACCCACCAGCCGATGAAGCTGGCGCTCACCGGCGTGGCGCTTTCTGCAT                                                                <  2:983363/98‑1
                                            CGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCAC‑‑GCCGATGAAGCTGGCGCTCACCGGCGTGGCGCTTTCTGCATGCTGGG                                                          <  1:1001304/100‑1
                                                                                                    CCGATGAAGCTGGCGCTCACCGGCGTGGCGCTTTCTGCATGCTGGGTCAGCCTGACGGATTATCTGATGCTCTCGCGCCCGCAGGATGTGAACAACGCCC    <  1:306491/100‑1
                                                                                                    CCGATGAAGCTGGCGCTCACCGGCGTGGCGCTTTCTGCATGCTGGGCCAGCCTGACGGATTATCTGATGCTCTCGCGCCCGCAGGATGTGAACAACGCac    >  2:2206388/1‑98
                                                                                                    CCGATGAAGCTGGCGCTCACCGGCGTGGCGCTTTCTGCATGCTGGGCCAGCCTGACGGATTATCTGATGCTCTCGCGCCCGCAGGATGTGAACAACGCCC    <  2:1686153/100‑1
                                                                                                     CGATGAAGCTGGCGCTCACCGGCGTGGCGCTTTCTGCATGCTGGGCCAGCCTGACGGATTATCTGATGCTCTCGCGCCCGCCGGATGTGAACAACcc      >  2:46730/1‑95
                                                                                                     CGATGAAGCTGGCGCTCACCGGCGTGGCGCTTTCTGCATGCTGGGCCAGCCTGACGGATTATCTGATGCTCTCGCGCCCGCAGGATGTGAACAACGCCCT   >  2:619978/1‑100
                                                                                                      GATGAAGCTGGCGCTCACCGGCGTGGCGCTTTCTGCATGCTGGGCCAGCCTGACGGATTATCTGATGCTCTCGCGCCCGCAGGATGTGAACAACGCCCTG  >  2:1271529/1‑100
                                                                                                      GATGAAGCTGGCGCTCACCGGCGTGGCGCTTTCTGCATGCTGGGCCAGCCTGACGGATTATCTGATGCTCTCGCGCCCGCAGGATGTGAACAACGCCCTG  >  2:2665760/1‑100
                                                                                                      GATGAAGCTGGCGCTCACCGGCGTGGCGCTTTCTGCATGCTGGGCCAGCCTGACGGATTATCTGATGCTCTCGCGCCCGCAGGATGTGAACAACGCCCTG  >  1:2899409/1‑100
                                                                                                      GATGAAGCTGGCGCTCACCGGCGTGGCGCTTTCTGCATGCTGGGCCAGCCTGACGGATTATCTGATGCTCTCGCGCCCGCAGGATGTGAACAACGCCCTG  <  1:2034378/100‑1
                                                                                                                                                                                                          
AGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCAC‑‑GCCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4511892‑4511992
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gccgATGAAGCTGGCGCTCACCGGCGTGGCGCTTTCTGCATGCTGGGCCAGCCTGACGGATTATCTGATGCTCTCGCGCCCGCAGGATGTGAACAACGCCCTG  <  NC_000913/4512009‑4511911

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