New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 376746621 (0.650)4 (0.130) 4/234 3.6 15.4% pseudogene (246/945 nt) rhsJ Rhs‑family protein
?NC_000913 = 3767490 24 (0.780)pseudogene (270/945 nt) rhsJ Rhs‑family protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAGTGATGAGAACCCGCATCATCTGCAGCAGCTCATTCGCCTGCCGGGTCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGA  >  NC_000913/3767344‑3767611
                                                                                                                          |                                                                                                                                                 
tACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAGTGATGAGACCCCGCATCATCTg                                                                                                                           >  7:62175/1‑147 (MQ=37)
                   tCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAGTGATGAGAACCCGca                                                                                                                                  <  2:75855/121‑1 (MQ=35)
                   tCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAGTGATGAGAACCCGca                                                                                                                                  >  1:75855/1‑121 (MQ=35)
                                     aCCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAGTGATGAGAACCCGCATCATCTGCAGCAGCTCATTCGCCTGCCGGGTCAGCAGTATGATgag                                                                                    <  8:62175/149‑1 (MQ=255)
                                       cATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAGTGATGAGAACCCGc                                                                                                                                   >  8:126877/1‑100 (MQ=37)
                                       cATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAGTGATGAGAACCCGc                                                                                                                                   <  7:126877/100‑1 (MQ=37)
                                                 tGCCGCTGGCGCTCATCAGCACGGAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAGTGATGAGAACCCGCATCATCTGCAGCAGCTCATTCGCCTGCCGGGTcagca                                                                                              <  8:85868/127‑1 (MQ=255)
                                                 tGCCGCTGGCGCTCATCAGCACGGAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAGTGATGAGAACCCGCATCATCTGCAGCAGCTCATTCGCCTGCCGGGTcagca                                                                                              >  7:85868/1‑127 (MQ=255)
                                                 tGCCGCTGGCGCTCATCAGCACGGAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAGTGATGAGAACCCGCATCATCTGCAGCAGCTCATTCGCCTGCCGGGTcagca                                                                                              <  6:161413/127‑1 (MQ=255)
                                                 tGCCGCTGGCGCTCATCAGCACGGAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAGTGATGAGAACCCGCATCATCTGCAGCAGCTCATTCGCCTGCCGGGTcagca                                                                                              >  5:161413/1‑127 (MQ=255)
                                                   ccGCTGGCGCTCATCAGCACGGAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAGTGATGAGAACCCGCATCATCTGCAGCAGCTCATTCGCCTGCCGGGTCAGCAGTATGATGAGGAGTCCGGCCTGTa                                                                      >  8:65045/1‑149 (MQ=255)
                                                           gcTCATCAGCACGGAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAGTGATGAGAACCCGCATCATCTGCAGCAGCTCATTCGCCTGCCGGGTCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAAcc                                                              <  7:65045/149‑1 (MQ=255)
                                                                                        ggtgCGCAGAATATGATGAATGGGGCAACCTGCTGAGTGATGAGAACCCGCATCATCTGCAGCAGCTCATTCGCCTGCCGGGTCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCTgggg                                 >  8:17525/1‑149 (MQ=255)
                                                                                          tgCGCAGAATATGATGAATGCGCCAACCTGCTGAGAGATGAGAACCAGCATCATCTGCAGCAGCTCATTCCCCAGCCGGGTCAGCAGTATGATGAGGAGTCCGCCCTCTATTACAAACGCCACCGCTATTATGACCCGCAGCTAGGGCg                               >  8:19400/1‑149 (MQ=17)
                                                                                               aGAATATGATGAATGGGGCAACCTGCTGAGTGATGAGAACCCGCATCATCTGCAGCAGCTCATTCGCCTGCCGGGTCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCTGGGGCGGtata                          <  7:19400/149‑1 (MQ=255)
                                                                                                        tgaATGGGGCAACCTGCTGAGTGATGAGAACCCGCATCATCTGCAGCAGCTCATTCGCCTGCCGGGTCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGc                                                            <  4:135385/106‑1 (MQ=255)
                                                                                                        tgaATGGGGCAACCTGCTGAGTGATGAGAACCCGCATCATCTGCAGCAGCTCATTCGCCTGCCGGGTCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGc                                                            >  3:135385/1‑106 (MQ=255)
                                                                                                          aaTGGGGCAACCTGCTGAGTGATGAGAACCCGCATCATCTGCAGCAGCTCATTCGCCTGCCGGGTCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGAt               >  8:113360/1‑149 (MQ=255)
                                                                                                          aaTGGGGCAACCTGCTGAGTGATGAGAACCCGCATCATCTGCAGCAGCTCATTCGCCTGCCGGGTCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGAt               >  4:54198/1‑149 (MQ=255)
                                                                                                                     ctgctgAGTGATGAGAACCCGCATCATCTGCAGCAGCTCATTCGCCTGCCGGGTCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACt    <  7:113360/149‑1 (MQ=255)
                                                                                                                        ctgAGTGATGAGAACCCGCATCATCTGCAGCAGCTCATTCGCCTGCCGGGTCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGa  <  7:17525/148‑1 (MQ=255)
                                                                                                                          |                                                                                                                                                 
TACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAGTGATGAGAACCCGCATCATCTGCAGCAGCTCATTCGCCTGCCGGGTCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGA  >  NC_000913/3767344‑3767611

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.