New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi 164559 =42 (1.060)6 (0.160) 5/268 4.8 14.1% intergenic (+25/‑171) hrpB/mrcB putative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
?NC_000913_3_bme_pgi = 3959482 34 (0.920)intergenic (+37/+50) ilvC/ppiC ketol‑acid reductoisomerase, NAD(P)‑binding/peptidyl‑prolyl cis‑trans isomerase C (rotamase C)
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AGGCGCAAAAGCCCGATTTTCTGTTCTTCTGTTACAAGACAGAAGAAGATATCTCTCAAAATTTGTGCTCGATCGGTCCCCTCTCCCTTGAGGGAGAGGGTTAGGGTGAGGGTGAACGTGCGCGCTGATGCCCTCACCCCGGCCCTCTCCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913_3_bme_pgi/164709‑164559
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccctctcccTGTGGGAGAGGGTCGGGGTGAGGGCATCAGCGCGCACTTAACCCGCAACAGCAATACGTTTCATATCTGTCATATAGCCGCGCAGTTTCTTACCTACCTGCTCAATCGCATGGCTGCGAATCGCTTCGTTCACATCACGCAGTTGCCC  <  NC_000913_3_bme_pgi/3959482‑3959335
                                                                                                                                                                                                                                                                                                           
AGGCGCAAAAGCCCGATTTTCTGTTCTTCTGTCACAAGACAGAAGAAGATATCTCTCAAAATTTGTGCTCGATCGGTCCCCTCCCCCTTGAGGGAGAGGGTTAGGGTGAGGGTGAACGTGCGCGCTGATGCCCTCACCCCGGCCCTCTC                                                                                                                                                        <  1:591602/149‑1
              GATTTTCTGTTCTTCTGTTACAAGACAGAAGAAGATATCTCTCAAAATTTGTGCTCGATCGGTCCCCTCTCCCTTGAGGGAGAGGGTTAGGGTGAGGGTGAACGTGCGCGCTGATGCCCTCACCCCGGCCCTCTCCCTGTGGGAGAGGG                                                                                                                                          <  2:224233/149‑1
                                                                                                                                       ACCCCGACCCTCTCCCACAGGGAGAGGGTCGGGGTGAGGGCATCAGCGCGCACTTAGCCCGCAACAGCAATACGTTTCATATCTGTCATATAGCAGCGCAGTTTCTTACCTACCTGCTCAAT                                            <  1:389064/122‑1
                                                                                                                                       ACCCCGACCCTCTCCCACAGGGAGAGGGTCGGGGTGAGGGCATCAGCGCGCACTTAGCCCGCAACAGCAATACGTTTCATATCTGTCATATAGCAGCGCAGTTTCTTACCTACCTGCTCAAT                                            >  2:389064/1‑122
                                                                                                                                         CCCGACCCTCTCCCACAGGGAGAGGGTCGGGGTGAGGGCATCAGCGCGCACTTAACCCGCAACAGCAATACGTTTCATATCTGTCATATAGCCGCGCAGTTTCTTACCTACCTGCTCAATCGCATGGCTGCGAATCGCTTCGTTCACAT               >  1:68053/1‑149
                                                                                                                                         CCCGACCCTCTCCCACAGGGAGAGGGTCGGGGTGAGGGCATCAGCGCGCACTTAACCCGCAACAGCAATACGTTT                                                                                         >  1:317097/1‑75
                                                                                                                                         CCCGACCCTCTCCCACAGGGAGAGGGTCGGGGTGAGGGCATCAGCGCGCACTTAACCCGCAACAGCAATACGTTT                                                                                         <  2:317097/75‑1
                                                                                                                                                      CTGTGGGAGAGGGTCGGGGTGAGGGCATCAGCGCGCACTTAACCCGCAACAGCAATACGTTTCATATCTGTCATATAGCCGCGCAGTTTCTTACCTACCTGCTCAATCGCATGGCTGCGAATCGCTTCGTTCACATCACGCAGTTGCCC  >  1:125327/1‑149
                                                                                                                                                                                                                                                                                                           
AGGCGCAAAAGCCCGATTTTCTGTTCTTCTGTTACAAGACAGAAGAAGATATCTCTCAAAATTTGTGCTCGATCGGTCCCCTCTCCCTTGAGGGAGAGGGTTAGGGTGAGGGTGAACGTGCGCGCTGATGCCCTCACCCCGGCCCTCTCCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913_3_bme_pgi/164709‑164559
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccctctcccTGTGGGAGAGGGTCGGGGTGAGGGCATCAGCGCGCACTTAACCCGCAACAGCAATACGTTTCATATCTGTCATATAGCCGCGCAGTTTCTTACCTACCTGCTCAATCGCATGGCTGCGAATCGCTTCGTTCACATCACGCAGTTGCCC  <  NC_000913_3_bme_pgi/3959482‑3959335

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