New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 2048837 =9 (0.410)3 (0.160) 3/222 NT 26.3% intergenic (+33/+106) PP_1821/PP_1822 glutathione S‑transferase family protein/hypothetical protein
?NC_002947 2048928 = 9 (0.470)intergenic (+124/+15) PP_1821/PP_1822 glutathione S‑transferase family protein/hypothetical protein

CTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCAGGTCACCTTGCCGGCTTCTTCAATGCGCAGGCATGGCACTTTCACCCGGCCGCCACCTTCCTGCAGGGCCTGGC  >  NC_002947/2048813‑2049076
                                                                                                                   |                                                                                                                                                    
cTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCt                                                                                                                    >  2:367773/1‑150 (MQ=255)
                         cttgtgggagcgggcgagcccgcgaagcaggcaccgcggtggttggcacgggctgcgcccgtgttcgcgggcacgcccgctcccacaaggGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCtt                                                                                           >  1:407792‑M2/91‑150 (MQ=255)
                         CTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCtt                                                                                           >  2:160272‑M1/1‑150 (MQ=255)
                              GGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACt                                                                                      >  2:535718‑M1/1‑150 (MQ=255)
                                  gcgggcgagcccgcgaagcaggcaccgcggtggttggcacgggctgcgcccgtgttcgcgggcacgcccgctcccacaaggGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATa                                                                                  >  2:253928‑M2/82‑150 (MQ=255)
                                     GGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACAt                                                                               >  1:224053‑M1/1‑150 (MQ=255)
                                        GTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCa                                                                            <  2:336121‑M1/150‑1 (MQ=255)
                                         tGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCAg                                                                           <  1:367773/150‑1 (MQ=255)
                                            ccGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCgtcagt                                                                                                                                       >  2:439272/1‑87 (MQ=255)
                                            ccGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCgtcagt                                                                                                                                       <  1:439272/87‑1 (MQ=255)
                                                         accgcggtggttggcacgggctgcgcccgtgttcgcgggcacgcccgctcccacaaggGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCAGGTCACCTTGCCGGctt                                                           <  2:407792‑M2/92‑1 (MQ=255)
                                                                     cAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCAGGTCACCTTGCCGGCTTCTTCAATGCGCa                                               >  1:142949/1‑150 (MQ=255)
                                                                         caccGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCAGGTCACCTTGCCGGCTTCTTCAAt                                                    <  2:149660/141‑1 (MQ=255)
                                                                         caccGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATAATTGCCTTGGACTCATACATCCAGGTCACCTTGCCGGCTTCTTCAAt                                                    >  1:149660/1‑141 (MQ=255)
                                                                                     tGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATggcg                                                                                                                              >  2:135015/1‑55 (MQ=255)
                                                                                     tGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATggcg                                                                                                                              <  1:135015/55‑1 (MQ=255)
                                                                                                  cgcccgctcccacaaggGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCAGGTCACCTTGCCGGCTTCTTCAATGCTCAGGCATGGCACTTTCACCCGGCCGCCACCt                  <  1:336121‑M2/133‑1 (MQ=255)
                                                                                                            cacaaggGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCAg                                                                           <  2:231097‑M2/76‑1 (MQ=255)
                                                                                                            cacaaggGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCAg                                                                           >  1:231097‑M2/8‑83 (MQ=255)
                                                                                                                  gGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCAGGTCACCTTGCCGGCTTCTTCAATGCGCAGGCATGGCACTTTCACCCGGCCGCCACCTTCCTGCAGGGCCTGGc  <  2:142949‑M2/149‑1 (MQ=255)
                                                                                                                   |                                                                                                                                                    
CTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCAGGTCACCTTGCCGGCTTCTTCAATGCGCAGGCATGGCACTTTCACCCGGCCGCCACCTTCCTGCAGGGCCTGGC  >  NC_002947/2048813‑2049076

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