New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 396867210 (0.450)3 (0.160) 3/112 1.8 21.9% coding (757/1104 nt) wecA UDP‑GlcNAc:undecaprenylphosphate GlcNAc‑1‑phosphate transferase
?NC_000913 = 3968690 13 (0.680)coding (775/1104 nt) wecA UDP‑GlcNAc:undecaprenylphosphate GlcNAc‑1‑phosphate transferase

CCTGCTCGAAACGACCCAGGGCAAAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATGGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTC  >  NC_000913/3968605‑3968738
                                                                   |                                                                  
ccTGCTCGAAACGACCCAGGGCAAACCCCATCCCATCAGCCCGGATACCGCTTTGTGGATAATCGACA                                                                    >  1:491040‑M1/1‑68 (MQ=255)
     tCGAATCGACCCAGGGCATAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCg                                                               <  2:491041/68‑1 (MQ=255)
           cGACCCAGGGCAAAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATg                                                         <  2:403035/68‑1 (MQ=255)
               ccAGGGCAAAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATGGata                                                     >  2:246039/1‑68 (MQ=255)
                    gCAAAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATGGATATGGtg                                                <  2:487286‑M2/68‑3 (MQ=255)
                         aCCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCAccatatccatta                                                        <  2:835053‑M1/55‑13 (MQ=255)
                         aCCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCAccatatccatta                                                        >  1:835052‑M1/1‑43 (MQ=255)
                         aCCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATGGATATGGtggcgat                                           >  1:126048‑M2/1‑61 (MQ=255)
                         aCCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATGGATATGGtggcgat                                           <  2:576042‑M2/68‑8 (MQ=255)
                         aCCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATGGATATGGtggcgat                                           <  2:284832‑M2/68‑8 (MQ=255)
                         aCCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATGGATATGGtggcgat                                           <  1:422160‑M2/68‑8 (MQ=255)
                         aCCCATCCCATCAGCCCGGTTACCGCCTTGTGGATAATCGCCATTCCGCTAATGGATATGGtggcgat                                           >  2:113089‑M2/1‑61 (MQ=255)
                         aCCCATCCCATCAGCCCGGATACCGCTTTGTGGATAATCGCCATTCCGCTAATGGATAAGGtggcgat                                           <  2:144861‑M2/68‑8 (MQ=255)
                           ccatcccatcAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATGGATATGGtggcgatta                                         >  2:595036‑M2/1‑59 (MQ=255)
                                      gCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATGGATATGGtggcgattatccacaaagcg                              >  2:38338‑M2/1‑48 (MQ=255)
                                       cccGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATGGATATGGTGGCGATTATGTACCGTCGcc                             >  2:250820/1‑68 (MQ=255)
                                                cGCTTTGTGGATAATCGCCAccatatccattagcggaatggcgattatccacaaagcggtaaccgggc                    >  2:38338‑M1/1‑20 (MQ=255)
                                                 gCTTTGTGGATAATCGCCATTCCGCTAATGGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAgg                   >  1:656241/1‑68 (MQ=255)
                                                 gCTTTGTGGATAATCGCCATTCCGCTAATGGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAgg                   >  1:177863/1‑68 (MQ=255)
                                                  cTTTGTGGATAATCGCCATTCCGCTAATGGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGc                  >  2:465846/1‑68 (MQ=255)
                                                       tgGATAATCGCCATTCCGCTAATGGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAg             <  1:538926/68‑1 (MQ=255)
                                                           tAATCGCCAccatatccattagcggaatggcgattatccacaaagcggtaaccgggctgatgggatgg         >  2:595036‑M1/1‑9 (MQ=255)
                                                             aTCGCCAccttatccattagcggaatggcgattatccacaaagcggtatccgggctgatgggatgggt       <  2:144861‑M1/68‑62 (MQ=255)
                                                             aTCGCCAccatatccattagcggaatggcgattatccacaaggcggtaaccgggctgatgggatgggt       >  2:113089‑M1/1‑7 (MQ=255)
                                                             aTCGCCAccatatccattagcggaatggcgattatccacaaagcggtaaccgggctgatgggatgggt       <  2:284832‑M1/68‑62 (MQ=255)
                                                             aTCGCCAccatatccattagcggaatggcgattatccacaaagcggtaaccgggctgatgggatgggt       <  1:422160‑M1/68‑62 (MQ=255)
                                                             aTCGCCAccatatccattagcggaatggcgattatccacaaagcggtaaccgggctgatgggatgggt       <  2:576042‑M1/68‑62 (MQ=255)
                                                             aTCGCCAccatatccattagcggaatggcgattatccacaaagcggtaaccgggctgatgggatgggt       >  1:126048‑M1/1‑7 (MQ=255)
                                                                 ccATTCCGCTAATGGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATtctct   >  1:216874/1‑68 (MQ=255)
                                                                  caccatatccattagcggaatggcgattatccacaaagcggtaaccgggctgatgggatgggttttgc  <  2:487286‑M1/68‑67 (MQ=255)
                                                                  cATTCCGCTAATGGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATtctctc  <  2:280454/68‑1 (MQ=255)
                                                                   |                                                                  
CCTGCTCGAAACGACCCAGGGCAAAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATGGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTC  >  NC_000913/3968605‑3968738

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