New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4513364177 (1.260)7 (0.050) 6/164 NT 4.0% coding (46/999 nt) fecC ferric citrate ABC transporter permease
?NC_000913 = 4513381 176 (1.370)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                                                                                                           >  2:3422324/1‑100
 CGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGT                                                                                                          >  2:99730/1‑100
 CGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGT                                                                                                          <  2:1253552/100‑1
  GTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTG                                                                                                         >  2:1346801/1‑100
  GTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTG                                                                                                         >  1:1401865/1‑100
   TTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCCGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGC                                                                                                        >  1:1701966/1‑100
     GGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTG                                                                                                      >  2:3802550/1‑100
     GGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTG                                                                                                      >  2:2175355/1‑100
      GCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAGACAGACTCAGCCAGAAAATAATAATAAGTGCTGC                                                                                                     <  1:1116364/100‑1
                                GTTGCATCTGCTGCGGAAACAGGAATGGCCGAGTAGCAAAACATACTCAGCCAGAAAATAATAATAAGTGCTGCTGTGGGGGCTTCCCGTTGCAGCACTT                                                                           <  1:3075359/100‑1
                                 TTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCTGTGGGGGCTTCCCGTTGCAGCACTTA                                                                          <  2:422707/100‑1
                                  TGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCTGTGGGGGCTTCCCGTTGCAGCACTTAT                                                                         <  1:1149930/100‑1
                                  TGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCTGTGGGGGCTTCCCGTTGCAGCACTTAT                                                                         <  1:1263751/100‑1
                                     ATCTGCTCCGGAAACAGGAATGGCCGAGTAGCGAAGCAGACTCAGTCAGAAAATAATAATAAGTGCTGCTGTGGGGGCTTCCCGTTGCAGCACTTATTAT                                                                      <  1:2170843/100‑1
                                      taTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCTGTGGGGGCTTCCCGTTGCAGCACTTATTATT                                                                     <  1:1285037/98‑1
                                          CTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTT                                                                 <  1:150542/100‑1
                                                                                                   GTGCTGCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCC        >  1:2571199/1‑100
                                                                                                    TGCTGCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCG       <  2:3435925/100‑1
                                                                                                    TGCTGCTGTGGGGGCTTCCCGTTGCACCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCG       <  2:2719304/100‑1
                                                                                                     GCTGCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGC      >  1:3860648/1‑100
                                                                                                      CTGCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGCG     >  2:2465300/1‑100
                                                                                                      CTGCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGCG     >  2:3546387/1‑100
                                                                                                      CTGCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGCG     >  1:3476870/1‑100
                                                                                                      CTGCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGCG     <  2:3842375/100‑1
                                                                                                        GCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGCGCC   <  2:2108293/100‑1
                                                                                                        GCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGCGCC   >  2:201624/1‑100
                                                                                                        GCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGCGCC   >  2:2712828/1‑100
                                                                                                         CTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGCGCCC  <  1:3461729/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 < 22 ≤ 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.