New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_pgi = 316787921 (1.210)3 (0.190) 3/254 2.1 13.5% coding (368/396 nt) mqsA antitoxin for MqsR toxin; transcriptional repressor
?NC_000913_3_pae_pgi = 3167883 19 (1.170)coding (364/396 nt) mqsA antitoxin for MqsR toxin; transcriptional repressor
Rejected: Frequency below/above cutoff threshold.

CCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913_3_pae_pgi/3167761‑3167879
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gttctggatAAGCATCCAGAACTATTGAATGAAATCCGTTAATTAATGGAAAAGCGGGGGCAACTCGCTTTTTTCTTTACCGGAGTGTGAGCCAAATCACTCGCGATAAGTTTGAGCGAACAATGCCACAGGCGTATTGTTGTTGGGCTTATAATT  <  NC_000913_3_pae_pgi/3167883‑3167737
                                                                                                                                                                                                                                                                          
CCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACGCTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATAAGCATCCAGAACTATTGAATGAAATCCGTT                                                                                                                      <  3:94315/150‑1
                                                                   AGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATAAGCATCCAGAACTATTGAATGAAATCCGTTAA                                                                                                                    >  3:58920/1‑85
                                                                   AGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATAAGCATCCAGAACTATTGAATGAAATCCGTTAA                                                                                                                    <  4:58920/85‑1
                                                                                                                    GATAAGCATCCAGAACTATTGAATGAAATCCGTTAATTAATGGAAAAGCGGGGGCAACTCGCTTTTTTCTTTACCGGAGTGTGAGCCAAATCACTCGCGATAAGTTTGAGCGAACAATGACACAGGCGTATTGTTGTTGGGCTTATAATT  >  2:35078/1‑150
                                                                                                                                                                                                                                                                          
CCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913_3_pae_pgi/3167761‑3167879
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gttctggatAAGCATCCAGAACTATTGAATGAAATCCGTTAATTAATGGAAAAGCGGGGGCAACTCGCTTTTTTCTTTACCGGAGTGTGAGCCAAATCACTCGCGATAAGTTTGAGCGAACAATGCCACAGGCGTATTGTTGTTGGGCTTATAATT  <  NC_000913_3_pae_pgi/3167883‑3167737

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.