New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 164825996 (0.700)11 (0.090) 7/212 NT 10.8% coding (326/408 nt) dicA Qin prophage; putative regulator for DicB
?NC_000913 = 1648265 92 (0.710)coding (332/408 nt) dicA Qin prophage; putative regulator for DicB

TTGCAATGCTCACCAACATGGATTCTATTTGGCGATGAAGACAAGCAACCAACACCACCTGTTGAGAAGCCAGTTGCCTTATCCCCCAAAGAACTAGAGCTCCTTGAGCTGTTTAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/1648108‑1648259
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑agctcgcATTTCGGCGAGCTGGGTATCCTGTTCTGATTCTGGCAGTGCATTAAACAGCTCAAGGAGCTCTAGTTCTTTGGGGGATAAGGCAACTGGCTTCTCAACAGGTGGTGTTGGTTGCTTGTCT  <  NC_000913/1648265‑1648146
                                                                                                                                                                                                                                                                                
TTGCAATGCTCACCAACATGGATTCTATTTGGCGATGAAGACAAGCAACCAACACCACCTGTTGAGAAGCCAGTTGCCTTATCCCCCAAAGAACTAGAGCTCCTTGAGCTGTTTAATGCACTGCCAGAATCAGAACAGGATACCCAGCT                                                                                                                             <  3:595678/149‑1
   CAATGCTCACCAACATGGATTCTATTTGGCGATGAAGACAAGCAACCAACACCACCTGTTGAGAAGCCAGTTGCCTTATCCCCCAAAGAACTAGAGCTCCTTGAGCTGTTTAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGC                                                                                                                          >  4:626051/1‑149
   CAATGCTCACCAACATGGATTCTATTTGGCGATGAAGACAAGCAACCAACACCACCTGTTGAGAAGCCAGTTGCCTTATCCCCCAAAGAACTAGAGCTCCTTGAGCTGTTTAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGC                                                                                                                          <  2:1069173/149‑1
             CAACATGGATTCTATTTGGCGATGAAGACAAGCAACCAACACCACCTGTTGAGAAGCCAGTTGCCTTATCCCCCAAAGAACTAGAGCTCCTTGAGCTGTTTAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGCATTTCGGCGA                                                                                                                <  4:1174331/149‑1
                                             CAACCAACGCCACCTGTTGAGAAGCCAGTTGCCTTATCCCCCAAAGAACTAGAGCTCCTTGAGCTGTTTAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGCATTTCGGCGAGCT                                                                                                             <  3:875697/120‑1
                                             CAACCAACACCACCTGTTGAGAAGCCAGTTGCCTTATCCCCCAAAGAACTAGAGCTCCTTGAGCTGTTTAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGCATTTCGGCGAGCT                                                                                                             >  4:875697/1‑120
                                                        ACCTGTTGAGAAGCCAGTTGCCTTATCCCCCAAAGAACTAGAGCTCCTTGAGCTGTTTAATGCACTGCCAGAATCAGAACAGGATACCCA                                                                                                                                >  1:923210/1‑90
                                                        ACCTGTTGAGAAGCCAGTTGCCTTATCCCCCAAAGAACTAGAGCTCCTTGAGCTGTTTAATGCACTGCCAGAATCAGAACAGGATACCCA                                                                                                                                <  2:923210/90‑1
                                                                 GAAGCCAGTTGCCTTATCCCCCAAAGAACTAGAGCTCCTTGAGCTGTTTAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGCATTTCGGCGAGCTGG                                                                                                           <  1:866410/102‑1
                                                                 GAAGCCAGTTGCCTTATCCCCCAAAGAACTAGAGCTCCTTGAGCTGTTTAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGCATTTCGGCGAGCTGG                                                                                                           >  2:866410/1‑102
                                                                                                            CTGTTTAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGCATTTCGGCGAGCTGG                                                                                                           <  1:98987/59‑1
                                                                                                            CTGTTTAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGCATTTCGGCGAGCTGG                                                                                                           >  4:180865/1‑59
                                                                                                            CTGTTTAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGCATTTCGGCGAGCTGG                                                                                                           >  2:98987/1‑59
                                                                                                            CTGTTTAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGCATTTCGGCGAGCTGG                                                                                                           <  3:180865/59‑1
                                                                                                               TTTAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGCATTTCGGCGAGCTGG                                                                                                           >  3:10302/1‑56
                                                                                                               TTTAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGCATTTCGGCGAGCTGG                                                                                                           <  4:10302/56‑1
                                                                                                                                                   CTCGCATTTCGGCGAGCTGGGTATCCTGTTCTGATTCTGGCAGTGCATTAAA                                                                           <  1:1377657/52‑1
                                                                                                                                                   CTCGCATTTCGGCGAGCTGGGTATCCTGTTCTGATTCTGGCAGTGCATTAAA                                                                           >  2:1377657/1‑52
                                                                                                                                                    TCGCATTTCGGCGAGCTGGGTATCCTGTTCTGATTCTGGCAGTGCATTAAACAGCTCAAGGAGCTCTAGTTCTTTGGGGGATAAGGCAACTGGCTTCTCAACAGGTGGTGTTGGTTGCTTGTC   >  3:1144155/1‑123
                                                                                                                                                    TCGCATTTCGGCGAGCTGGGTATCCTGTTCTGATTCTGGCAGTGCATTAAACAGCTCAAGGAGCTCTAGTTCTTTGGGGGATAAGGCAACTGGCTTCTCAACAGGTGGTGTTGGTTGCTTGTC   <  4:1144155/123‑1
                                                                                                                                                    TCGCATTTCGGCGAGCTGGGTATCCTGTTCTGATTCTGGCAGTGCATTAAACAGCTCAAGGAGCTCTAGTTCTTTGGGGGATAAGGCAACTGGCTTCTC                           >  2:331552/1‑99
                                                                                                                                                    TCGCATTTCGGCGAGCTGGGTATCCTGTTCTGATTCTGGCAGTGCATTAAACAGCTCAAGGAGCTCTAGTTCTTTGGGGGATAAGGCAACTGGCTTCTC                           <  1:331552/99‑1
                                                                                                                                                     CGCATTTCGGCGAGCTGGGTATCCTGTTCTGATTCTGGCAGTGCATTAAACAGCTCAAGGAGCTCTAGTTCTTTGGGGGATAAGGCAACTGGCTTCTCAACAGGTGGTGTT              >  2:314332/1‑111
                                                                                                                                                     CGCATTTCGGCGAGCTGGGTATCCTGTTCTGATTCTGGCAGTGCATTAAACAGCTCAAGGAGCTCTAGTTCTTTGGGGGATAAGGCAACTGGCTTCTCAACAGGTGGTGTT              <  1:314332/111‑1
                                                                                                                                                      GCATTTCGGCGAGCTGGGTATCCTGTTCTGATTCTGGCAGTGCATTAAACAGCTCAAGGAGCTCTAGTTCTTTGGGGGATAAGGCAACTGGCTTCTCAACAGGTGGTGTTGGTTGCTTGTCT  <  4:652539/122‑1
                                                                                                                                                      GCATTTCGGCGAGCTGGGTATCCTGATCTGATTCTGGCAGTGCATTAAACAGCTCAAGGAGCTCTAGTTCTTTGGGGGATAAGGCAACTGGCTTCTCAACAGGTGGTGTTGGTTGCTTGTCT  >  3:652539/1‑122
                                                                                                                                                                                                                                                                                
TTGCAATGCTCACCAACATGGATTCTATTTGGCGATGAAGACAAGCAACCAACACCACCTGTTGAGAAGCCAGTTGCCTTATCCCCCAAAGAACTAGAGCTCCTTGAGCTGTTTAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/1648108‑1648259
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑agctcgcATTTCGGCGAGCTGGGTATCCTGTTCTGATTCTGGCAGTGCATTAAACAGCTCAAGGAGCTCTAGTTCTTTGGGGGATAAGGCAACTGGCTTCTCAACAGGTGGTGTTGGTTGCTTGTCT  <  NC_000913/1648265‑1648146

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