New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 33520234 (0.190)6 (0.290) 6/262 NT 60.0% coding (1003/2337 nt) arcB sensory histidine kinase ArcB
?NC_000913 3365556 = NA (NA)noncoding (1195/1195 nt) IS5 repeat region

CAATCCTTTTTGTTGCGCCTGCAATGCGGAGAGATTTTCCAGATCGGCAAGGAAGCTGGTGAAATCAACCGGTTGATTATCAAGCTGGACCTTGCGCCGTTCCATCTTATCCATGTCGATAATATCGTTAAAGATATTCCCCAGCGTGACGGCCGAAACATGGATGGTCTTGAGATATTTTTCCTGCTCGGCGGTGAGTTCGGTATCCAGCAGAATGCGGCTCAGACCGACGATACCGTTCAGCGGTGTACGCAATTCGTGACTGATGGT  >  NC_000913/3351892‑3352161
                                                                                                                                   |                                                                                                                                          
cAATCCTTTTGGTTGCGCCTGCAATGCGGAGAGATTTTCCAGATCGGCAAGGAAGCTGGTGAAATCAACCGGTTGATTATCAAGCTGGACCTTGCGCCGTTCCATCTTATCCATGTCGATAATATCGTTAAAGATATTCCCCa                                                                                                                                 <  1:2827/143‑1 (MQ=255)
                            gagagaTTTTCCAGATCGGCAAGTAAGCTGGTGAAATCAACCGGTTGATTATGAAGCTGGACCTTGCGCCGTTCTATCTTATCCATGTCGATAATCTCGGTAAAggaaggtgcgaataagcggggaaattcttctcggctgac                                                                                                     <  2:376252‑M1/143‑40 (MQ=255)
                                      ccAGATCGGCAAGGAAGCTGGTGAAATCAACCGGTTGATTATCAAGCTGGACCTTGCGCCGTTCCATCTTATCCATGTCGATAATATCGTCAAAggaaggtgcgaataagcggggaaattcttctcggctgactcagccattg                                                                                           >  2:237914‑M1/1‑94 (MQ=255)
                                                  tgcaGCTGGTCGACTCAACCGGTTAATGATCAAGCGGGACCTTGCGCCGTTCCATCTTATCCATGTCGATAATATCGTTAAAGATATTCCCCAGCGTGACGGCCGAAACATGGATGGTCTTGAGATATTTTTCCTGCTcggcg                                                                               <  2:223027/140‑1 (MQ=255)
                                                                   aCCGGTGGATTATCAAGCTGGACCTTGCGCCGTTCCATCTTATCCATGTCGATAATATCGTTAAAggaaggtgcgaataagcggggaaattcttctcggctgactccgtcatttcatttcttcatgtttgagccgattttttc                                                              <  1:256378‑M1/143‑79 (MQ=255)
                                                                             tATCAAGCTGGACCTTGCGCCGTTCCATCTTATCCATGTCGATAATATCGTTAAAGATATTCCCCAGCGTGACGGCCGAAACATGGATGGTCTTGAGATATTTTTCCTGCTCGGCGGTGAGTTCGGTATCCAGCAGAATGCgg                                                    <  1:1040/143‑1 (MQ=255)
                                                                               tCAAGCTGGACCTTGCGCCGTTCCATCTTATCCATGTCGATAATATCGTTAACggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgactttttctcccgtaaatgc                                                  >  2:221279‑M1/1‑53 (MQ=255)
                                                                                                    tCCATCTGCTCCATGTCGATAATATCGTTAAAGATATTCCCCAGCGTGACGGCCGAAACATGGATGGTCTTGAGATATTTTTCCTGCTCGGCGGTGAGTTCGGTATCCAGCAGAATGCGGCTCAGACCGACGATACCGTTCAg                             <  2:358587/143‑1 (MQ=255)
                                                                                                             tCCATGTCGATAATATCGTGAAAggaaggtgcgaagaagcggggaaattcttcgcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcagcctatttagaccgtttcttc                    <  1:58751‑M1/143‑121 (MQ=255)
                                                                                                                           atCGTTAAAggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtcaatgccttgaatcagcctatttagaccgtttcttcgccatttaaggcgt      >  1:221940‑M1/1‑9 (MQ=255)
                                                                                                                               ttAAAggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcagcctatttagaccgtttcttcgccatttaaggcgttatc  >  1:338524‑M1/1‑5 (MQ=255)
                                                                                                                                   |                                                                                                                                          
CAATCCTTTTTGTTGCGCCTGCAATGCGGAGAGATTTTCCAGATCGGCAAGGAAGCTGGTGAAATCAACCGGTTGATTATCAAGCTGGACCTTGCGCCGTTCCATCTTATCCATGTCGATAATATCGTTAAAGATATTCCCCAGCGTGACGGCCGAAACATGGATGGTCTTGAGATATTTTTCCTGCTCGGCGGTGAGTTCGGTATCCAGCAGAATGCGGCTCAGACCGACGATACCGTTCAGCGGTGTACGCAATTCGTGACTGATGGT  >  NC_000913/3351892‑3352161

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