New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 4599018 =47 (1.080)7 (0.170) 6/244 NT 12.8% coding (837/3627 nt) PP_4071 membrane protein
?NC_002947 4599066 = 50 (1.200)coding (789/3627 nt) PP_4071 membrane protein

GCCAAGCGTCGAAGGTGTCGACAGCATCCAGCTTGAACGTGAATCCAAACAGGCTATTGCGCTTGGCGGCCGAAAGCTTGCTGTAGAATTGGTCGAAGCCGTCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAGGGCCGTGCGTTGCTCGACGGTACCGTGCAGTAGCGCAGGCAGGTCGACAACCAGCACCAGGCCGTTGAGGGCTCTCTGGCTTCGATTGC  >  NC_002947/4598870‑4599157
                                                                                                                                                    |                                                                                                                                           
ggtaccgtctggaacccaagccaaactgtggcagcatctgttggattggctgcagcgcaatcgaagccagagagccctcaacggcctggtgctggttgtcgacctgcctgcgctactgcacggtaccgtcgagcaacgcacggcccttGc                                                                                                                                            >  1:702358‑M1/149‑150 (MQ=255)
  cAAGCGTCGAAGGTGTCGACAGCATCCAGCTTGAACGTGAATCCAAACAGGCTATTGCGCTTGGCGGCCGAAAGCTTGCTGTAGAATTGGTCGAAGCCGTCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGccc                                                                                                                                          <  2:171418/150‑1 (MQ=255)
        tCGCAGGTGTCGACAGCATCCAGCTTGAACGTGAATCCAAACAGGCTATTGCGCTTGGCGGCCGAAAGCTTGCTGTAGAATTGGTCGAAGCCGTCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTgg                                                                                                                                    <  2:128928/150‑1 (MQ=255)
         cGAAGGTGTCGACAGCATCCAGCTTGAACGTGAATCCAAACAGGCTATTGCGCTTGGCTGCCGAAAGCTTGCTGTAGAATTGGTCGAAGCCGTCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGgc                                                                                                                                   <  2:522298/150‑1 (MQ=255)
                    gcaaactgtggcagcatctgttggattggctgcagcgcaatcgaagccagagagccctcaacggcctggtgctggttgtcgacctgcctgcgctactgcacggtaccgtcgagcaacgcacggcccttGCCCAACTGGCTGCTCACCTCg                                                                                                                        <  1:880486‑M1/22‑1 (MQ=255)
                         aTCCAGCTTGAACGTGAATCCAAACAGGCTATTGCGCTTGGCGGCCGAAAGCTTGCTGTAGAATTGGTCGAAGCCGTCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAg                                                                                                                   >  1:832463/1‑150 (MQ=255)
                         aTCCAGCTTGAACGTGAATCCAAACAGGCTATTGCGCTTGGCGGCCGAAAGCTTGCTGTAGAATTGGTCGAAGCCGTCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAg                                                                                                                   >  2:1101673/1‑150 (MQ=255)
                          tCCAGCTTGAACGTGAATCCAAACAGGCTATTGCTCTTGGCGGCCGAAAGCTTGCTGTAGAATTGGTCGAAGCCGTCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGa                                                                                                                  <  1:973184/150‑1 (MQ=255)
                           ccAGCTTGAACGTGAATCCAAACAGGCTATTGCGCTTGGCGGCCGAAAGCTTGCTGTAGAATTGGTCGAAGCCGTCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGAc                                                                                                                 >  1:531471/1‑150 (MQ=255)
                             aGCTTGAACGTGAATCCAAACAGGCTATTGCGCTTGGCGGCCGAAAGCTTGCTGTAGAATTGGTCGAAGCCGTCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACgt                                                                                                               <  2:345339/150‑1 (MQ=255)
                              gcagcatctgttggattggctgcagcgcaatcgaagccagagagccctcaacggcctggtgctggttgtcgacctgcctgcgctactgcacggtaccgtcgagcaacgcacggcccttGCCCACCTGCTGCTCACCTCGTAAAGACgtgt                                                                                                             <  2:521789‑M1/32‑1 (MQ=255)
                              gCTTGAACGTGAATCCAAACAGGCTATTGCGCTTGGCGGCCGAAAGCTTGCTGTAGAATTGGTCGAAGCCGTCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGccc                                                                                                                                          >  2:259736/1‑122 (MQ=255)
                              gCTTGAACGTGAATCCAAACAGGCTATTGCGCTTGGCGGCCGAAAGCTTGCTGTAGAATTGGTCGAAGCCGTCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGccc                                                                                                                                          <  1:259736/122‑1 (MQ=255)
                                  gAACGTGAATCCAAACAGGCTATTGCGCTTGGCGGCCGAAAGCTTGCTGTAGAATTGGTCGAAGCCGTCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTCAAGACGTGTgcg                                                                                                          <  2:45434/150‑1 (MQ=255)
                                    tctgttggattggctgcagcgcaatcgaagccagagagccctcaacggcctggtgctggttgtcgacctgcctgcgctactgcacggtaccgtcgagcaacgcacggcccttGCCCACCTGCTGCTCACCTCGTAAAGACGTGTGCgcag                                                                                                       >  1:521789‑M1/113‑150 (MQ=255)
                                     cgggAATCCAAACAGGCTATGGCGCTTGGCGGCCGAAAGCGTGCGTTAGAATTGGTCGAAGCCGTACAGCAGATCTAACGTGCTGACCACGACATACAGCTGCAGGCGCGAGCCCAACTGGCTGCTCACATCGTAAAGGGGGGTGCgcag                                                                                                       <  2:493578/147‑1 (MQ=255)
                                         aaTCCAAACAGGCTATTGCGCTTGGCGGCCGAAAGCTTGCTGTAGAATTGGTCGAAGCCGTCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACGGGTGCGCGGCAgg                                                                                                   >  1:272111/1‑150 (MQ=255)
                                            ccAAACAGGCTATTGCGCTTGGCGGCCGAAAGCTTGCTGTAGAATTGGTCGAAGCCGTCCAGCAGATCGAAGTTGCATAATGCGACATTCGGCGGGAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTgg                                                                                                >  1:979778/1‑150 (MQ=25)
                                                  tgcagcgcaatcgaagccagagagccctcaacggcctggtgctggttgtcgacctgcctgcgctactgcacggtaccgtcgagcaacgcacggcccttGCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAgg                                                                                          >  2:310845‑M1/99‑150 (MQ=255)
                                                  aGGCTATTGCGCTTGGCGGCCGAAAGCTTGCTGTAGAATTGGTCGAAGCCGTCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAgg                                                                                          >  1:49950/1‑150 (MQ=255)
                                                    gcgaTTGCGCTTGGCGGCCGAAAGCTTGCTGTAGAATTGGTCGAAGCCGTCCAGCAGATCGAACTTGTTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAGGGc                                                                                        <  2:647/147‑1 (MQ=255)
                                                       aTTGCGCTTGGCGGCCGAAAGCTTGCTGTAGAATTGGTCGAAGCCGTCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAGGGCCGt                                                                                     >  1:687617/1‑150 (MQ=255)
                                                            tcgaagccagagagccctcaacggcctggtgctggttgtcgacctgcctgcgctactgcacggtaccgtcgagcaacgcacggcccttGCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAGGGCCGTGCGtt                                                                                >  2:880486‑M1/89‑150 (MQ=255)
                                                            gcTTGGCGGCCGAAAGCTTGCTGTAGAATTGGTCGAAGCCGTCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAGGGCCGTGCGtt                                                                                <  2:707209/150‑1 (MQ=255)
                                                              ttGGCGGCCGAAAGCTTGCTGTAGAATTGGTCGAAGCCGTCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAGGGCCGTGCGTTGc                                                                              >  1:7961/1‑150 (MQ=255)
                                                              gaagccagagagccctcaacggcctggtgctggttgtcgacctgcctgcgctactgcacggtaccgtcgagcaacgcacggcccttGCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAg                                                                                           <  2:537994‑M1/51‑1 (MQ=255)
                                                              gaagccagagagccctcaacggcctggtgctggttgtcgacctgcctgcgctactgcacggtaccgtcgagcaacgcacggcccttGCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAg                                                                                           >  1:537994‑M1/87‑137 (MQ=255)
                                                                    gcCGAAAGCTGGCTGTCGAATTGCTCCAAGCCGTCCCGCAGATCGAACTTACTAACCGCGACATACAGAGGCAAGCGCGAGCCCAAATGCCTGCTCACCTCGTAAAGACGCGGGCGCCGCGGCGGGGCAAGGGCCGTGCGTTGCTCGCCg                                                                        >  2:1121861/1‑150 (MQ=255)
                                                                     ccGAAAGCTTGCTGTAGAATTGGTCGAAGCCGTCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAGGGCCGTGCGTTGCTCGACgg                                                                       >  2:807923/1‑150 (MQ=255)
                                                                             ttGCTGTAGAATTGGTCGAAGCCGTCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAGGGCCGTGCGTTGCTCGACGGTACCGTGc                                                               >  2:189081/1‑150 (MQ=255)
                                                                                 tGTAGAATTGGTCGAAGCCGGCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTCAAGAGGTGTGCGCAGCAGGTGGGCAAGGGCCGTGCGTTGCTCGACGGTACCGTGCAGTa                                                           >  1:562328/1‑150 (MQ=255)
                                                                                   tAGAATTGGTCGAAGCCGTCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAGGGCCGTGCGTTGCTCGACGGTACCGTGCCGTAgc                                                         >  1:639554/1‑150 (MQ=255)
                                                                                   tAGAATTGGTCGAAGCCGTCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAGGGCCGTGCGTTGCTCGACGGTACCGTGCAGTAgc                                                         >  1:1016854/1‑150 (MQ=255)
                                                                                      aaTTGGTCGAAGCCGTCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAGGGCCGTGCGTTGCTCGACGGTACCGTGCAGTAGCgca                                                      <  2:666492/150‑1 (MQ=255)
                                                                                        ttGGTCGAAGCCGTCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACgtgt                                                                                                             >  2:188259/1‑93 (MQ=255)
                                                                                        ttGGTCGAAGCCGTCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACgtgt                                                                                                             <  1:188259/93‑1 (MQ=255)
                                                                                        gtGGTCGAAGCCGTCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAGGGCCGTGCGTTGCTCGACGGTACCGTGCAGTAGCgcagg                                                    <  2:429901/149‑1 (MQ=255)
                                                                                         tGGTCGAAGCCGTCCAGCAGATCGAACTTGCTCAACACGAAATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAGGGCCGTGCGTTGCTCGACGGTACCGTGCAGTAGCgcaggc                                                   <  1:61970/150‑1 (MQ=255)
                                                                                          ggTCGAAGCCGTCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGTAGGTGGGCAAGGGCCGTGCGTTGCTCGACGGTACCGTGCAGTAGCgcaggca                                                  <  2:49950/150‑1 (MQ=255)
                                                                                           gTCGAAGCCGTCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACgt                                                                                                               <  2:997905/88‑1 (MQ=255)
                                                                                           gTCGAAGCCGTCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACgt                                                                                                               >  1:997905/1‑88 (MQ=255)
                                                                                           gTCGAAGCCGTCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAGGGCCGTGCGTTGCTCGACGGTACCGTGCAGTa                                                           >  2:1119292/1‑140 (MQ=255)
                                                                                           gTCGAAGCCGTCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAGGGCCGTGCGTTGCTCGACGGTACCGTGCAGTa                                                           <  1:1119292/140‑1 (MQ=255)
                                                                                           gTCGAAGCCGTCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAGGGCCGTGCGTTGCTCGACGGTACCGTGCAGTAGCgcaggcag                                                 >  1:97793/1‑150 (MQ=255)
                                                                                             cGAAGCCGTCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAGGGCCGTGCGTTGCTCGACGGTACCGTGCAGTAGCGCAGGCAGGt                                               <  2:7961/150‑1 (MQ=255)
                                                                                                    gTCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAGGGCCGTGCGTTGCTCGACGGTACCGTGCAGTAGCGCAGGCAGGTCGACAAc                                        <  2:531471/150‑1 (MQ=255)
                                                                                                      ccAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAGGGCCGTGCGTTGCTCGACGGTAcc                                                                   >  1:107565/1‑121 (MQ=255)
                                                                                                      ccAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAGGGCCGTGCGTTGCTCGACGGTAcc                                                                   <  2:107565/121‑1 (MQ=255)
                                                                                                        agcagATCGAACTTGCTCAACACGCCATCCCGCGGCAAGCGCGAGCCCAACTGGTTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAGGGCCGTGCGTTGCCCGACGGTACCGTGCAGTAGCGCAGGCAGGTCGACAACCAGc                                    <  2:272111/150‑1 (MQ=255)
                                                                                                               cGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAGGGCCGTGCGTTGCTCGACGGTACCGTGCAGTAGCGCAGGCAGGTCGACAACCAGCACCAGGc                             <  2:97793/150‑1 (MQ=255)
                                                                                                                          aaCACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAGGGCCGTGCGTTGCTCGACGGTACCGTGCAGTAGCGCAGGCAGGTCGACAACCAGCACCAGGCCGTTGAGGGct                  >  1:655896/1‑150 (MQ=255)
                                                                                                                             acGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAGGGCCGTGCGTTGCTCGACGGTACCGTGCAGTAGCGCAGGCAGGTCGACAACCAGCACCAGGCCGTTGAGGGctctc               >  2:293093/1‑150 (MQ=255)
                                                                                                                                 cATACAGCGGCAAGCGCGATCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAGGGCCGTGCGTTGCTCGACGGTACCGTGCAGTAGCGCAGGCAGGTCGACAACCAGCACCAGGCCGTTGAGGGCTCTCTGGc           <  1:674879/150‑1 (MQ=255)
                                                                                                                                  aTACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAGGGCCGTGCGTTGCTCGACGGTACCGTGCAGTAGCGCAGGCAGGTCGACAACCAGCACCAGGCCGTTGAGGGCTCTCTGGCt          >  1:887924/1‑150 (MQ=255)
                                                                                                                                    aCAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAGGGCCGTGCGTTGCTCGACGGTACCGTGCAGTAGCGCAGGCAGGTCGACAACCAGCACCAGGCCGTTGAGGGCTCTCTGGCTTc        >  2:713291/1‑150 (MQ=255)
                                                                                                                                      ggCGGCAAGCGCGAGCCCAAGTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAGGGCCGTGCGTGGCTCGCCGGTACCGTGCAGTAGCGCAGGCAGGTCGACAACCAGCACCAGGCCGTTGAGGGCTCTCTGGCTTAGa      <  1:1121861/149‑1 (MQ=255)
                                                                                                                                        gggCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAGGGCCGTGCGTTGCTCGACGGTACCGTGCAGTAGCGAAGGCAGGTCGACAACCAGCACCAGGCCGTTGAGGGCTCTCTGGCTTCGAtt    <  2:562328/149‑1 (MQ=255)
                                                                                                                                          gCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAGGGCCGTGCGTTGCTCGACGGTACCGTGCAGTAGCGCAGGCAGGTCGACAACCAGCACCAGGCCGTTGAGGGCTCTCTGGCTTCGATTgc  <  2:832463/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                           
GCCAAGCGTCGAAGGTGTCGACAGCATCCAGCTTGAACGTGAATCCAAACAGGCTATTGCGCTTGGCGGCCGAAAGCTTGCTGTAGAATTGGTCGAAGCCGTCCAGCAGATCGAACTTGCTCAACACGACATACAGCGGCAAGCGCGAGCCCAACTGGCTGCTCACCTCGTAAAGACGTGTGCGCAGCAGGTGGGCAAGGGCCGTGCGTTGCTCGACGGTACCGTGCAGTAGCGCAGGCAGGTCGACAACCAGCACCAGGCCGTTGAGGGCTCTCTGGCTTCGATTGC  >  NC_002947/4598870‑4599157

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