New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 36291051 (0.040)29 (1.210) 17/170 0.2 96.7% coding (905/1053 nt) yhiM DUF2776 domain‑containing protein
?NZ_CP009273 3645397 = NA (NA)pseudogene (747/747 nt) BW25113_RS18155 hypothetical protein

AAAAGTCTGGCTACTGGCACTGGTATGGAGACGTACCTGTTCGTTAGCCAACAGAATACCGATGATTCCCGTCTTCACCTGCCTGTTTTGCCTTTTCCTGGCATCGTTTCTTGCGGAAATGGCGCAGACCGACATGGGATATTTTATTCCTTCGCGAGTTCTGGTCGGTTTGGGAGCGG  >  NZ_CP009273/3629016‑3629194
                                                                                         |                                                                                         
aaaaGTCTGGCTACTGGCACTGGTATGGAGACGTACCTGTTCGTTAGCCAACAGAATACCGATGATTCCCGTCTTCACCTGCCTGTTTTG                                                                                           >  1:339247‑M1/1‑90 (MQ=255)
   aGTCTGGCTACTGGCACTGGTATGGAGACGTACCTGTTCGTTAGCCAACAGAATACCGATGATTCCCGTCTTCACCTGCCTGTTTTGgaa                                                                                        <  1:318537‑M1/90‑4 (MQ=255)
       tGGCTACTGGCACTGGTATGGAGACGTACCTGTTCGTTAGCCAACAGAATACCGATGATTCCCGTCTTCACCTGCCTGTTTTGgaaggtg                                                                                    >  2:310817‑M1/1‑83 (MQ=255)
                  aCTGGTATGGAGACGTACCTGTTCGTTAGCCAACAGAATACCGATGATTCCCGTCTTCACCTGCCTGTTTTGgaaggtgcgaataagcgg                                                                         >  2:369466‑M1/1‑72 (MQ=255)
                  aCTGGTATGGAGACGTACCTGTTCGTTAGCCAACAGAATACCGATGATTCCCGTCTTCACCTGCCTGTTTTGgaaggtgcgaataagcgg                                                                         >  1:589281‑M1/1‑72 (MQ=255)
                        aTGGAGACTTCCCTGTTGGTTAGCCAACAAAATCCCGATGATTCCCGTCTCCCCCTGCCTGTTTTGgaaggtgcgaataagcggggaaat                                                                   <  1:595436‑M1/90‑25 (MQ=255)
                               cGTACCTGTTCGTTAGCCAACAGAATACCGATGATTCCCGTCTTCACCTGCCTGTTTTGgaaggtgcgaataagcggggaattcttcttg                                                            >  1:616627‑M1/1‑59 (MQ=255)
                               cGTACCTGTTCGTTAGCCAACAGAATACCGATGATTCCCGTCTTCACCTGCCTGTTTTGgaaggtgcgaataagcggggaaattcttctc                                                            >  1:223560‑M1/1‑59 (MQ=255)
                               cGTACCTGTTCGTTAGCCAACAGAATACCGATGATTCCCGTCTTCACCTGCCTGTTTTGgaaggtgcgaataagcggggaaattcttctc                                                            >  2:89268‑M1/1‑59 (MQ=255)
                               cGTACCTGTTCGTTAGCCAACAGAATACCGATGATTCCCGTCTTCACCTGCCTGTTTTGgaaggtgcgaataagcggggaaattcttctc                                                            >  2:628910‑M1/1‑59 (MQ=255)
                               cGTACCTGTTCGTTAGCCAACAGAATACCGATGATTCCCGTCTTCACCTGCCTGTTTTGgaaggtgcgaataagcggggaaattcttctc                                                            >  2:577765‑M1/1‑59 (MQ=255)
                               cGTACCTGTTCGTTAGCCAACAGAATACCGATGATTCCCGTCTTCACCTGCCTGTTTTGgaaggtgcgaataagcggggaaattcttctc                                                            >  1:628201‑M1/1‑59 (MQ=255)
                               cGTACCTGTTCGTTAGCCAACAGAATACCGATGATTCCCGTCTTCACCTGCCTGTTTTGgaaggtgcgaataagcggg                                                                        <  1:401548‑M1/78‑20 (MQ=255)
                               cGTACCTGTTCGTTAGCCAACAGAATACCGATGATTCCCGTCTTCACCTGCCTGTTTTGgaaggtgcgaataagcggg                                                                        >  2:401548‑M1/1‑59 (MQ=255)
                               cGTACCTGTTCGTTAGCCAACAGAATACCGATGATTCCCGTCTTCACCTGCCTGTTTTGCCTTTTCCTGGCATCGTTTCTTGCGGAAATg                                                            >  2:228207/1‑90 (MQ=255)
                                gTACCTGTTCGTTAGCCAACAGAATACCGATGATTCCCGTCTTCACCTGCCTGTTTTGgaaggtgcgaataagcggggaaattcttctcg                                                           >  1:401208‑M1/1‑58 (MQ=255)
                                          gTTAGCCAACAGAATACCGATGATTCCCGTCTTCACCTGCCTGTTTTGgaaggtgcgaataagcggggaaattcttctcggctgactcag                                                 <  2:161295‑M1/90‑43 (MQ=255)
                                           ttAGCCAACAGAATACCGATGATTCCCGTCTTCACCTGCCTGTTTTGgaaggtgcgaataagcggggaaattcttctcggctgactcagt                                                >  2:435896‑M1/1‑47 (MQ=255)
                                            tAGCCAACAGAATACCGATGATTCCCGTCTTCACCTGCCTGTTTTGgaaggtgcgaataagcggggaaattcttctcggctgactcagtc                                               <  2:628201‑M1/90‑45 (MQ=255)
                                                       aTACCGATGATTCCCGTCTTCACCTGCCTGTTTTGgaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttct                                    >  1:493148‑M1/1‑35 (MQ=255)
                                                        tACCGATGATTCCCGTCTTCACCTGCCTGTTTTGgaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttctt                                   >  2:578017‑M1/1‑34 (MQ=255)
                                                            gatgatTCCCGTCTTCACCTGCCTGTTTTGgaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatg                               >  1:425424‑M1/1‑30 (MQ=255)
                                                            gatgatTCCCGTCTTCACCTGCCTGTTTTGgaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatg                               >  2:60907‑M1/1‑30 (MQ=255)
                                                            gatgatTCCCGTCTTCACCTGCCTGTTTTGgaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatg                               >  2:286017‑M1/1‑30 (MQ=255)
                                                            gatgatTCCCGTCTTCACCTGCCTGTTTTGgaaggtgcgaataagcggggaaat                                                                   <  2:5700‑M1/54‑25 (MQ=255)
                                                            gatgatTCCCGTCTTCACCTGCCTGTTTTGgaaggtgcgaataagcggggaaat                                                                   >  1:5700‑M1/1‑30 (MQ=255)
                                                                           cACCTGCCTGTTTTGgaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttt                >  2:430512‑M1/1‑15 (MQ=255)
                                                                           cACCTGCCTGTTTTGgaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttt                >  2:174822‑M1/1‑15 (MQ=255)
                                                                             cctgcctgTTTTGgaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttc              >  1:679162‑M1/1‑13 (MQ=255)
                                                                                gcctgTTTTGgaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcc           >  2:107881‑M1/1‑10 (MQ=255)
                                                                                    gTTTTGgaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgta       >  2:84932‑M1/1‑6 (MQ=255)
                                                                                         ggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgc  >  2:535519‑M1/1‑1 (MQ=255)
                                                                                         ggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgc  >  1:10964‑M1/1‑1 (MQ=255)
                                                                                         |                                                                                         
AAAAGTCTGGCTACTGGCACTGGTATGGAGACGTACCTGTTCGTTAGCCAACAGAATACCGATGATTCCCGTCTTCACCTGCCTGTTTTGCCTTTTCCTGGCATCGTTTCTTGCGGAAATGGCGCAGACCGACATGGGATATTTTATTCCTTCGCGAGTTCTGGTCGGTTTGGGAGCGG  >  NZ_CP009273/3629016‑3629194

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