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.

CCGCTGGCAAGCTTGTATTACTTTGGTCGCCCGGTAGAAGAGCTAAGCCTCGACCAGCAGGCGCTGTTAGTCGGTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCCCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAG  >  NC_000913/165711‑165992
                                                                                                                                               |                                                                                                                                          
ccGCTGGCAAGCTTGTATTACTTTGGTCGACCGGTAGAAGAGCTAAGCCTCGACCAGCAGGCGCTGTTAGTCGGTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCAAAAACTGGCGCTGGAGCGACGTAATCTGGtgct                                                                                                                                       >  5:99935/1‑149 (MQ=255)
                  tACTTTGGTCGCCCGGTAGAAGAGCTAAGCCTCGACCAGCAGGCGCTGTTAGTCGGTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTgcaacagc                                                                                                                      >  8:144020/1‑148 (MQ=255)
                          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)
                                               ccTCGACCAGCAGGCGCTGTTAGTCGGTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTgcaacagc                                                                                                                      <  2:17530/119‑1 (MQ=255)
                                               ccTCGACCAGCAGGCGCTGTTAGTCGGTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTgcaacagc                                                                                                                      >  1:17530/1‑119 (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                                                                   >  4:167777/1‑79 (MQ=255)
                                                                                                                                          aaTCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCCCGTCCGCTgg                                                                   <  3:167777/79‑1 (MQ=255)
                                                                                                                                           aTCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCCCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAg  <  4:45207/143‑1 (MQ=255)
                                                                                                                                           aTCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCCCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAg  >  3:45207/1‑143 (MQ=255)
                                                                                                                                           aTCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCAATGACATGTTGAGTGCCCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGTGATCTCTCCTCAGCCAGCCtt                        <  4:23822/120‑1 (MQ=255)
                                                                                                                                           aTCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCAATGACATGTTGAGTGCCCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGTGATCTCTCCTCAGCCAGCCtt                        >  3:23822/1‑120 (MQ=255)
                                                                                                                                               |                                                                                                                                          
CCGCTGGCAAGCTTGTATTACTTTGGTCGCCCGGTAGAAGAGCTAAGCCTCGACCAGCAGGCGCTGTTAGTCGGTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCCCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAG  >  NC_000913/165711‑165992

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.