New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_pgi 730508 =18 (1.110)4 (0.260) 4/240 1.6 19.0% coding (926/4194 nt) rhsC Rhs protein with putative toxin domain; putative neighboring cell growth inhibitor
?NC_000913_3_pae_pgi 730544 = 17 (1.100)coding (962/4194 nt) rhsC Rhs protein with putative toxin domain; putative neighboring cell growth inhibitor
Rejected: Frequency below/above cutoff threshold.

ATACCCTGCCGGGTTACACCGAATATGGCCGGGACAACGGCATCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTT  >  NC_000913_3_pae_pgi/730406‑730663
                                                                                                                                          |                                                                                                                       
aTACCCTGCCGGGTTACACCGAATATGGCCGGGACAACGGCATCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGcggcgg                                                                                                              <  1:102513/150‑1 (MQ=35)
                 acCGAATATGGCCGGGACAACGGCATCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAAt                                                                                             <  1:163945/150‑1 (MQ=255)
                  ccGAATATGGCCGGGACAACGGCATCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATa                                                                                            <  1:175123/150‑1 (MQ=255)
                    gccaccatccggccccggtatttatcatcgtaagtaaagctgcgcacctgtgtattgctgcggtcatacaccgccgccagttCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGaa                                                                                                                        >  4:14751‑M1/83‑120 (MQ=255)
                    gccaccatccggccccggtatttatcatcgtaagtaaagctgcgcacctgtgtattgctgcggtcatacaccgccgccagttCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGaa                                                                                                                        <  3:14751‑M1/38‑1 (MQ=255)
                         catccggccccggtatttatcatcgtaagtaaagctgcgcacctgtgtattgctgcggtcatacaccgccgccagttCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGt                                                                                     >  1:119778‑M1/78‑150 (MQ=255)
                                   aaCGGCATCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCttt                                                                           <  1:60442/150‑1 (MQ=255)
                                             gtctgtctGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGTGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGAcc                                                                                                     >  4:24397/1‑114 (MQ=255)
                                             gtctgtctGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGTGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGAcc                                                                                                     <  3:24397/114‑1 (MQ=38)
                                             gtctgtctGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTa                                                                                                          <  2:67306/109‑1 (MQ=38)
                                             gtctgtctGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTa                                                                                                          >  1:67306/1‑109 (MQ=38)
                                              tctgtctgCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTttact                                                                        >  2:98261/1‑142 (MQ=255)
                                              tctgtctgCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTttact                                                                        <  1:98261/142‑1 (MQ=255)
                                                tgtctgCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATa                                                              >  2:32248/1‑150 (MQ=255)
                                                tgtctgCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATa                                                              <  3:42869/150‑1 (MQ=255)
                                                                 acctgtgtattgctgcggtcatacaccgccgccagttcgccgcggggcgtccagccatagcgcaccagcggcgAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGAtg                                             >  1:119778‑M2/74‑150 (MQ=255)
                                                                     aCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGAtggtgg                                         <  2:184526/150‑1 (MQ=255)
                                                                       gACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGgcg                                       >  1:24616/1‑150 (MQ=255)
                                                                        aCCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGgcgc                                      <  1:10481/150‑1 (MQ=255)
                                                                                                  gccgcgccgcTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATc            >  2:64617/1‑150 (MQ=255)
                                                                                                    ttcgccgcggggcgtccagccatagcgcaccagcggcgAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGgc                                        <  3:14751‑M2/82‑1 (MQ=255)
                                                                                                    ttcgccgcggggcgtccagccatagcgcaccagcggcgAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGgc                                        >  4:14751‑M2/39‑120 (MQ=255)
                                                                                                         cgcggggcgtccagccatagcgcaccagcggcgAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTAcc     <  2:119778‑M2/117‑1 (MQ=255)
                                                                                                            tGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGtt  <  1:64617/150‑1 (MQ=255)
                                                                                                                                          |                                                                                                                       
ATACCCTGCCGGGTTACACCGAATATGGCCGGGACAACGGCATCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTT  >  NC_000913_3_pae_pgi/730406‑730663

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