New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 522020 =33 (1.140)6 (0.220) 6/278 2.7 16.1% coding (707/1428 nt) PP_0435 M23/M37 family peptidase
?NC_002947 522075 = 31 (1.120)coding (652/1428 nt) PP_0435 M23/M37 family peptidase
Rejected: Frequency below/above cutoff threshold.

TAGGTTTTGCCGCGGTTGGTGAAGCGCGCAGAGAGGATGTTACCGGTGCCGACGACCTTGCCATCCATCACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAGCTCTTGATCACCCCGTGGGCATAGGCGGTGCGTATCACCGGCTTGCTGATTTCACGGTTGAAGGTGTAGCCCTT  >  NC_002947/521882‑522165
                                                                                                                                          |                                                                                                                                                 
tAGGTTTTGCCGCGGTTGGTGAAGCGCGCAGAGAGGATGTTACCGGTGCCGACGACCTTGCCATCCATCACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCAcgc                                                                                                                                          >  5:15030/1‑148 (MQ=255)
acagcaaagtcagcaatctcgaaaccattcgcctgcccaaaaccgccaagggctacacctttaaccgtgaaatcagcaagccggtgatacgcaccgcctatgcccacggggtgatcaagagctcgctgtctgcttccaCCCAGCAcgc                                                                                                                                          <  4:104234‑M1/10‑1 (MQ=255)
                  tcgaaaccattcgcctgaccaaaaccgccaagggctacaccttcaaccgtgaaatcagcaagccggtgatacgcaccgcctatgcccacggggtgatcaagagctcgctgtctgcttccgCCCAGCACGCGGGCCATGTCCATGGTCa                                                                                                                        <  7:22867‑M1/28‑1 (MQ=255)
                     aaGCGCGCAGAGAGGATGTTACCGGTGCCGACGACCTTGCCATCCATCACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGc                                                                                                                     >  1:150135/1‑148 (MQ=255)
                          cggagagagGATGGTCCCGGTGCCGACGACCTTTCCATCCATCAGCTTCTGTTCGTAGACCACGTCGAACTCGTCACCAGGTCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTgg                                                                                                                <  8:132582/145‑1 (MQ=255)
                               agagGCTGTTACCGGTGCCGACGACCTTGCCATCCATCACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAg                                                                                                           <  7:71772/148‑1 (MQ=255)
                                 tgaccaaaaccgccaagggctacaccttcaaccgtgaaatcagcaagccggtgatacgcaccgcctatgcccacggggtgatcaagagctcgctgtctgcttccgCCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGc                                                                                                         <  2:51925‑M1/43‑1 (MQ=255)
                                      gTTACCGGTGCCGACGACCTTGCCATCCATAACCTTCTGTTCGTAGACCACGTCGAACTCGACAACAGGGCGAATATCCTGTCCGAAGTCGATGTCCTAACCCAGCACCCGGGCCATGTCAAAGGTCATGCTCTGGTACAGTCCGgcg                                                                                                    >  6:4422/1‑148 (MQ=255)
                                          ccgccaagggctacaccttcaaccgtgaaatcagcaagccggtgatacgcaccgcctatgcccacggggtgatcaagagctcgctgtctgcttccgCCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCt                                                                                                >  1:51925‑M1/97‑148 (MQ=255)
                                            ggTGCCGACGACCTTGCCATCCATCACCTTCTGTTCGTAGACCACGTCGAACTCGGCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTgg                                                                                              <  7:53390/148‑1 (MQ=255)
                                            gccaagggctacaccttcaaccgtgaaatcagcaagccggtgatacgcaccgcctatgcccacggggtgatcaagagctcgctgtctgcttccgCCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTgg                                                                                              >  8:22867‑M1/95‑148 (MQ=255)
                                             gTGCCGACGACCTTGCCATCCATCACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTggg                                                                                             >  4:83387/1‑148 (MQ=255)
                                                agggctacacctttaaccgtgaaatcagcaagccggtgatacgcaccgcctatgcccacggggtgatcaagagctcgctgtctgcttccaCCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCgg                                                                                          >  3:104234‑M1/91‑148 (MQ=255)
                                                  tacgacCTTGCCATCCATCACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGAGTTCGTAACCCAGCACGCGGGCCATGTCCAGGGTCATGCTGTGTGACAGGCCGGCGCGCTGGTCTGaa                                                                                        >  4:121126/2‑148 (MQ=255)
                                                  gacgacCTTGCCATCCATCACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGaa                                                                                        >  5:132033/1‑148 (MQ=255)
                                                  gaccacCTTGCCCTCCATCCCCTTCTGTTCGTAGACCACGTCGAACTCGGCCCCAGGGCGAATATCCTGTGCGAAGTCGATGGCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGaa                                                                                        <  1:154938/148‑1 (MQ=255)
                                                         ttGCCATCCATCACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACa                                                                                 >  8:121241/1‑148 (MQ=255)
                                                           gCCATCCATCACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCtg                                                                                                                   >  4:86107/1‑112 (MQ=255)
                                                            ccatccatCACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCtgt                                                                                                                  <  3:86107/112‑1 (MQ=255)
                                                              atccatcACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAg                                                                            >  1:120241/1‑148 (MQ=255)
                                                              atccatcACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAg                                                                            >  5:136937/1‑148 (MQ=255)
                                                                  atcACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAGctct                                                                        <  8:6024/148‑1 (MQ=255)
                                                                    cACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAGCTCTTg                                                                      <  2:120241/148‑1 (MQ=255)
                                                                       cTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGCGCTCTTGATc                                                                   >  2:42111/1‑148 (MQ=255)
                                                                                     ccACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAGCTCTTGATCACCCCGTGGGCATa                                                     >  5:43163/1‑148 (MQ=255)
                                                                                         gTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAGCTCTTGATCACCCCGTGGGCATAGGCg                                                 >  5:9791/1‑148 (MQ=255)
                                                                                             aaCTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAGCTCTTGATCACCCCGTGGGCATAGGCGGTGc                                             <  6:4971/148‑1 (MQ=255)
                                                                                               cTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAGCTCTTGATCACCCCGTGGGCATAGGCGGTGCGt                                           <  6:50389/148‑1 (MQ=255)
                                                                                                    ccccAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAGCTCTTGATCACCCCGTGGGCATAGGCGGTGCGTATCAc                                      >  2:114748/1‑148 (MQ=255)
                                                                                                          ggCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAGCTCTTGATCACCCCGTGGGCATAGGCGGTGCGTATCACCGGCtt                                >  6:81603/1‑148 (MQ=255)
                                                                                                             gAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAGCTCTTGATCACCCCGTGGGCATAGGCGGTGCGTATCACCGGCTTGc                              >  1:19141/1‑147 (MQ=255)
                                                                                                               atatCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAGCTCTTGATCACCCCGTGGGCATAGGCGGTGCGTATCACCGGCTTGCTGa                           <  4:149789/148‑1 (MQ=255)
                                                                                                                     tgtgCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAGCTCTTGATCACCCCGTGGGCATAGGCGGTGCGTATCACCGGCTTGCTGATTTCAc                     >  3:154575/1‑148 (MQ=255)
                                                                                                                             gTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAGCTCTTGATCACCCCGTGGGCATAGGCGGTGCGTATCACCGGCTTGCTGATTTCACGGTTGAAg             >  3:2804/1‑148 (MQ=255)
                                                                                                                                gATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAGCTCTTGATCACCCCGTggg                                                         >  1:10191/1‑101 (MQ=255)
                                                                                                                                gATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAGCTCTTGATCACCCCGTGGGCATAGGCGGTGCGTATCACCGGCTTGCTGATTTCACGGTTGAAGgt           >  6:154621/1‑147 (MQ=255)
                                                                                                                                 aTGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAGCTCTTGATCACCCCGTGGGc                                                        <  2:10191/101‑1 (MQ=255)
                                                                                                                                      gTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAGCTCTTGATCACCCCGTGGGCATAGGCGGTGCGTATCACCGGCTTGCTGATTTCACGGTTGAAGGTGTAGccc    >  8:51389/1‑148 (MQ=255)
                                                                                                                                        aaCCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAGCTCTTGATCACCCCGTGGGCATAGGCGGTGCGTATCACCGGCTTGCTGATTTCACGGTTGAAGGTGTAGCCCtt  >  5:79203/1‑148 (MQ=255)
                                                                                                                                          |                                                                                                                                                 
TAGGTTTTGCCGCGGTTGGTGAAGCGCGCAGAGAGGATGTTACCGGTGCCGACGACCTTGCCATCCATCACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAGCTCTTGATCACCCCGTGGGCATAGGCGGTGCGTATCACCGGCTTGCTGATTTCACGGTTGAAGGTGTAGCCCTT  >  NC_002947/521882‑522165

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.