New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 5911055 =23 (0.990)4 (0.180) 4/246 NT 15.0% coding (1159/1359 nt) spuB glutamylpolyamine synthetase
?NC_002947 5911078 = 23 (1.020)coding (1136/1359 nt) spuB glutamylpolyamine synthetase

AGTTTTCATGCTCGGCGCGCTTGGTAGCGACATACCCGGTGATGAATTTCTTGCCCAGGTACTGCACCAGCGCACGGCTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAGGCTGGCGGCGATGGCCAGGTACGGGTTGGCGTCGGCGCCCGGCAGGCGGTTCTCGACGCGCCGGCTTTGCGGGCCGGCATCCGGTAC  >  NC_002947/5910908‑5911202
                                                                                                                                                   |                                                                                                                                                   
aaaccgccccgtcggcgtgcgcgtaccggatgccggcccgcaaagccggcgcgtcgagaaccgcctgccgggcgccgacgccaacccgtacctggccatcgccgccagcctactgtgtggctacatcggcatggtcgaaggtatcgaGGc                                                                                                                                                   >  2:417514‑M1/148‑150 (MQ=255)
           tCGGCGCGCTTGGTAGCGACATACCCGGTGATGAATTTCTTGCCCAGGTACTGCACCAGCGCACGGCTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTgg                                                                                                                                         >  2:292629/1‑149 (MQ=255)
           tCGGCGCGCTTGGTAGCGACATACCCGGTGATGAATTTCTTGCCCAGGTACTGCACCAGCGCACGGCTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTgg                                                                                                                                         <  1:292629/149‑1 (MQ=255)
           tCGGCGCGCTTGGTAGCGACATACCCGGTGATGAATTTCTTGCCCAGGTACTGCACCAGCGCACGGCTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGgc                                                                                                                                        <  2:29683/150‑1 (MQ=255)
                     tGGTAGCGACATACCCGGTGATGAATTTCTTGCCCAGGTACTGCACCAGCGCACGGCTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCt                                                                                                                              <  1:90503/150‑1 (MQ=255)
                                    cGGTGATGAATTTCTTGCCCAGGTACTGCACCAGCGCACGGCTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCa                                                                                                               >  2:379606/1‑150 (MQ=255)
                                      gTGATGAATTTCTTGCCCAGGTACTGCACCAGCGCACGGCTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATg                                                                                                             >  2:13013/1‑150 (MQ=255)
                                              tttCTTGCCCAGGTACTGCACCAGCGCACGGCTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTa                                                                                                     >  2:350626/1‑150 (MQ=255)
                                                      ccAGGTACTGCACCAGCGCACGGCTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCcacaca                                                                                             >  1:101901/1‑150 (MQ=255)
                                                                cACCAGCGCACGGCTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGCCGCAGATTGCGCCGTTCGTAGCCACGGCCCTGCACTGGAGCGCTGCCCTCGATACCTGCGACCATGCCGATGAAGCGACAAAGTAGGCTggc                                                                                   >  2:107284/1‑150 (MQ=255)
                                                                 aCCAGCGCACGGCTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAGGCTggcg                                                                                  >  1:76550/1‑150 (MQ=255)
                                                                   cAGCGCACGGCTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAGGCTggc                                                                                   <  1:130661/147‑1 (MQ=255)
                                                                   cAGCGCACGGCTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAGGCTggc                                                                                   >  2:130661/1‑147 (MQ=255)
                                                                             acgccaacccgtacctggccatcgccgccagcctactgtgtggctacatcggcatggtcgaaggtatcgaGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTc                                                                                                                    <  1:570941‑M1/34‑1 (MQ=255)
                                                                             acgccaacccgtacctggccatcgccgccagcctactgtgtggctacatcggcatggtcgaaggtatcgaGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTc                                                                                                                    >  2:570941‑M1/71‑104 (MQ=255)
                                                                               gTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAGGCTGGCGGCGATGGCCAGGTa                                                                    <  1:379606/150‑1 (MQ=255)
                                                                                          gTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAGGCTGGCGGCGATGGCCAGGTACGGGTTGGCGt                                                         >  2:400610/1‑150 (MQ=255)
                                                                                                  tcgccgccagcctactgtgtggctacatcggcatggtcgaaggtatcgaGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGa                                                                                                                  >  1:145941‑M1/50‑85 (MQ=255)
                                                                                                  tcgccgccagcctactgtgtggctacatcggcatggtcgaaggtatcgaGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGa                                                                                                                  <  2:145941‑M1/36‑1 (MQ=255)
                                                                                                              gTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAGGCTGGCGGCGATGGCCAGGTACGGGTTGGCGTCGGCGCCCGGCAGGCGGTtc                                     >  1:546639/1‑150 (MQ=255)
                                                                                                               tCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAGGCTGGCGGCGATGGCCAGGTACGGGTTGGCGTCGGCGCCCGGCAGGCGGTtct                                    >  2:77086/1‑150 (MQ=255)
                                                                                                                 aGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAGGCTGGCGGCGATGGCCAGGTACGGGTTGGCGTCGGCGCCCGGCCGGCCGTtcgcc                                  >  2:284966/1‑147 (MQ=255)
                                                                                                                   gggCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAGGCTGCATGCGATGGCCAGGTACGGGTTGGCGTCGGCGCCCGGCAGGCGGTTCTCGAc                                <  1:350626/150‑1 (MQ=255)
                                                                                                                        ggCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAGGCTGGCGGCGATGGCCAGGTACGGGTTGGCGTCGGCGCCCGGCAGGCGGTTCTCGACGCGcc                           <  2:546639/150‑1 (MQ=255)
                                                                                                                              ggTTGCGGCGTTCGTAGCCAATGTCGTTCACTGGCTCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAGGCTGGCGGCGATGGCCAGGTACGGGTTGGCGTCGGCGCCCGGCAGGCGGTTCTCGACGCGCCGGCttt                     <  1:107284/150‑1 (MQ=37)
                                                                                                                                ttGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAg                                                                                                    <  1:431892/69‑1 (MQ=255)
                                                                                                                                ttGCGGCGTTCGTAGCAACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAg                                                                                                    >  2:431892/1‑69 (MQ=255)
                                                                                                                                             aGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAGGCTGGCGGCGATGGCCAGGTACGGGTTGGCGTCGGCGCCCGGCAGGCGGTTCTCGACGCGCCGGCTTTGCGGGCCGGCATCCg      <  1:13013/150‑1 (MQ=255)
                                                                                                                                              gCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAGGCTGGCGGCGATGGCCAGGTACGGGTTGGCGTCGGCGCCCGGCAGGCGGTTCTCGACGCGCCGGCTTTGCGGGCCGGCATCCgg     <  1:77086/150‑1 (MQ=255)
                                                                                                                                               ccACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAGGCTGGCGGCGATGGCCAGGTACGGGTTGGCGTCGGCGCCCGGCAGGCGGTTCTCGACGCGCCGGCTTTGCGGGCCGGCATCCGGt    >  2:268922/1‑150 (MQ=255)
                                                                                                                                                 aCGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAGGCTGGCGGCGATGGCCAGGTACGGGTTGGCGTCGGCGCCCGGCAGGCGGTTCTCGACGCGCCGGCTTTGCGGGCCGGCATCCGGTAc  >  1:297617/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGTTTTCATGCTCGGCGCGCTTGGTAGCGACATACCCGGTGATGAATTTCTTGCCCAGGTACTGCACCAGCGCACGGCTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAGGCTGGCGGCGATGGCCAGGTACGGGTTGGCGTCGGCGCCCGGCAGGCGGTTCTCGACGCGCCGGCTTTGCGGGCCGGCATCCGGTAC  >  NC_002947/5910908‑5911202

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.