New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 5739035 =16 (0.590)6 (0.240) 6/270 2.1 29.5% intergenic (‑56/+101) proY/hutH proline (histidine) APC transporter/histidine ammonia‑lyase
?NC_002947 5739114 = 14 (0.560)intergenic (‑135/+22) proY/hutH proline (histidine) APC transporter/histidine ammonia‑lyase

AGGCCGGTGCCGATGGCGGAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTA  >  NC_002947/5738894‑5739181
                                                                                                                                              |                                                                                                                                                   
aGGCCGGTGCCGATGGCGGAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTcgcg                                                                                                                                            <  2:599120/150‑1 (MQ=255)
  ggtgccgcactacgaccgtgaccgtttcttcgcgccggatatcgacaaggccgtggacctgttggccaaggggagtctgaccggcttgctgccggcgggtgtactgccaagcatgtaacgcccctggggccgctgggcggCCCTTCGCgg                                                                                                                                           >  1:287868‑M1/141‑150 (MQ=255)
  ggtgccgcactacgaccgtgaccgtttcttcgcgccggatatcgaaaaggccgtggacctgttggccaagggtagtctgaccggcttgctgccggcgggtgtactgccaagcctgtaacgcccctggggccgctgggcggCCCTTCGCgg                                                                                                                                           >  1:287858‑M1/141‑150 (MQ=255)
     ggTGCCGATGGCGGAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCAc                                                                                                                                       <  2:228094/150‑1 (MQ=255)
             tGGCGGAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTc                                                                                                                               <  2:94776/150‑1 (MQ=255)
             tGGCGGAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTc                                                                                                                               <  2:94775/150‑1 (MQ=255)
                 ggAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTAGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGACCCTTCGCGGGCACGCCCGCTCCcac                                                                                                                           >  1:401139/1‑150 (MQ=255)
                       gAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGGACGCCCGCGCCCCCCggggg                                                                                                                     >  1:405981/1‑149 (MQ=255)
                             cATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTGTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGgtgt                                                                                                               <  2:510212/150‑1 (MQ=255)
                                gAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGtttt                                                                                                            <  1:144846/150‑1 (MQ=255)
                                      gATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGAGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTGTTGCACCt                                                                                                      >  2:549296/1‑150 (MQ=255)
                                           tcgaaaaggccgtggacctgttggccaagggtagtctgaccggcttgctgccggcgggtgtactgccaagcctgtaacgcccctggggccgctgggcggCCCTTCGCGGGCGCGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACtt                                                                                                  >  1:605554‑M1/100‑150 (MQ=255)
                                                   gccgtggacctgttggccaagggtagtctgaccggcttgctgccggcgggtgtactgccaagcctgtaacgccccgggggctgctgcgcagCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTggg                                                                                           >  1:136057‑M1/92‑150 (MQ=255)
                                                   cTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGGATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACGTCTGTGGGAg                                                                                         >  2:258313/1‑150 (MQ=255)
                                                   cTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGa                                                                                          <  2:204632/150‑1 (MQ=255)
                                                             tgttggccaagggtagtctgaccggcttgctgccggcgggtgtactgccaagcctgtaacgcccctggggccgctgggcggCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCtgt                                                                                              <  1:95647‑M1/55‑1 (MQ=255)
                                                             tgttggccaagggtagtctgaccggcttgctgccggcgggtgtactgccaagcctgtaacgcccctggggccgctgggcggCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCtgt                                                                                              >  2:95647‑M1/82‑136 (MQ=255)
                                                                 gagaCCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCcgcg                                                                           >  2:360716/1‑150 (MQ=255)
                                                                                   cggcttgctgccggcgggtgtactgccaagcctgtaacgcccctggggccgctgggcggCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCgg                                                                                      <  2:409780‑M1/63‑1 (MQ=255)
                                                                                   cggcttgctgccggcgggtgtactgccaagcctgtaacgcccctggggccgctgggcggCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCgg                                                                                      >  1:409780‑M1/60‑122 (MQ=255)
                                                                                                          ttGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCa                                  <  2:401139/150‑1 (MQ=255)
                                                                                                                tCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACCCCCGCGCCCCACCGGGATGGTGTTTTGAAACTACGTCGGTGGGGGCGGGCGCGCCCGCGAATGGCCGCCCCGCGGGCCCAGGGGGGGGAACGGCGTGGCAGTacac                            >  1:422191/1‑150 (MQ=255)
                                                                                                                  gTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATTGTGTTTTTGACCTAGTTGTTTGGGTGCGGGCGCGCCCGCGAAGGGGCGCCCAGCGGGCCCAGGGGCGTTAAATGCTTTGCAGTAAAccc                          >  2:203007/1‑150 (MQ=37)
                                                                                                                                           aGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTa  <  2:541159/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
AGGCCGGTGCCGATGGCGGAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTA  >  NC_002947/5738894‑5739181

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