New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 1397238NA (NA)8 (0.800) 6/322 0.3 100% noncoding (1195/1195 nt) IS5 repeat region
?NC_000913 = 1978494 0 (0.000)intergenic (‑297/+24) flhD/insB5 DNA‑binding transcriptional dual regulator FlhD/IS1 family protein InsB

CCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCG  >  NC_000913/1978323‑1978723
                                                                                                                                                                           |                                                                                                                                                                                                                                     
ccccctttTTCGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACGACACGCACATACAACGGAGGGGGGCTGCGATTGTCAATAATGCGTGATGCAGATCACACAAAACACTCAAGTACTTAAggaaggtgcgaataagcggggaaattcttcgcggctgactcagtcatttcagttcttcatgtttgagccgattttttct                                                                                                                                                        >  1:96640‑M2/4‑172 (MQ=255)
             ttAATATGTCCTTACAAATAGAAATGGGTCTTGACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGTATGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaatgtgcga                                                                                                                                                                                                                            >  1:7009‑M2/1‑159 (MQ=255)
                                 gAAATGGGGCTTGACACTTATCTAAGATTTTTCCTAAAGCGACGCAACGGGACTCGTCACGACACGCACATACAACGGAGGGGGGCTGCGATTTGCAAGAATGCGTGAGGCAGATCACACAAAACACGCAATTACTTAAggaaggtgcgaataagcggggaaattcttctcggctn                                                                                                                                                                                                  >  1:24918‑M2/1‑139 (MQ=255)
                                       ngtCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcagcctatttagaccgtttct                                                                                                                 <  1:129975‑M2/250‑119 (MQ=255)
                                                                tCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGAGCACACAAAACACTCAAGTACTTAAggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgcc                                                                                                                                            >  1:156033‑M2/1‑108 (MQ=255)
                                                                                                         cAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcagcctatttagaccgtttcttcgccatttaaggcgttatccccagtttttagtgagatctctcccactgacgtatcatttggtccg                                               <  1:197395‑M2/251‑185 (MQ=255)
                                                                                                                                         gTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaataagcgggtaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcagcctatttagaccgtttctgcgccatttaaggcgttatccccagtttttagtgagatctctcccactgacgtatcatttgttccgcccgaaacaggttggccagcgtgaataacatc               >  1:67379‑M2/1‑35 (MQ=255)
                                                                                                                                                      cacaAAACACTCAATTACTTAAggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcagcctatttagaccgtttcttcgccatttaaggcgttatcccagtttttagtgagatctctcccgctgacgtatcatttggtccgcccgaaacaggttggccagcgtgaataacagcgccagtttgttatc  >  1:76934‑M2/1‑22 (MQ=255)
                                                                                                                                                                           |                                                                                                                                                                                                                                     
CCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCG  >  NC_000913/1978323‑1978723

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