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

GCAAAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATGGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCA  >  NC_000913/3968625‑3968757
                                                                 |                                                                   
gCAAAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATGGATATGGtg                                                                   <  2:487286‑M2/68‑3 (MQ=255)
     aCCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATGGATATGGtggcgat                                                              >  1:126048‑M2/1‑61 (MQ=255)
     aCCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATGGATATGGtggcgat                                                              <  2:284832‑M2/68‑8 (MQ=255)
     aCCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATGGATATGGtggcgat                                                              <  2:576042‑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)
                             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)
                                             ccATTCCGCTAATGGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATtctct                      >  1:216874/1‑68 (MQ=255)
                                              cATTCCGCTAATGGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATtctctc                     <  2:280454/68‑1 (MQ=255)
                                                      tAATGGATATGGtggcgattatccacaaagcggtaaccgggctgatgggatgggt                          <  2:835053‑M2/55‑44 (MQ=255)
                                                      tAATGGATATGGtggcgattatccacaaagcggtaaccgggctgatgggatgggt                          >  1:835052‑M2/1‑12 (MQ=255)
                                                         tGGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAg          <  2:422161/68‑1 (MQ=255)
                                                            atatGGTGGCGATTATGTACCGTCGCCTGCGTAAACGCATGAGCCCATTCTCTCCTGACCGTCAGCat       >  2:614734/1‑68 (MQ=255)
                                                                ggtggCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTc   >  2:389697/1‑68 (MQ=255)
                                                                ggtggCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTc   >  1:144101/1‑68 (MQ=255)
                                                                 gtggCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTcc  >  1:614400/1‑67 (MQ=255)
                                                                 gtggCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCa  >  2:324665/1‑68 (MQ=255)
                                                                 |                                                                   
GCAAAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATGGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCA  >  NC_000913/3968625‑3968757

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