New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 591424026 (0.820)4 (0.160) 4/192 NT 17.7% intergenic (+76/+77) spuI/argA glutamylpolyamine synthetase/amino acid N‑acetyltransferase
?NC_002947 = 5914280 17 (0.690)intergenic (+116/+37) spuI/argA glutamylpolyamine synthetase/amino acid N‑acetyltransferase

TCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCTTCGAATTGCGCTGGTAGTTGTACAGCGAAGCCCGCGCCGCCGGCAGCCGCTCCACGCCACTCGGCGCAAAGCCGCGCT  >  NC_002947/5914133‑5914417
                                                                                                                                                   |                                                                                                                                         
tCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgg                                                                                                                                         <  1:554647‑M2/150‑3 (MQ=255)
gcgAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgg                                                                                                                                         >  2:339717‑M2/2‑148 (MQ=255)
          cTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggcgt                                                                                                                               >  2:638109‑M2/1‑138 (MQ=255)
                                  gtagaCGCCGTCGCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCGCGCCCGCTCCCGCAGGTTTATTTTCAACGCTGAAGCCTGTAACTTACTTCGGGGAGCGGGCGTACCCGCGAAGAGGGCGGCCCAGGc                                                                                                       >  2:426598/5‑150 (MQ=255)
                                         ccgccCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAgg                                                                                                <  2:554647/150‑1 (MQ=255)
                                           gccCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGTGCACGCTCGCCCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCGCGCGAAGAGTGCGGCCCAGGCTTACAgggg                                                                                              >  2:88234/1‑150 (MQ=255)
                                                   aGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGTGCACGCCCGCTCCCGCAGGGTAGTTTCGACGGCTGAAGCTTGTGACTTATCTGTGTGCGCGGGCGTGCTTTTGGAGAGGGTGGGCTAGGCTTCCAAGGGCTTCTCGa                                                                                      >  1:684475/1‑150 (MQ=255)
                                                      gcggcgTTTTTGTGTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGaaaa                                                                                   <  1:72385/150‑1 (MQ=255)
                                                                  tttgAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAgg                                                                                                        >  1:314456/1‑117 (MQ=255)
                                                                  tttgAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAgg                                                                                                        <  2:314456/117‑1 (MQ=255)
                                                                          cGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGTGATTGGCCACACGCCAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCTTCGAATTGCGCTGGTAg                                                               >  1:444265/1‑150 (MQ=255)
                                                                            gtgtCGCCCTCTCCGCGGGCACGCCAGCTCTCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACATGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGTGCTTCTCGAAAATCTTCGAATTGCGCTGGTAGtt                                                             <  1:426598/150‑1 (MQ=255)
                                                                                 gCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCTTCGAATTGCGCTGGTAGTTGTACa                                                        >  2:19673/1‑150 (MQ=255)
                                                                                 gCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGCGCTTCTCGAACATCTTCTAATTGCGCTGGTAGTTGTACa                                                        <  1:19673/150‑1 (MQ=255)
                                                                                 gCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGATGTCGCGTTGGACTTACCTGGgggagcgggggtgcccgcgaagcggccaacaccggtctcaaacaaaaacgccgcccttgaggggcggcgtgttccattacacg                                                        <  2:204120‑M2/150‑84 (MQ=255)
                                                                                  cccTCTTCGCGAGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTGACCTGTGGTGACCGGCGGGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCTTCGAATTGCGCTGGTAGATCTACAg                                                       >  1:275543/1‑150 (MQ=255)
                                                                                     tctTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCTTCGAATTGCGCTGGTAGTTGTACAGCGa                                                    >  1:339677/1‑150 (MQ=255)
                                                                                                   cccGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggcgtgcccgcgaagagggcgacaccggtctcaaacaaaaacgccgcccttgaggggcggcgttttcgtttacacggtatgcaggtaccagttg                                      <  1:846849‑M2/150‑102 (MQ=255)
                                                                                                    ccGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggcgtgcccgcgaagagggcgacaccggtctcaaacaaaaacgccgcccttgaggggcggcgttttcgtttacacggtatgcaggtaccagttgt                                     <  2:163716‑M2/150‑103 (MQ=255)
                                                                                                     cGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggcgtgcccgcgaagagggcgacaccggtctcaaacaaaaacgccgcccttgaggggcggcgttttcgtttacacggtatgcaggtaccagttgta                                    >  2:846849‑M2/1‑47 (MQ=255)
                                                                                                                         gAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCTTCGAATTGCGCTGGTAGTTGTACAGCGAAGCCCGCGCCGCCGGCAGCCGCTCCACGCCACTCgg                >  2:191586/1‑150 (MQ=255)
                                                                                                                                gAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCTTCGAATTGCGCTGGTAGTTGTACAGCGAAGCCCGCGCCGCCGGCAGCCGCTCCACGCCACTCGGCGCAAAg         >  2:622811/1‑150 (MQ=255)
                                                                                                                                      tgtgACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCTTCGAATTGCGCTGGTAGTTGTACAGCGAAGCCCGCGCCGCCGGCAGCCGCTCCACGCCACTCGGCGCAAAGCCgcgc   >  1:49991/1‑150 (MQ=255)
                                                                                                                                       gtgACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCTTCGAATTGCGCTGGTAGTTGTACAGCGAAGCCCGCGCCGCCGGCAGCCGCTCCACGCCACTCGGCGCAAAGCCGCGCt  >  1:258149/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                         
TCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCTTCGAATTGCGCTGGTAGTTGTACAGCGAAGCCCGCGCCGCCGGCAGCCGCTCCACGCCACTCGGCGCAAAGCCGCGCT  >  NC_002947/5914133‑5914417

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