New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 3387807 =10 (0.520)3 (0.170) 3/242 NT 23.6% coding (545/795 nt) PP_2989 short‑chain dehydrogenase/reductase family oxidoreductase
?NC_002947 3387849 = 10 (0.550)coding (587/795 nt) PP_2989 short‑chain dehydrogenase/reductase family oxidoreductase

GAGTTCGGCATGTCCATCTACGGCGCTACCAAAGCCTTCGTGCTGTTTCTGTCACAGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTCGATGCTGCACTTGTCGGTTTCGACCGTCGCGAGTTGGTAACCATCCCGCCACTTCATGTTGCCGAG  >  NC_002947/3387701‑3387961
                                                                                                                                                    |                                                                                                                
acatgaagtggcgggatggttaccaactcgcgacggtcgaaacggaccagtgcagcatcgaccagctcaccgacctccatcacctccgggagtgttttgacgtcaaTGACGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCAtt                                                                                                                 >  2:60943‑M1/107‑150 (MQ=255)
   ttCGGCATGTCCATCTACGGCGCTACCAAAGCCTTCGTGCTGTTTCTGTCACAGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAACGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGAc                                                                                                              >  1:61095/1‑150 (MQ=255)
                                         gCTGTTTCTGTCACAGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGAGATGGAGGTCGGTGAGCt                                                                        >  1:69200/1‑150 (MQ=255)
                                             tttCTGTCACAGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCGCGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTc                                                                    >  1:30260/1‑150 (MQ=255)
                                                                          gTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTCGATGCTGCACTTGTCGGTTTCGACCGTcg                                       >  1:162369/1‑150 (MQ=255)
                                                                             aCCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTCGATGCTGCACTTGTCGGTTTCGACCGTCGCGa                                    >  2:283166/1‑150 (MQ=255)
                                                                               ccAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAACGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTCGATGCTGCACTTGTCGGTTTCGACCGTCGCGAGt                                  <  2:61095/150‑1 (MQ=255)
                                                                                 aaGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTCGATGCTGCACTTGTCGGTTTCGACCGTCGCGAGTTg                                <  1:171958/150‑1 (MQ=255)
                                                                                              tGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTCGATGCTGCACTTGTCGGTTTCGACCGTCGCGAGTTGGTAACCATCCCGc                   >  1:16092/1‑150 (MQ=255)
                                                                                                   gCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTCGATGCTGCACTTGTCGGTTTCGACCGTCGCGAGTTGGTAACCATCCCGCCACtt              <  2:16092/150‑1 (MQ=255)
                                                                                                        aatgccggcacgcgcccagatctcggtgcgggtcgctgccgtcaTTGACGTCAAAACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTCGATGCTGCACTGGTCCGTTTCGACCGTCGCGAGTTGGTAACCATCCCGCCACTTCATGt         >  2:60943‑M2/45‑150 (MQ=255)
                                                                                                          tgccggcacgcgcccagatctcggtgcgggtcgctgccggcaTTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGa                                                                           <  1:349435‑M2/40‑1 (MQ=255)
                                                                                                          tgccggcacgcgcccagatctcggtgcgggtcgctgccggcaTTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGa                                                                           >  2:349435‑M2/43‑82 (MQ=255)
                                                                                                          TGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTCGATGCTGCACTTGTCGGTTTCGACCGTCGCGAGTTGGTAACCATCCCGCCACTTCATGTTg       <  1:417418‑M1/150‑1 (MQ=255)
                                                                                                            CCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTCGATGCTGCACTTGTCGGTTTCGACCGTCGCGa                                    <  1:430895‑M1/119‑1 (MQ=255)
                                                                                                            CCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTCGATGCTGCACTTGTCGGTTTCGACCGTCGCGa                                    >  2:430895‑M1/1‑119 (MQ=255)
                                                                                                               GCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTCGATGCTGCACTTGTCGGTTTCGACCGTCGCGAGTTGGTAACCATCCCGCCACTTCATGTTGCCGAg  <  1:428499‑M1/150‑1 (MQ=255)
                                                                                                                        cGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGt                                                                                 <  1:362713/62‑1 (MQ=255)
                                                                                                                        cGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGt                                                                                 >  2:362713/1‑62 (MQ=255)
                                                                                                                                                    |                                                                                                                
GAGTTCGGCATGTCCATCTACGGCGCTACCAAAGCCTTCGTGCTGTTTCTGTCACAGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTCGATGCTGCACTTGTCGGTTTCGACCGTCGCGAGTTGGTAACCATCCCGCCACTTCATGTTGCCGAG  >  NC_002947/3387701‑3387961

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