New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 1076410 =13 (0.780)3 (0.190) 3/274 1.8 21.9% intergenic (+4/‑217) mreB/mreC cell wall structural actin‑like protein MreB/cell shape‑determining protein MreC
?NC_002947 1076429 = 9 (0.560)intergenic (+23/‑198) mreB/mreC cell wall structural actin‑like protein MreB/cell shape‑determining protein MreC

TGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGAGCGCCCGGTTTACAGGTAGTACGTACAGTGCTACCTGTATGCGCTGGGCTGGCGCCCGTTTGGGGCGCCGTATCCGTCAACCCGCAACCAGGCTGGTGCGTTGTTCCATGTCCAATGACCTCCTGAGTCGTCCA  >  NC_002947/1076265‑1076559
                                                                                                                                                 |                                                                                                                                                     
tGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCcgctc                                                                                                                                                   <  1:178911/150‑1 (MQ=255)
  aCCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCcgctcgc                                                                                                                                                 <  1:357163/150‑1 (MQ=255)
   ccTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCcgctcgct                                                                                                                                                >  1:188933/1‑150 (MQ=255)
           aaCTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGa                                                                                                                                        <  1:232842/150‑1 (MQ=255)
               ggacgactcaggaggtcattggacatggaacaacgcaccagcctggttgcgggttgacggatacggcgccccaaacgggcgccagcccagcgcatacaggtagcactgtacgtactacctgtaaaccgggCGCTCGCTGTTCACGAgcgc                                                                                                                                    >  1:244751‑M1/131‑150 (MQ=255)
                              cGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGAGCGCCCGGTTTACAGgtag                                                                                                                     <  2:165789/150‑1 (MQ=255)
                                                 gCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGAGCGCCCGGTTTACAGGTAGTACGTACAGTGCTACCTGt                                                                                                  <  2:201946/150‑1 (MQ=255)
                                                                           ccGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGAGCGCCCGGTTTACAGGTAGTACGTACAGTGCTACCTGTATGCGCTGGGCTGGCGCCCGTTTggg                                                                        <  2:6660/150‑1 (MQ=255)
                                                                                       aaacgggcgccagcccagcgcatacaggtagcactgtacgtactacctgtaaaccgggCGCTCGCTGTTCACGAGCGCCCGGtt                                                                                                                              <  2:69510‑M1/26‑1 (MQ=255)
                                                                                       aaacgggcgccagcccagcgcatacaggtagcactgtacgtactacctgtaaaccgggCGCTCGCTGTTCACGAGCGCCCGGtt                                                                                                                              >  1:69510‑M1/59‑84 (MQ=255)
                                                                                                ggAGTTGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGAGCGGGCGGTTTACAGGTAGTACGTACAGTGGTAGCTGTATGCGCTGGGCTGGCGCCCGTTTGGGGCGCCGTATCCGTCAACCCGc                                                   >  2:145144/1‑150 (MQ=255)
                                                                                                ggAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGAGCGCCCGGTTTACAGGTAGTACGTACAGTGCTACCTGTATGCGCTGGGCTGGCGCCCGTTTGGGGCGCCGTATCCGTCAACCCGc                                                   >  2:186483/1‑150 (MQ=255)
                                                                                                                            cTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGAGCGCCCGGTTTACAGGTAGTACGTACAGTGCTa                                                                                                       <  2:71118/70‑1 (MQ=255)
                                                                                                                            cTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGAGCGCCCGGTTTACAGGTAGTACGTACAGTGCTa                                                                                                       >  1:71118/1‑70 (MQ=255)
                                                                                                                                             aTCCCGCTCGCTGTTCACGAGCGCCCGGTTTACAGGTAGTACGTACAGTGCTACCTGTATGCGCTGGGCTGGCGCCCGTTTGGGGCGCCGTATCCGTCAACCCGCAACCAGGCTGGTGCGTTGTTCCATGTc                        >  2:194616/1‑132 (MQ=255)
                                                                                                                                             aTCCCGCTCGCTGTTCACGAGCGCCCGGTTTACAGGTAGTACGTACAGTGCTACCTGTATGCGCTGGGCTGGCGCCCGTTTGGGGCGCCGTATCCGTCAACCCGCAACCAGGCTGGTGCGTTGTTCCATGTc                        <  1:194616/132‑1 (MQ=255)
                                                                                                                                                 CGCTCGCTGTTCACGAGCGCCCGGTTTACAGGTAGTACGTACAGTGCTACCTGTATGCGCTGGGCTGGCGCCCGTTTGGGGCGCCGTATCCGTCAACCCGCAACCAGGCTGGTGCGTTGTTCCATGTCCAATGACCTCCTGAGTCGTCCa  <  1:46671‑M1/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGAGCGCCCGGTTTACAGGTAGTACGTACAGTGCTACCTGTATGCGCTGGGCTGGCGCCCGTTTGGGGCGCCGTATCCGTCAACCCGCAACCAGGCTGGTGCGTTGTTCCATGTCCAATGACCTCCTGAGTCGTCCA  >  NC_002947/1076265‑1076559

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