New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4153821 =27 (1.040)5 (0.220) 4/114 1.8 20.3% coding (1027/1152 nt) argE acetylornithine deacetylase
?NC_000913 4154127 = 16 (0.720)coding (721/1152 nt) argE acetylornithine deacetylase

ATCCATATGCAACTCACAGCAAGCGCAAATACGGTTAGAAGCGTCGCCACCGTGAATATGCCCGAGGTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAAT  >  NC_000913/4154062‑4154194
                                                                 |                                                                   
atCCTTATGCAACTCACAGCAAGCGCTATTACGGTTAGAAGCGTCGCCACCTTGAATATGCCCGAGGt                                                                   <  2:529597/68‑1 (MQ=255)
atCCATATGCAACTCACAGCAAGCGCAAATACGGTTAGAAGCGTCGCCACCGTGAATATGCCCGAGGt                                                                   >  2:868390/1‑68 (MQ=255)
    atatGCAAATCACAGAAAGCGCAAATACGGTTAGAAGAGTCGCCACCGTGAATATGACAGAGGTTGAg                                                               >  1:962132/1‑68 (MQ=255)
     tatGCATCTCACAGCTAGCGCAAATACGGCTAGAAGCGTCGCCACCGTGAATATGCCCGAGGTTGAGc                                                              <  2:12177/68‑1 (MQ=255)
     tatGCAACTCACAGCAAGCGCAAATACGGTTAGAAGCGTCGCCACCGTGAATATGCCCGAGGTTGAGc                                                              >  1:911336/1‑68 (MQ=255)
      atGAAACTAACAGCAAGCGCAAAAACGGTTAGAAGCGTCGACACAGTTAATATAACAGAGATTGAGCg                                                             >  2:488949/1‑68 (MQ=255)
          aaCTCACAGCAAGCGCAAATACGGTTAGAAGCGTCGCCACCGTGAATATGCCCGAGGTTGAGCGTAgg                                                         <  2:750983/68‑1 (MQ=255)
          aaCTCACAGCAAGCGCAAATACGGCTAGAAGCGTCGCCACCGTGAATATGCCCGAGGTTGAGCGTAgg                                                         <  2:884899/68‑1 (MQ=255)
                 tggtgaactactgtaccgaagcgccgtttattcaaacgttatgcccgaGGTTGAGCGTAGGGTATGGc                                                  >  1:863637‑M2/49‑68 (MQ=255)
                     gaactactgtaccgtagcgccgtttattcaaacgttatgcccgaGGTTGTGCGTAGGGTATGGCACgg                                              <  2:452861‑M2/24‑1 (MQ=255)
                     gaactactgtaccgaagcgccgtttattcaaacgttatgcccgaGGTTGAGCGTAGGGTATGGCACgg                                              <  1:47331‑M2/24‑1 (MQ=255)
                     aGCGCAAATACGGTTAGAAGCGTCGCCACCGTGAATATGCCCGAGGTTGAGCGTAGTGTATGGCACgg                                              <  2:91309/68‑1 (MQ=255)
                     aGCGCAAATACGGCTAGAAGCGTCGCCACCGTGAATATGCCCGAGGTTGAGCGTAGGGTATGGCACgt                                              <  1:212441/68‑2 (MQ=255)
                                  gaagcgccgtttattcaaacgttatgcccgaGGTTGAGCGTAGGGTATGGCACg                                               >  1:818628‑M2/32‑54 (MQ=255)
                                  gaagcgccgtttattcaaacgttatgcccgaGGTTGAGCGTAGGGTATGGCACg                                               <  2:818625‑M2/23‑1 (MQ=255)
                                               cACCGTTAATATGCCCGAGGTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCt                    <  2:840172/68‑1 (MQ=255)
                                               cACCGTGAATTTGCCCGAGGTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCt                    <  1:696988/68‑1 (MQ=255)
                                               cACCGTGAATATGCCCGAGGTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCt                    >  2:36256/1‑68 (MQ=255)
                                                 ccGTGAATATGCCCGAGGTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCttt                  <  2:593848/68‑1 (MQ=255)
                                                       atatGCCCGAGGTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGtt            >  1:616810/1‑68 (MQ=255)
                                                       atatGCCCGAGGTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGtt            >  1:552774/1‑68 (MQ=255)
                                                        tatgcccgaGGTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTa           >  2:527910‑M2/10‑68 (MQ=255)
                                                            cccgaGGTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATcgc       >  2:909762‑M2/6‑68 (MQ=255)
                                                                 ggTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAAt  >  1:501206/1‑68 (MQ=255)
                                                                 |                                                                   
ATCCATATGCAACTCACAGCAAGCGCAAATACGGTTAGAAGCGTCGCCACCGTGAATATGCCCGAGGTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAAT  >  NC_000913/4154062‑4154194

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