New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 16585427 (0.740)5 (0.150) 5/226 3.2 15.5% coding (1125/2535 nt) mrcB fused glycosyl transferase and transpeptidase
?NC_000913 = 165878 29 (0.850)coding (1149/2535 nt) mrcB fused glycosyl transferase and transpeptidase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TCGCCCGGTAGAAGAGCTAAGCCTCGACCAGCAGGCGCTGTTAGTCGGTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCCCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAA  >  NC_000913/165737‑166004
                                                                                                                                             |                                                                                                                              
tCGCCCGGTAGAAGAGCTAAGCCTCGACCAGCAGGCGCTGTTAGTCGGTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGatta                                                                                                                         <  3:38212/149‑1 (MQ=255)
 cGCCCGGTAGAAGAGCTAAGCCTCGACCAGCAGGCGCTGTTAGTCGGTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGattat                                                                                                                        <  8:121274/149‑1 (MQ=255)
      ggTAGAAGAGCTAAGCCTCGACCAGCAGGCGCTGTTAGTCGGTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATc                                                                                                                   <  7:144020/149‑1 (MQ=255)
                  aaGCCTCGACCAGCAGGCGCTGTTAGTCGGTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGaa                                                                                                              >  6:145085/1‑142 (MQ=255)
                  aaGCCTCGACCAGCAGGCGCTGTTAGTCGGTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGaa                                                                                                              <  5:145085/142‑1 (MQ=255)
                                  gcgcTGTTAGTCGGTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGCGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTAt                                                                                                        >  8:98114/1‑132 (MQ=255)
                                  gcgcTGTTAGTCGGTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTAt                                                                                                        >  8:95708/1‑132 (MQ=255)
                                  gcgcTGTTAGTCGGTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTAt                                                                                                        <  7:98114/132‑1 (MQ=255)
                                  gcgcTGTTAGTCGGTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTAt                                                                                                        <  7:95708/132‑1 (MQ=255)
                                            tCGGTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCCCGTCCGCTgggg                                                                             >  3:1006/1‑149 (MQ=255)
                                                         aGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCCCGTCCGCTGGGGGTTCAGCcgcgcg                                                                >  8:21905/1‑149 (MQ=255)
                                                           gggCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCCCGTCCGCTGGGGGTTCAGCCGCGCGGt                                                              >  1:33960/1‑149 (MQ=255)
                                                                   aTCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGa                                                                                                      >  1:17923/1‑101 (MQ=255)
                                                                   aTCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGa                                                                                                      <  2:17923/101‑1 (MQ=255)
                                                                           cccGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCCCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTctcctc                                              <  7:21905/149‑1 (MQ=255)
                                                                                 gCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCCCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTcagcca                                         <  2:33960/148‑1 (MQ=255)
                                                                                              cTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCCCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACt                           <  6:99935/149‑1 (MQ=255)
                                                                                                                aaTCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCCCGTCCGCTgg                                                                               <  3:167777/79‑1 (MQ=255)
                                                                                                                aaTCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCCCGTCCGCTgg                                                                               >  4:167777/1‑79 (MQ=255)
                                                                                                                 aTCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCCCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAg              >  3:45207/1‑143 (MQ=255)
                                                                                                                 aTCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCCCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAg              <  4:45207/143‑1 (MQ=255)
                                                                                                                 aTCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCAATGACATGTTGAGTGCCCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGTGATCTCTCCTCAGCCAGCCtt                                    >  3:23822/1‑120 (MQ=255)
                                                                                                                 aTCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCAATGACATGTTGAGTGCCCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGTGATCTCTCCTCAGCCAGCCtt                                    <  4:23822/120‑1 (MQ=255)
                                                                                                                       tgctgcGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCCCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCaaaa  <  4:1006/149‑1 (MQ=255)
                                                                                                                              tCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCCCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGAtctctc                                                 >  4:85778/1‑95 (MQ=255)
                                                                                                                              tCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCCCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGAtctctc                                                 <  3:85778/95‑1 (MQ=255)
                                                                                                                                        cagcaacagATTATTGATCAAGAACTCTATGACATGTTGAGTGCCCGt                                                                                      <  4:74889/48‑1 (MQ=255)
                                                                                                                                        cagcaacagATTATTGATCAAGAACTCTATGACATGTTGAGTGCCCGt                                                                                      >  3:74889/1‑48 (MQ=255)
                                                                                                                                           caacagATTATTGATCAAGAACTCTATGACATGTTGAGAGCCCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATg                                >  3:60325/1‑99 (MQ=255)
                                                                                                                                           caacagATTATTGATCAAGAACTCTATGACATGTTGAGAGCCCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATg                                <  4:60325/99‑1 (MQ=255)
                                                                                                                                             |                                                                                                                              
TCGCCCGGTAGAAGAGCTAAGCCTCGACCAGCAGGCGCTGTTAGTCGGTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCCCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAA  >  NC_000913/165737‑166004

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