New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 4615591 =NA (NA)5 (0.110) 5/250 NT 15.3% coding (1962/4593 nt) PP_4085 rhs family protein
?NC_002947 4615615 = 29 (0.590)coding (1986/4593 nt) PP_4085 rhs family protein

GGTCTATAAGCACATCGACCCTCAGGGCAACTGCTGGCTGGCGGAAGTGATCGCCGGACAGTTGATCAAGCGCATCGACCCGCAAGGAGGCGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGAACACATCAAATACCTACGACATTGGGCGCTACCCATCCAGGTCACCGATACAGCCGGACGGACTCGGCAATACGGCTACGACAGCCACGGCAATCTGCTCTGGCA  >  NC_002947/4615468‑4615731
                                                                                                                                                   |                                                                                                                    
tcctgctgccagagcagattgccgtggctgtcgtagccgtattgccgagtccgtccggctgtatcggtgacctggatgggtagcgcccaatgtcgtaggtatttgatgtgttcgctgcggccgaggggatcgcgggtttcgatcaggcg                                                                                                                     >  4:105922‑M2/148‑149 (MQ=255)
    tataAGCACATCGACCCTCAGGGCAACTGCTGGCTGGCGGAAGTGATCGCCGGACAGTTGATCAAGCGCATCGACCCGCAAGGAGGCGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCg                                                                                                                 >  4:126148/1‑149 (MQ=255)
     ctgccagagcagattgccgtggctgtcgtagccgtattgccgagtccgtccggctgtatcggtgacctggatgggtagcgcccaatgtcgtaggtatttgatgtgttcgctgcggccgaggggatcgcgggtttcgatcaggCGGcc                                                                                                                  <  5:238474‑M2/5‑1 (MQ=255)
           ccctCGACCCTCAGGTCAACTGCTGGCTGGCGTAAGCGATCGCCGGACAGTTTATCAAGCGCATCGACCCGCAAGGAGTCTAATGGCGCTACATCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGaa                                                                                                          <  8:202819/146‑1 (MQ=18)
                 aCCCTCAGGGCAACTGCTGGCTGGCGGAAGTGATCGCCGGACAGTTGATCAAGCGCATCGACCCGCAAGGAGGCGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGAACAGAtg                                                                                                    >  7:139650/1‑148 (MQ=255)
                 aCCCTCAGGGCAACTGCTGGCTGGCGGAAGTGATCGCCGGACAGTTGATCAAGCGCATCGACCCGCAAGGAGGCGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGAACACATc                                                                                                    <  3:126148/149‑1 (MQ=255)
                     tCAGGGCAACTGCTGGCTGGCGGAAGTGATCGCCGGACAGTTGATCAAGCGCATCGACCCGCAAGGAGGCGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGaa                                                                                                          <  2:10254/139‑1 (MQ=255)
                     tCAGGGCAACTGCTGGCTGGCGGAAGTGATCGCCGGACAGTTGATCAAGCGCATCGACCCGCAAGGAGGCGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGaa                                                                                                          >  1:10254/1‑139 (MQ=255)
                     tCAGGGCAACTGCTGGCTGGCGGAAGTGATCGCCGGACAGTTGATCAAGCGCATCGACCCGCAAGGAGGCGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGAACACATCaaa                                                                                                 >  2:98228/1‑148 (MQ=255)
                      cAGGGAAACTGCTGGCTGGCGGAAGTGATCGACGGACAGTTGATCAAGCGCATCGACCCGCAAGGAGGCGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAATCCCGCGATCCCCTCGGCCGCAGCGAACACATCAAATa                                                                                               <  6:224702/149‑1 (MQ=255)
                               tgctgGCTGGCGGAAGTGATCGCCGGACAGTTGATCAAGCGCATCGACCCGCAAGGAGGCGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGAACACATCAAATACCTACGACa                                                                                      >  5:36450/1‑149 (MQ=255)
                                   ggctggcGGAAGTGATCGCCGGACAGTTGATCAAGCGCATCGACCCGCAAGGAGGCGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGAACACATCAAATACCTACGACAttcg                                                                                  >  2:80427/1‑147 (MQ=255)
                                        gcGGAAGTGATCGCCGGACAGTTGATCAAGCGCATCGACCCGCAAGGAGGCGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGAACACATCAAATACCTACGACAt                                                                                     <  3:188715/141‑1 (MQ=255)
                                        gcGGAAGTGATCGCCGGACAGTTGATCAAGCGCATCGACCCGCAAGGAGGCGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGAACACATCAAATACCTACGACAt                                                                                     >  4:188715/1‑141 (MQ=255)
                                          ggAAGTGATCGCCGGACAGTTGATCAAGCGCATCGACCCGCAAGGAGGCGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGAACACATCAAATACCTACGACATTGGGCGCTAc                                                                           >  8:264975/1‑149 (MQ=255)
                                            aaGTGATCGCCGGACAGTTGATCAAGCGCATCGACCCGCAAGGAGGCGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGcc                                                                                                                  <  4:107380/108‑1 (MQ=255)
                                                  tCGCCGGACAGTTGATCAAGCGCATCGACCCGCAAGGAGGCGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGAACACATCAAATACCTACGACATTGGGCGCTACCCATCCAg                                                                   <  3:71586/149‑1 (MQ=255)
                                                   cGCCGGACAGTTGATCAAGCGCATCGACCCGCAAGGAGGCGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGAACACATCAAATACCTACGACATTGGGCGCTACCCATCCAgg                                                                  >  6:99953/1‑149 (MQ=255)
                                                     ccGGACAGTTGATCAAGCGCATCGACCCGCAAGGAGGCGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCTCGATCCCCTTGGCCGCAGCGAACACATCAAATACCTACGACATTGGGCGCTACCCATCCAGGTc                                                                <  1:80427/149‑1 (MQ=255)
                                                       ggACAGTTGATCAAGCGCATCGACCCGCAAGGAGGCGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGAACACATCAAATACCTACGACATTGGGCGCTACCCATCCAGGTCAc                                                              <  7:264975/149‑1 (MQ=255)
                                                       ggACAGTTGAGCAAGCGCATCGACCCGCAAGGAGGCGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGAACACATCAAATACCTACGACATTGGGCGCTACCCATCCAGGTCAc                                                              <  7:227167/149‑1 (MQ=255)
                                                           aGTTGATCAAGCGCATCGACCCGCAAGGAGGCGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGAACACATCAACTACCTACGACATTGGGCGCTACCCATCCAGGTCACCGa                                                           >  5:115453/1‑148 (MQ=255)
                                                           aGTTGATCAAGCGCATCGACCCGCAAGGAGGCGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCCAAACCCGCGATCCCCTCGGCCGCAGCGAACACATCAAATACCTACGACATTGGGCGCTACCCATCCAGGTCACCGAt                                                          <  1:98228/149‑1 (MQ=255)
                                                                aTCAAGCGCATCGACCCGCAAGGAGGCGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATACCCTAGGCCGCAGCGAACACATCAAATACCTACGACATAGGGCGCTACCCATCCAGGTACCCCATACAGc                                                     >  1:281028/1‑149 (MQ=14)
                                                                 tCAAGCGCATCGACCCGCAAGGAGGCGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGAACACATCAAATACCTACGACATTGGGCGCTACCCATCCAGGTCACCGATACAGcc                                                    >  1:112875/1‑149 (MQ=39)
                                                                     gcgcATCGACCCGCAAGGAGGCGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGAACACATCAAATACCTACGACATTGGGCGCTACCCATCCAGGTCACCGATACAGCcggac                                                >  5:24148/1‑149 (MQ=39)
                                                                     gcgcATCGACCCGCAAGGAGGCGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGAACACATCAAATACCTACGACATTGGGCGCTACCCATCCAGGTCACCGATACAGCcggac                                                >  3:196740/1‑149 (MQ=39)
                                                                        cATCGACCCGCAAGGAGGCGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGAACACATCAAATACCTACGACATTGGGCGCTACCCATCCAGGTCACCGATACAGCcggacgga                                             <  4:63394/149‑1 (MQ=39)
                                                                            gACCCGCAAGGAGGCGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGAACACATCAAATACCTACGACATTGGGCGCTACCCATCCAGGTCACCGATACAGCCGGACGGACTc                                          <  5:99953/148‑1 (MQ=39)
                                                                                 gCAAGGAGGCGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGAACACATCAAATACCTACGACATTGGGCGCTACCCATCCAGGt                                                                 >  4:198737/1‑120 (MQ=39)
                                                                                 gCAAGGAGGCGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGAACACATCAAATACCTACGACATTGGGCGCTACCCATCCAGGt                                                                 <  3:198737/120‑1 (MQ=39)
                                                                                     ggaggCGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGAACACATCAAATACCTACGACATTGGGCGCTACCCATCCAGGTCACCGATACAGCCGGACGGACTCGGCAATACgg                                <  5:285233/149‑1 (MQ=34)
                                                                                                          tgtgttcgctgcggccgaggggatcgcgggtttcgatcaggCGGCCGCAGCGAACACATCAAATACCTACGACATTGGGCGCTACCCATCCAGGTCACCGATACAGCCGGACGGACTCGGCAATACGGCTACGACAGCCAc                   >  8:64548‑M2/42‑141 (MQ=255)
                                                                                                          tgtgttcgctgcggccgaggggatcgcgggtttcgatcaggCGGCCGCAGCGAACACATCAAATACCTACGACATTGGGCGCTACCCATCCAGGTCACCGATACAGCCGGACGGACTCGGCAATACGGCTACGACAGCCAc                   <  7:64548‑M2/100‑1 (MQ=255)
                                                                                                                  ctgcggccgaggggatcgcgggtttcgatcaggCGGCCGCAGCGAACACATCAAATACCTACGACATTGGGCGCTACCCATCCAGGTCACCGATACAGCCGGACGGACTCGGCAATACGGCTACGACAGCCACGGCAATCTGCTCTgg    <  4:194757‑M2/115‑1 (MQ=255)
                                                                                                                  ctgcggccgaggggatcgcgggtttcgatcaggCGGCCGCAGCGAACACATCAAATACCTACGACATTGGGCGCTACCCATCCAGGTCACCGATACAGCCGGACGGACTCGGCAATACGGCTACGACAGCCACGGCAATCTGCTCTgg    >  3:194757‑M2/34‑148 (MQ=255)
                                                                                                                   tgcggccgaggggatcgcgggtttcgatcaggCGGCCGCAGCGAACACATCAAATACCTACGACATTGGGCGCTACCCATCCAGGTCACCGATACAGCCGGACGGACTCGGCAATACGGCTACGACAGCCACGGCAATCTGCTCTggcc  >  6:274580‑M2/33‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                    
GGTCTATAAGCACATCGACCCTCAGGGCAACTGCTGGCTGGCGGAAGTGATCGCCGGACAGTTGATCAAGCGCATCGACCCGCAAGGAGGCGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGAACACATCAAATACCTACGACATTGGGCGCTACCCATCCAGGTCACCGATACAGCCGGACGGACTCGGCAATACGGCTACGACAGCCACGGCAATCTGCTCTGGCA  >  NC_002947/4615468‑4615731

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