New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4035360NA (NA)25 (0.290) 10/166 NT 100% intergenic (+207/‑171) hemG/rrsA protoporphyrin oxidase, flavoprotein/16S ribosomal RNA of rrnA operon
?NC_000913 = 4035384 0 (0.000)intergenic (+231/‑147) hemG/rrsA protoporphyrin oxidase, flavoprotein/16S ribosomal RNA of rrnA operon

GTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCGAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAA  >  NC_000913/4035281‑4035459
                                                                               |                                                                                                   
gTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCGAAG                                                                                 <  2:1935845‑M2/100‑1 (MQ=255)
   ctcCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCGAAGCGG                                                                              <  1:698582‑M2/100‑1 (MQ=255)
    tcCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCGAAGCGGC                                                                             <  1:1236802‑M2/100‑1 (MQ=255)
    tcCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCGAAGCGGC                                                                             <  2:2210298‑M2/100‑1 (MQ=255)
    tcCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCGAAGCGGC                                                                             <  1:1950861‑M2/100‑1 (MQ=255)
    tcCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCGAAGCGGC                                                                             >  1:334534‑M2/1‑100 (MQ=255)
                              aaaaTAAATGCTTGACTCTGTAGCGGGAAGGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCGAAGCGGCgcggggtgtgcataatacgccttccc                                                   <  1:818135‑M2/100‑27 (MQ=255)
                                                         aaaGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCgcggggtgtgcataatacgctttcccgctacagagtcaagcatttatttttgc                        <  1:1838240‑M2/100‑54 (MQ=255)
                                                            gCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCGAAGCGGCacggggtgtgcataatacgccttcccgctacagagtcaagcatttatttttgcttt                     <  1:2023993‑M2/100‑57 (MQ=255)
                                                              gTATTATGCACACCCCGCGCCGCTGAGAAAAAGCGAAGCGGCgcggggtgtgcataatacgccttcccgctacagagtcaagcatttatttttgctttct                   <  1:1144960‑M2/100‑59 (MQ=255)
                                                                attatGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCgcggggtgtgcataatacgctttcccgctacagagtcaagcatttatttttgctttctct                 <  1:2533059‑M2/100‑61 (MQ=255)
                                                                       gatctccgcgccGCTGAGAAAAAGCAAAGCGGCgcggggtgtgcataatacgctttcccgctacagagtcaagcatttatttttgctttctctgccggag          >  1:514328‑M2/6‑33 (MQ=255)
                                                                          ncccgCGCCGCTGAGAAAAAGCGAAGCGGCgcggggtgtgcataatacgccttcccgctacagagtccatcatttatttttgctttctctgccggagct        <  1:2066432‑M2/98‑70 (MQ=255)
                                                                          ccccGCGCCGCTGAGAAAAAGCGAAGCGGCgcggggtgtgcataatacgccttcccgctacagagtcaagcatttatttttgctttctctgccggagttc       <  2:1824441‑M2/100‑71 (MQ=255)
                                                                          ccccGCGCCGCTGAGAAAAAGCGAAGCGGCgcggggtgtgcataatacgccttcccgctacagagtcaagcatttatttttgctttctctgccggagttc       <  1:1105844‑M2/100‑71 (MQ=255)
                                                                          ccccGCGCCGCTGAGAAAAAGCGAAGCGGCgcggggtgtgcataatacgccttcccgctacagagtcaagcatttatttttgctttctctgccggagttc       <  1:711416‑M2/100‑71 (MQ=255)
                                                                          ccccGCGCCGCTGAGAAAAAGCAAAGCGGCgcggggtgtgcataatacgctttcccgctacagagtcaagcatttatttttgctttctctgccggagttc       <  1:532079‑M2/100‑71 (MQ=255)
                                                                           cccGCGCCGCTGAGAAAAAGCGAAGCGGCgcggggtgtgcataatacgccttcccgctacagagtcaagcatttatttttgctttctctgccggagttct      <  1:25583‑M2/100‑72 (MQ=255)
                                                                            ccGCGCCGCTGAGAAAAAGCGAAGCGGCgcggggtgtgcataatacgccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctc     <  2:1832441‑M2/100‑73 (MQ=255)
                                                                            ccGCGCCGCTGAGAAAAAGCGAAGCGGCgcggggtgtgcataatacgccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctc     <  1:1502691‑M2/100‑73 (MQ=255)
                                                                            ccGCGCCGCTGAGAAAAAGCAAAGCGGCgcggggtgtgcataatacgctttcccgctacagagtcaagcatttatttttgctttctctgccggagagatc     <  2:514322‑M2/100‑73 (MQ=255)
                                                                             ngcgccGCTGAGAAAAAGCGAAGCGGCgcggggtgtgcataatacgccttcccgctacagagtcaagcatttatttttgctttttctgccggagttctca    <  2:233906‑M2/99‑74 (MQ=255)
                                                                             cgcgCCGCTGAGAAAAAGCGAAGCGGCgcggggtgtgcataatacgccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctca    <  1:80826‑M2/100‑74 (MQ=255)
                                                                             cgcgCCGCTGAGAAAAAGCGAAGCGGCgcggggtgtgcataatacgccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctca    <  1:2249868‑M2/100‑74 (MQ=255)
                                                                              gcgcCGCTGAGAAAAAGCGAAGCGGCgcggggtgtgcataatacgccttcccgctacagagtcaagcatttatttttgctttctgtgccggagttctcag   <  1:2516794‑M2/100‑75 (MQ=255)
                                                                              gcgcCGCTGAGAAAAAGCGAAGCGGCgcggggtgtgcataatacgccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcag   <  2:2055048‑M2/100‑75 (MQ=255)
                                                                               cgccgcTGAGAAAAAGCGAAGCGGCgcggggtgtgcataatacgccttcccgctacagaggcaagcatttatttttgctttctctgccggagttctcag   <  2:114733‑M2/99‑75 (MQ=255)
                                                                               cgccgcTGAGAAAAAGCGAAGCGGCgcggggtgtgcataatacgccttcccgctacagaggcaagcatttatgtttgctttctctgccggagttctcagg  <  2:1935164‑M2/100‑76 (MQ=255)
                                                                               cgccgcTGAGAAAAAGCGAAGCGGCgcggggtgtgcataatacgccttaccgctacagagtcaagcatttatttttgctttctctgccggagttatcagg  <  2:1280488‑M2/100‑76 (MQ=255)
                                                                               |                                                                                                   
GTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCGAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAA  >  NC_000913/4035281‑4035459

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