New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 3976787 =NA (NA)11 (0.920) 9/228 0.2 100% intergenic (+177/‑203) PP_3501/PP_3502 transposase/transposase
?NC_002947 4578248 = 0 (0.000)intergenic (+63/‑74) PP_4054/glgX hypothetical protein/glycogen debranching enzyme

AAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGACTGTTCGCGCAGTACAGACCGTAGGAGAAATCACCCCTCCTCCTTTCAACACCCAAGCGC  >  NC_002947/3976655‑3976894
                                                                                                                                    |                                                                                                           
aaaaGCAGCTAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTa                                                                                                          >  1:42998/1‑136 (MQ=255)
           ctccgcgtgGCCGATTCGATACTTAGGCGCGGGTTAACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTTGACATTGGCTAGTAGGCCCGGCTCTTTCGCTGGTAAACCCGCTTGTATGGTAGGACtt                                                                                               <  2:71363/130‑1 (MQ=38)
                   gggCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGggag                                                                                       >  1:324816/1‑136 (MQ=255)
                                        cGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGCGAGGAGGAGGGCCATTACTCGCtt                                                                  >  2:329177/1‑136 (MQ=255)
                                             tcaCCCGCGAATACGGTTGCCGCGACAATGGTCAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGAc                                                             <  1:325729/134‑1 (MQ=255)
                                                  cgcgAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGACTGTTc                                                        <  2:324816/136‑1 (MQ=255)
                                                  cgccAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGACTGTTc                                                        <  1:289042/136‑1 (MQ=255)
                                                                                  gCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGACTGTTCGCGCAGTACAGACCGTAGGAGAAATCACCcct                        <  2:119367/136‑1 (MQ=255)
                                                                                                        gCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGACTGTTCGCGCAGTACAGACCGTAGGAGAAATCACCCCTCCTCCTTTCAACACCCAAgcgc  >  1:255967/1‑136 (MQ=35)
                                                                                                                                    |                                                                                                           
AAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGACTGTTCGCGCAGTACAGACCGTAGGAGAAATCACCCCTCCTCCTTTCAACACCCAAGCGC  >  NC_002947/3976655‑3976894

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