New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4513364155 (1.120)14 (0.110) 10/164 NT 8.4% coding (46/999 nt) fecC ferric citrate ABC transporter permease
?NC_000913 = 4513381 164 (1.290)coding (29/999 nt) fecC ferric citrate ABC transporter permease
Rejected: Frequency below cutoff threshold.

GCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4513259‑4513364
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gtgctgcTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGCGCCC  <  NC_000913/4513381‑4513283
                                                                                                                                                                                                             
GCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAG                                                                                                           <  1:1948127/100‑1
GCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAG                                                                                                           <  1:3341140/100‑1
 CGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGT                                                                                                          >  2:1224230/1‑100
   TTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGC                                                                                                        <  2:384295/100‑1
   TTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGC                                                                                                        <  1:3136853/100‑1
     GGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTG                                                                                                      <  1:1898356/100‑1
     GGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTG                                                                                                      >  1:3746541/1‑100
     GGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTG                                                                                                      >  1:2690606/1‑100
                            GCGGGTTGCATCTGCTCCGGGAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATACGTGCTGCTGTGGGGGCTTCCCGTTGCAGC                                                                               <  1:3673456/100‑1
                            GCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCTGTGGGGGCTTCCCGTTGCAGC                                                                               <  1:57677/100‑1
                               GGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCTGTGGGGGCTTCCCGTTGCAGCACT                                                                            <  2:3401876/100‑1
                                GTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCTGTGGGGGCTTCCCGTTGCAGCACTT                                                                           <  2:2780017/100‑1
                                GTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCTGTGGGGGCTTCCCGTTGCAGCACTT                                                                           <  2:715769/100‑1
                                GTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAACAATAATAAGTGCTGCTGTGGGGGCTTCCCGTTGCAGCACTT                                                                           <  1:2104146/100‑1
                                   GCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCTGTGGGGGCTTCCCGTTGCAGCACTTATT                                                                        <  1:3188066/100‑1
                                      gctaCTCCAGAAACAGGAATGGCCGATTAGCAAAGCAGCCTCAGCCAGAAAATAAGAATAAGTGCTGCTGTGGGGGCTTCCCGTTGCAGCACTTATTATT                                                                     <  1:3035114/96‑1
                                      TCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCTGTGGGGGCTTCCCGTTGCAGCACTTATTATT                                                                     <  1:1044277/100‑1
                                       CTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAACCAGAAAATAATAATAAGTTCTTCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTA                                                                    <  1:2709089/100‑1
                                              GGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTG                                                             <  1:3792260/100‑1
                                               GAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGG                                                            <  2:1740931/100‑1
                                                           GCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTT                                                 <  1:3171220/99‑1
                                                            CCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTT                                               <  2:2680208/100‑1
                                                                                                   GTGCTGCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCC        <  1:3077233/100‑1
                                                                                                     GCTGCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGC      >  1:431182/1‑100
                                                                                                     GCTGCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGC      >  2:88872/1‑100
                                                                                                      CTGCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGCG     >  2:909813/1‑100
                                                                                                      CTGCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGCG     <  1:2131175/100‑1
                                                                                                       TGCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGCGC    >  1:621874/1‑100
                                                                                                       TGCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGCGC    >  1:2752752/1‑100
                                                                                                        GCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGCGCC   >  2:1495623/1‑100
                                                                                                        GCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGCGCC   <  2:3110026/100‑1
                                                                                                         CTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGCGCCC  <  2:3296852/100‑1
                                                                                                         CTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGCGCCC  >  1:1950088/1‑100
                                                                                                         CTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGCGCCC  <  1:3493781/100‑1
                                                                                                                                                                                                             
GCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4513259‑4513364
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gtgctgcTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGCGCCC  <  NC_000913/4513381‑4513283

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