New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 516489815 (0.650)4 (0.180) 4/234 NT 16.7% coding (3030/3906 nt) hrpA ATP‑dependent helicase HrpA
?NC_002947 = 5164971 26 (1.200)coding (2957/3906 nt) hrpA ATP‑dependent helicase HrpA

TCACGCACAGTGCCGTTGTCTTCCACCAGTGCCGGGTACACGGTCATCGACAGGCCGGCGATCTTCTGCTGGGCGGTCTGTTTCACTTCGCTGAAGGCCTTGGCCTGCACCGGCTGCTCACCCTTGGCGTCACGCGGCACGGCCAGGGCAGCCTGGCTGGCGGCAGCAAAGCGCGCGGTCAGTTCCGCCAGGTCGCGGCCTTCGCCAAGGAACTTGCCCTGGCCATCGACCACTTCGATGTTCATGCGCAGGTGGCCCTCGACCTGGTTGACCGA  >  NC_002947/5164770‑5165044
                                                                                                                                |                                                                                                                                                  
tCACGCACAGTGCCGTTGTCTTCCACCAGTGCCGGGTACACGGTCATCGACAGGCCGGCGATCTTCTGCTGGGCGGTCTGTTTCACTTCGCTGAAGGCCTTGGCCTGCACCGGCTGCTCACCCTTGGCGaaggccgcgacctggcggaac                                                                                                                               >  1:18466‑M1/1‑129 (MQ=255)
                    ttCCACCAGTGCCGGGTACACGGTCATCGACAGGCCGGCGATCTTCTGCTGGGCGGTCTGTTTCACTTCGCTGAAGGCCTTGGCCTGCACCGGCTGCTCACCCTTGGCGTCACGCGGCACGGCCAGGGCAGCCTGGCTGGCGGCAGCaaa                                                                                                           >  1:7042/1‑150 (MQ=255)
                                           tCATCGACAGGCCGGCGATCTTCTGCTGGGCGGTCTGTTTCACTTCGCTGAAGGCCTTGGCCTGCACCGGCTGCTCACCCTTGGCGTCACGCGGCACGGCCAGGGCAGCCTGGCTGGCGGCAGCAAAGCGCGCGGTCAGTTCCGCCAGGt                                                                                    >  2:158690/1‑150 (MQ=255)
                                                        ggCGATCTTCTGCTGGGCGGTCTGTTTCACTTCGCTGAAGGCCTTGGCCTGCACCGGCTGCTCACCCTTGGCGTCACGCGGCACGGCCAGGGCAGCCTGGCTGGCGGCAGCAAAGCGCGCGGTCAGTTCCGCCAGGTCGCGGCCTTcgcc                                                                       <  2:466338‑M2/150‑5 (MQ=255)
                                                        ggCGATCTTCTGCTGGGCGGTCTGTTTCACTTCGCTGAAGGCCTTGGCCTGCACCGGCTGCTCACCCTTGGCGTCACGCGGCACGGCCAGGGCAGCCTGGCTGGCGGCAGCAAAGCGCGCGGTCAGTTCCGCCAGGTCGCGGCCTTcgcc                                                                       <  2:285372‑M2/150‑5 (MQ=255)
                                                            aTCTTCTGCTGGGCGGTCTGTTTCACTTCGCTGAAGGCCTTGGCCTGCACCGGCTGCTCACCCTTGGCGaaggccgcgacctggcggaactgaccgcgcgctttgctgccgccagccaggctgccctggccgtgccgcgtgacgccaagg                                                                   <  2:18466‑M1/150‑82 (MQ=255)
                                                                tctGCTGGGCGGTCTGTTTCACTTCGCTGAAGGCCTTGGCCTGCACCGGCTGCTCACCCTTGGCGTCGCGCGGCACGGCCAGGGCAGCCTGGCTGGCGGCAGCAAAGCGCGCGGTCAGTTCCGCCAGGTCGCGGCCTTCGCCAAGGAACt                                                               >  1:377085/1‑150 (MQ=255)
                                                                         cGGTTTGTTTCACTTCGCTGAAGGTCTTTGCCTGCACCGGCTGCTCACTCTTGGCGTCCCGCGGCACGGCCAGGGCAGCCTGGCTGGCGGCGGCAAAGCGCGCGGTCAGTTCCGCCAGGTCGCGGCCTTCGCCAAGGAACTTGCCCTGGc                                                      <  1:517394/150‑1 (MQ=255)
                                                                                tttCACTTCGCTGAAGGCCTTGGCCTGCACCGGCTGCTCACCCTTGGCGTCACGCGGCACGGCCAGGGCAGCCTGGCTGGCGGCAGCAAAGCGCGCGGTCAGTTCCGCCAGGTCGCGGCCTTCGCCAAGGAACTTGCCCTGGCCATCGAc                                               >  1:557057/1‑150 (MQ=255)
                                                                                  tCACTTCGCTGAAGGCCTTGGCCTGCACCGGCTGCTCACCCTTGGCGaaggccgcgacctggcggaactgaccgcgcgctttgctgccgccagccaggctgccctggccgtgccgcgtgacgccaaggg                                                                  >  2:566482‑M1/1‑47 (MQ=255)
                                                                                  tCACTTCGCTGAAGGCCTTGGCCTGCACCGGCTGCTCACCCTTGGCGaaggccgcgacctggcggaactgaccgcgcgctttgctgccgccagccaggctgccctggccgtgccgcgtgacgccaaggg                                                                  <  1:566482‑M1/129‑83 (MQ=255)
                                                                                      ttCGCTGAAGGCCTTGGCCTGCACCGGCTGCTCACCCTTGGCGTCACGCGGCACGGCCAGGGCAGCCTGGCTGGCGGCAGCAAAGCGCGCGGTCAGTTCCGCCAGGTCGCGGCCTTCGCCAAGGAACTTGccc                                                          >  1:33089/1‑133 (MQ=255)
                                                                                      ttCGCTGAAGGCCTTGGCCTGCACCGGCTGCTCACCCTTGGCGTCACGCGGCACGGCCAGGGCAGCCTGGCTGGCGGCAGCAAAGCGCGCGGTCAGTTCCGCCAGGTCGCGGCCTTCGCCAAGGAACTTGccc                                                          <  2:33089/133‑1 (MQ=255)
                                                                                                   ttGGCCTGCACCGGCTGCTCACCCTTGGCGTCACGCGGCACGGCCAGGGCAGCCTGGCTGGCGGCAGCAAAGGGCGCGGTCGGTTCCGCCAGGTCGCGGCCTTCGCCAAGGAACTTGCCCTGGCCATCGACCACTTCGATGTTCATgcgc                            >  1:230675/1‑150 (MQ=255)
                                                                                                    tGGCCTGCACCGGCTGCTCACCCTTGGCGTCACGCGGCACGGCCAGGGCAGCCTGGCTGGCGGCAGCAAAGCGCGCGGTCAGTTCCGCCAGGTCGCGGCCTTCGCCAAGGAACTTGCCCTGGCCATCGACCACTTCGATGTTCATGCGCa                           <  1:158690/150‑1 (MQ=255)
                                                                                                      gCCTGCACCGGCTGCTCACCCTTGGCGTCACGCGGCACGGCCAGGGCAGCCTGGCTGGCGGCAGCAAAGCGCGCGGTCAGTTCCGCCAGGTCGCGGCCTTCGCCAAGGAACTTGCCCTGGCCATCGACCACTTCGATGTTCATGCGCAgg                         >  2:75631/1‑150 (MQ=255)
                                                                                                        cTGCACCGGCTGCTCACCCTTGGCGTCACGCGGCACGGCCAGGGCAGCCTGGCTGGCGGCAGCAAAGCGCGCGGTCAGTTCCGCCAGGTCGCGGCCTTCGCCAAGGAACTTGCCCTGGCCATCGACCACTTCGATGTTCATGCGCAGGTg                       >  2:233076/1‑150 (MQ=255)
                                                                                                                    ctCACCCTTGGGGTCACGCGGCAGGGCCGGGGCAGCCTGGCTGGCGGCAGCAAAGCGCGCGGTCAGTTCCGCCAGGTCGCGGCCTTCGCCAAGGAACTTGCCCTGGCCATCGACCACTTCGATGTTCATGCGCAGGTGGCCCTCGACCTg           <  2:310255/150‑1 (MQ=255)
                                                                                                                     tCACCCTTGGCGTCACGCGGCACGGCCAGGGCAGCCTGGCTGGCGGCAGCAAAGCGCGCGGTCAGTTCCGCCAGGTCGCGGCCTTCGCCAAGGAACTTGCCCTGGCCATCGACCACTTCGATGTTCATGCGCAGGTGGCCCTCGACCTgg          <  1:613842/150‑1 (MQ=255)
                                                                                                                     tCACCCTTGGCGTCACGCGGCACGGCCAGGGCAGCCTGGCTGGCGGCAGCAAAGCGCGCGGTCAGTTCCGCCAGGTCGCGGCCTTCGCCAAGGAACTTGCCCTGGCCATCGACCACTTCGATGTTCATGCGCAGGTGGCCCTCGACCTgg          >  1:383610/1‑150 (MQ=255)
                                                                                                                        cccTTGGCGTCACGCGGCACGGCCAGGGCAGCCTGGCTGGCGGCAGCAAAGCGCGCGGTCAGTTCCGCCAGGTCGCGGCCTTcgccaagggtgagcagccggtgcaggccaaggccttcagcgaagtga                            <  1:566482‑M2/129‑48 (MQ=255)
                                                                                                                        cccTTGGCGTCACGCGGCACGGCCAGGGCAGCCTGGCTGGCGGCAGCAAAGCGCGCGGTCAGTTCCGCCAGGTCGCGGCCTTcgccaagggtgagcagccggtgcaggccaaggccttcagcgaagtga                            >  2:566482‑M2/1‑82 (MQ=255)
                                                                                                                         ccTTGGCGTCACGCGGCACGGCCAGGGCAGCCTGGCTGGCGGCAGCAAAGCGCGCGGTCAGTTCCGCCAGGTCGCGGCCTTcgccaagggtgagcagccggtgcaggccaaggccttcagcgaagtgaaacagaccgcccagcagaagat      <  2:18466‑M2/150‑70 (MQ=255)
                                                                                                                         ccTTGGCGTCACGCGGCACGGCCAGGGCAGCCTGGCTGGCGGCAGCAAAGCGCGCGGTCAGTTCCGCCAGGTCGCGGCCTTCGCCAAGGAACTTGCCCTGGCCATCGACCACTTCGATGTTCATGCGCAGGTGGCCCTCGACCTGGTTGa      >  1:126623/1‑150 (MQ=255)
                                                                                                                             ggCGaaggccgcgacctggcggaactgaccgcgcgctttgctgccgccagccaggctgccctggccgtgccgcgtgacgccaagggtgagcagccggtgcaggccaaggccttcagcgaagtgaaacagaccgcccagcagaagatcgcc  <  2:285372‑M1/150‑147 (MQ=255)
                                                                                                                             ggCGaaggccgcgacctggcggaactgaccgcgcgctttgctgccgccagccaggctgccctggccgtgccgcgtgacgccaagggtgagcagccggtgcaggccaaggccttcagcgaagtgaaacagaccgcccagcagaagatcgcc  <  2:466338‑M1/150‑147 (MQ=255)
                                                                                                                             ggCGTCACGCGGCACGGCCAGGGCAGCCTGGCTGGCGGCAGCAAAGCGCGCGGTCAGTTCCGCCAGGTCGCGGCCTTCGCCAAGGAACTTTCCCTGGCCATCGACCACTTCGATGTTCATGCGCAGGTGGCCCTCGACCTGGTTGACCGa  <  1:75631/150‑1 (MQ=255)
                                                                                                                                |                                                                                                                                                  
TCACGCACAGTGCCGTTGTCTTCCACCAGTGCCGGGTACACGGTCATCGACAGGCCGGCGATCTTCTGCTGGGCGGTCTGTTTCACTTCGCTGAAGGCCTTGGCCTGCACCGGCTGCTCACCCTTGGCGTCACGCGGCACGGCCAGGGCAGCCTGGCTGGCGGCAGCAAAGCGCGCGGTCAGTTCCGCCAGGTCGCGGCCTTCGCCAAGGAACTTGCCCTGGCCATCGACCACTTCGATGTTCATGCGCAGGTGGCCCTCGACCTGGTTGACCGA  >  NC_002947/5164770‑5165044

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.