New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_pgi = 453339911 (0.680)3 (0.200) 3/232 1.8 19.9% intergenic (‑86/+291) yjhX/yjhZ UPF0386 family protein/pseudogene, rimK paralog, C‑terminal fragment
?NC_000913_3_pae_pgi = 4533420 14 (0.940)intergenic (‑107/+270) yjhX/yjhZ UPF0386 family protein/pseudogene, rimK paralog, C‑terminal fragment
Rejected: Frequency below/above cutoff threshold.

AGCGAGAACGTGTAAGGTATGTTGTTCCTGACGGGATAAATTCATACTGAAGAACCTATTTAATCATCATGTGCAAAACGTGCAAACACACCGCGGTGTCCGCATTCGATTTCGGCGCATTGATAATCAGTCCGGCCTGAAAAGGTCGGGGAACTGATTATCGGATGATTACATTCTCCAGCATCAAAGCCTCGGGTTGAGTTGAAAGGTATTTACGGAGCGAATATTAACACGTAGTTATAGGAGGGATATGTTTCAGTAGAAAAAGTGTGATGT  >  NC_000913_3_pae_pgi/4533269‑4533544
                                                                                                                                  |                                                                                                                                                 
aGCGAGAACGTGTAAGGTATGTTGTTCCTGACGGGATAAATTCATACTGAAGAACCTATTTAATCATCATGTGCAAAACGTGCAAACACACCGCGGTGTCCGCATTCGATTTCGGCGCATTGATAATCAGTCCGGCCTGAAAAGGTCggg                                                                                                                                <  2:66278/150‑1 (MQ=255)
                 tATGTTGTTCCTGACGGGATAAATTCATACTGAAGAACCTATTTAATCATCATGTGCAAAACGTGCAAACACACCGCGGTGTCCGCATTCGATTTCGCCGCATTGATAATCAGTtccccgaccttttcaggccggactgcttatcaatgc                                                                                                               >  1:154192‑M1/1‑114 (MQ=255)
                                        ttCATACTGAAGAACCTATTTAATCATCATGTGCAAAACGTGCAAACACACCGCGGTGTCCGCATTCGATTTCGGCGCATTGATAATCAGTCCGGCCTGAAAAGGTCGGGGAACTGATTATCGGATGATTACATTCTCCAGCATGaa                                                                                           >  1:52441/1‑147 (MQ=255)
                                        ttCATACTGAAGAACCTATTTAATCATCATGTGCAAAACGTGCAAACACACCGCGGTGTCCGCATTCGATTTCGGCGCATTGATAATCAGTCCGGCCTGAAAAGGTCGGGGAACTGATTATCGGATGATTACATTCTCCAGCATGaa                                                                                           <  2:52441/147‑1 (MQ=255)
                                                                          aaaaCGTGCAAACACACCGCGGTGTCCGCATTCGATTTCGGCGCATTGATAATCAGTCCGGCCTGAAAAGGTCGGGGAACTGATTATCGGATGATTACATTCTCCAGCATCAAAGCCTCGGGTTGAGTTGAAAGGTATTTACGGAGCGaa                                                      >  1:102714/1‑150 (MQ=255)
                                                                                 gCAAACACACCGCGGTGTCCGCATTCGATTTCGGCGCATTGATAATCAGTCCGGCCTGAAAAGGTCGGGGAACTGATTATCGGATGATTACATTCTCCAGCATCAAAGCCTCGGGTTGAGTTGAAAGGTATTTACGGAGCGAATATTAac                                               >  1:162101/1‑150 (MQ=255)
                                                                                 gCAAACACACCGCGGTGTCCGCATTCGATTTCGGCGCATTGATAATCAGTCCGGCCTGAAAAGGTCGGGGAACTGATTATCGGATGATTACATTCTCCAGCATCAAAGCCTCGGGTTGAGTTGAAAGGTATTTACGGAGCGAATATTAac                                               <  2:102714/150‑1 (MQ=255)
                                                                                     acacacCGCGGTGTCCGCATTCGATTTCGGCGCATTGATAATCAGTtccccgaccttttcaggccggactgatt                                                                                                                       >  4:89282‑M1/1‑46 (MQ=255)
                                                                                     acacacCGCGGTGTCCGCATTCGATTTCGGCGCATTGATAATCAGTtccccgaccttttcaggccggactgatt                                                                                                                       <  3:89282‑M1/74‑29 (MQ=255)
                                                                                         acCGCGGTGTCCGCATTCGATTTCGGCGCATTGATAATCAGTCCGGCCTGAAAAGGTCGGGGAACTGATTATCGGATGATTACATTCTCCAGCATCAAAGCCTCGGGTTGAGTTGAAAGGTATTTACGGAGCGAATATTAACACGTAGtt                                       >  1:15120/1‑150 (MQ=255)
                                                                                           cgcgGTGTCCGCATTCGATTTCGGCGCATTGATAATCAGTCCGGCCTGAAAAGGTCGGGGAACTGATTATCGGATGATTACa                                                                                                         <  3:22430/82‑1 (MQ=255)
                                                                                           cgcgGTGTCCGCATTCGATTTCGGCGCATTGATAATCAGTCCGGCCTGAAAAGGTCGGGGAACTGATTATCGGATGATTACa                                                                                                         >  4:22430/1‑82 (MQ=255)
                                                                                           cgcgGTGTCCGCATTCGATTTCGGCGCATTGATAATCAGTCCGGCCTGAAAAGGTCGGGGAACTGATTATCGGATGATTACATTCTCCAGCATCAAAGCCTCGGGTTGAGTTGAAAGGTATTTACGGAGCGAATATTAACACGTAGTtat                                     <  2:126452/150‑1 (MQ=255)
                                                                                                                    gcATTGATAAGCAGTCCGGCCTGAAAAGGTCGGGGAactgattatcaatgcggcgaaatcgaatgcggacaccgcggtgtgtttgcacgttttgcacatgatgattaaataggttcttcagtatgaatttatcccgtcaggaacaacata            >  1:154192‑M2/1‑36 (MQ=255)
                                                                                                                            aaTCAGTCCGGCCTGAAAAGGTCGGGGAactgattatcaatgcgccgaaatcgaatgcggacaccgcggtgtgt                                                                                <  3:89282‑M2/74‑47 (MQ=255)
                                                                                                                            aaTCAGTCCGGCCTGAAAAGGTCGGGGAactgattatcaatgcgccgaaatcgaatgcggacaccgcggtgtgt                                                                                >  4:89282‑M2/1‑28 (MQ=255)
                                                                                                                              tCAGTCCGGCCTGAAAAGGTCGGGGAACTGATTATCGGATGATTACATTCTCCAGCATCAAAGCCTCGGGTTGAGTTGAAAGGTATTTACGGAGCGAATATTAACACGTAGTTATAGGAGGGATATGTTTCAGTAGAAAAAGTGTGATGt  >  2:75899/1‑150 (MQ=255)
                                                                                                                                  |                                                                                                                                                 
AGCGAGAACGTGTAAGGTATGTTGTTCCTGACGGGATAAATTCATACTGAAGAACCTATTTAATCATCATGTGCAAAACGTGCAAACACACCGCGGTGTCCGCATTCGATTTCGGCGCATTGATAATCAGTCCGGCCTGAAAAGGTCGGGGAACTGATTATCGGATGATTACATTCTCCAGCATCAAAGCCTCGGGTTGAGTTGAAAGGTATTTACGGAGCGAATATTAACACGTAGTTATAGGAGGGATATGTTTCAGTAGAAAAAGTGTGATGT  >  NC_000913_3_pae_pgi/4533269‑4533544

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