New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 464759517 (0.530)3 (0.130) 3/212 3.6 20.4% intergenic (+48/+83) fusB/PP_4112 elongation factor G 2/sulfate transporter
?NC_002947 = 4647625 11 (0.470)intergenic (+78/+53) fusB/PP_4112 elongation factor G 2/sulfate transporter
Rejected: Coverage evenness skew score above cutoff.

AGGGCCGGGCCAGCTTCTCGATGGAGTTCACCCGCTATGCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCAAACATTGTGCGCCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGACCGCAACAGTTTTCAATGCTGCTCGCCCGCCCGCTTCAACAGCTTCTTGCACCGCTCCGACAAATGCACCACGCGC  >  NC_002947/4647453‑4647741
                                                                                                                                              |                                                                                                                                                  
aGGGCCGGGCCAGCTTCTCGATGGAGTTCACCCGCTATGCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGc                                                                                                                                               <  2:106748/148‑1 (MQ=255)
            gCTTCTCGATGGAGTTCACCCGCTATGCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCaaa                                                                                                                                   <  7:141942/148‑1 (MQ=255)
                    aTGGAGTTCACCCGCTATGCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCAAACATTGTgc                                                                                                                           <  3:19167/148‑1 (MQ=255)
                         gTTCACCCGCTATGCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCAAACATTGTGCGcccc                                                                                                                      >  1:124705/1‑148 (MQ=255)
                                 gCTATGCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTTAACCCGCTCCTAAAGGGTGGCCGCAGGCCCACACCTTGTGCGCCCCCCGGtt                                                                                                                >  5:72129/1‑145 (MQ=38)
                                      gCCGAGGCCCGGGCGAGGATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCAAACATTGTGCGCCCCCCTGTAGGAGCgg                                                                                                         <  6:155594/148‑1 (MQ=255)
                                      gCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCAAACATTGTGCGCCCCCCTGTAGGAGc                                                                                                           <  2:158140/146‑1 (MQ=255)
                                        cGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCAAACATTGTGCGCCCCCCTGTAGGAGCGGGt                                                                                                       <  2:124705/148‑1 (MQ=255)
                                         gAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCAAACATTGTGCGCCCCCCTGTAGGAGCGGGtt                                                                                                      <  2:76416/148‑1 (MQ=255)
                                            gCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCAAACATTGTGCGCCCCCCTGTAGGAGCGGGTTTAc                                                                                                   <  6:154230/148‑1 (MQ=255)
                                                                                         gtgTAAGCGGTAGTCCTTTGCCGGCCTCTTCGGGGGTAAACCCGCTCCTACAGGAGGGCCGCCGGCCCAAACTTTGTGCGCCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGACCGCAACAGTTTTCAATGCTGCTc                                                      <  6:72129/146‑1 (MQ=37)
                                                                                                              cGGCTTCTTCCCGTGAAAACCCGCGCCTACCGGAGGCCCGCAGGCCCAACATTTGTGCGCCCCCCTGTAGGAGCGGTTTAACCCGCGAAGAGGCCGGCACAGACCGCAACAGTTTTCAATGCTGCTCGCCCGCCCGCTTCAACAGctt                                 <  2:20906/148‑1 (MQ=37)
                                                                                                                  ctctTCGCGGTTAAACCCGCTCCTACAGGATGGCCGCAGGCCCAAACATTGTGCGCCCCCCTGTAGGAGGGGTTTTACGCGCGAAGAGGCCGGCACAGACCGCAACAGTTTTCAATGCTGCTCGCCCGCCCGCTTCAACAGCTTCTTg                             <  6:32927/148‑1 (MQ=255)
                                                                                                                   tctTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCACACATTGTGCGCCCCCCTGTAGGAGCGGGTTTTCCCGCGAAGAGGCCGGCACAGACCGGAACAGTTTTCCATGCTGCTCGCCCGCCCGCGTCAACAGCTTCTTGc                            >  5:84322/1‑148 (MQ=39)
                                                                                                                        gcgGGTAACCCCGCGCCTACAGTAGGCCCGCAGGCCCAAACATTGTCCGCCCCCCTGTCGGAGCGGTTTTACCCGCGAAGAGGCCGGCACAGACCGCAACAGTTTTCAATGCTGCTCGCCCGCCCGCTTCAACAGCTTCTTGCAc                          <  3:37720/145‑1 (MQ=38)
                                                                                                                            gTAAACCCGCTCCTACAGGATGGCCGCAGGCCCAAACATTGTGCGCCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGACCGCAACAGTTTTCAATGCTGCTCGCCCGCCCGCTTCACCAGCTTCTTGCACCGCTCCg                   >  5:38002/1‑148 (MQ=255)
                                                                                                                                              gATGGCCGCAGGCCCAAACATTGTGCGCCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGACCGCAACAGTTTTCAATGCTGCTCGCCCGCCCGCTTCAACAGCTTCTTGCACCGCTCCGACAAATGCACCACgcgc  >  4:45845/1‑147 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
AGGGCCGGGCCAGCTTCTCGATGGAGTTCACCCGCTATGCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCAAACATTGTGCGCCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGACCGCAACAGTTTTCAATGCTGCTCGCCCGCCCGCTTCAACAGCTTCTTGCACCGCTCCGACAAATGCACCACGCGC  >  NC_002947/4647453‑4647741

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.