New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4516891 =104 (0.880)5 (0.050) 3/282 12.9 4.9% coding (827/954 nt) fecR anti‑sigma transmembrane signal transducer for ferric citrate transport; periplasmic FecA‑bound ferric citrate sensor and cytoplasmic FecI ECF sigma factor activator
?NC_000913 4516925 = 98 (0.880)coding (793/954 nt) fecR anti‑sigma transmembrane signal transducer for ferric citrate transport; periplasmic FecA‑bound ferric citrate sensor and cytoplasmic FecI ECF sigma factor activator
Rejected: Position hash score below cutoff.

TGCCCCCGCGCAAAAACGCATCGTGAACGCTGGTGAAAGCCTGCAGTTCAGCGCCTCTGAGTTTGGCGCAGTGAAACCGCTGGATGACGAGAGTACAAGCTGGACGAAGGACATCCTGAGCTTCAGCGATAAACCGCTGGGTGAGGTGATAGCCACGCTAACCCGTTACCGCAACGGCGTGCTGCGCTGCGATCCCGCCGTTGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4517094‑4516891
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgccgttgcGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAG  >  NC_000913/4516925‑4517116
                                                                                                                                                                                                                                                                                                                                                                                                            
TGCCCCCGCGCAAAAACGCATCGTGAACGCTGGTGAAAGCCTGCAGTTCAGCGCCTCTGAGTTTGGCGCAGTGAAACCGCTGGATGACGAGAGTACAAGCTGGACGAAGGACATCCTGAGCTTCAGCGATAAACCGCTGGGTGAGGTGATAGCCACGCTAACCCGTTACCGCAACGGCGTGCTGCGCTGCGATCCCGCCG                                                                                                                                                                                                      <  1:138922/200‑1
                                                          GAGTTTGGCGCAGTGAAACCGCTGGATGACGAGAGTACAAGCTGGACGAAGGACATCCTGAGCTTCAGCGATAAACCGCTGGGTGAGGTGATAGCCACGCTAACCCGTTACCCCAACGGCGTGCTGCGCTGCGATCCCGCC                                                                                                                                                                                                       >  2:184124/1‑141
                                                                  CGCAGTGAAACCGCTGGATGACGAGAGTACAAGCTGGACGAAGGACATCCTGAGCTTCAGCGATAAACCGCTGGGTGAGGTGATAGCCACGCTAACCCGTTACCGCAACGGCGTGCTGCGCTGCGATCCCGC                                                                                                                                                                                                        <  3:377762/132‑1
                                                                      GTGAAACCGCTGGATGACGAGAGTACAAGCTGGACGAAGGACATCCTGAGCTTCAGCGATAAACCGCTGGGTGAGGTGATAGCCACGCTAACCCGTTCCCGCAACGGCGTGCTGCGCTGCGCTCCCGCCGTT                                                                                                                                                                                                    >  3:311321/1‑132
                                                                      GTGAAACCGCTGGATGACGAGAGTACAAGCTGGACGAAGGACATCCTGAGCTTCAGCGATAAACCGCTGGGTGAGGTGATAGCCACGCTAACCCGTTACCGCAACGGCGTGCTGCGCTGCGATCCCGCCGTT                                                                                                                                                                                                    >  3:615781/1‑132
                                                                       TGAAACCGCTGGAGGACGAGAGTACAAGCTGGACGAAGGACATCCTGAGCTTCAGCGATAAACCGCTGGGTGAGGTGATAGCCACGCTAACCCGTTACCGCAACGGCGTGCTGCGCTGCGATCCCGCCGTTG                                                                                                                                                                                                   <  3:530627/132‑1
                                                                        GAAACCGCTGGATGACGAGAGTACAAGCTGGACGAAGGACATCCTGAGCTTCAGCGATAAACCGCTGGGTGAGGTGATAGCCACGCTAACCCGTTACCGCAACGGCGTGCTGCGCTGCGATCCCGCCGTTGCGGTAACGGG                                                                                                                                                                                         <  2:285136/141‑1
                                                                                    TGAGGAGAGTACAAGCTGGACGAAGGACACCCCGAGCTTCAGCGATAAACCGCTGTGTGAGGTGATAGCCACGCTAACCCGTTACCACAACGGCGTGCTGCGCTGCGATCCCGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTC                                                                                                                  <  1:383921/200‑1
                                                                                                                        CTTCAGCGATAAACCGCTGGGTGAGGTGATAGCCACGCTAACCCGTTACCGCAACGGCGTGCTGCGCTGCGATCCCGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGA                                                                                                                                         <  2:859745/141‑1
                                                                                                                                 TAAACCGCTGGGTGGGGTGCTAGCCACGCTAACCCGTTACCGCAACGGCGTGCTGCGCTGCGATCCCGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGA                                                                                                                                         <  3:146508/132‑1
                                                                                                                                      CGCTGGGTGAGGTGATAGCCTCGCTAACCCGTTACCGCAACGGCGTGCTGCGCTGGGATCCCGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTC                                                                                                                  <  4:515453/150‑1
                                                                                                                                                                                                    GCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAG  <  1:733527/200‑1
                                                                                                                                                                                                    GCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGACTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAG  <  1:564243/200‑1
                                                                                                                                                                                                     CCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTT                                                            >  2:88161/1‑141
                                                                                                                                                                                                     CCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAAC                                                                     >  3:763957/1‑132
                                                                                                                                                                                                          GCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACC                                                       >  2:547660/1‑141
                                                                                                                                                                                                          GCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAG                                                                <  3:178857/132‑1
                                                                                                                                                                                                                                                                                                                                                                                                            
TGCCCCCGCGCAAAAACGCATCGTGAACGCTGGTGAAAGCCTGCAGTTCAGCGCCTCTGAGTTTGGCGCAGTGAAACCGCTGGATGACGAGAGTACAAGCTGGACGAAGGACATCCTGAGCTTCAGCGATAAACCGCTGGGTGAGGTGATAGCCACGCTAACCCGTTACCGCAACGGCGTGCTGCGCTGCGATCCCGCCGTTGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4517094‑4516891
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgccgttgcGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAG  >  NC_000913/4516925‑4517116

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