New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 35124490 (0.000)10 (0.810) 10/166 0.3 100% intergenic (‑26/+375) aroK/hofQ shikimate kinase AroK/DNA uptake porin HofQ
?NZ_CP009273 3576788 = NA (NA)coding (391/417 nt) BW25113_RS17825 YrhA family protein

TTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAACCTCATAAGATATTGCGAGAGTCAGACTGAAAAT  >  NZ_CP009273/3512373‑3512534
                                                                            |                                                                                     
ttttCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGggtgatgctgcca                                                                          >  2:57787‑M1/1‑77 (MQ=255)
              tAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGggtgatgctgccaacttactgatttag                                                            <  2:100539‑M1/90‑28 (MQ=255)
                         aGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGggtgatgctgccaacttactgatttagtgtatgatggt                                                 <  1:104860‑M1/90‑39 (MQ=255)
                                 atTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttga                                         >  2:126217‑M1/1‑44 (MQ=255)
                                    taaaacactatcaataagttggagtcattaccGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGa                                      >  2:307560‑M2/33‑90 (MQ=255)
                                    gCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggt                                      <  1:71864‑M1/90‑50 (MQ=255)
                                       aacactatcaataagttggagtcattaccGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAAc                                   <  1:307560‑M2/61‑1 (MQ=255)
                                       aacactatcaataagttggagtcattaccGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAAc                                   <  1:76159‑M2/61‑1 (MQ=255)
                                          actatcaataagttggagtcattaccGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAACCTc                                >  2:359454‑M2/27‑90 (MQ=255)
                                              tcaataagttggagtcattaccGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAACCTCATaa                            <  1:150923‑M2/68‑1 (MQ=255)
                                                     tttCGGTACTACTAAGACTATTCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgt                     <  1:126217‑M1/90‑67 (MQ=255)
                                                       tggagtcattaccGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAACCTCATAAGATATTGCg                   >  1:113215‑M2/14‑90 (MQ=255)
                                                             cattaccGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATg                                            <  1:93281‑M2/52‑1 (MQ=255)
                                                             cattaccGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATg                                            >  2:267785‑M2/8‑59 (MQ=255)
                                                             cattaccGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATg                                            <  1:267785‑M2/52‑1 (MQ=255)
                                                             cattaccGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATg                                            >  2:93281‑M2/8‑59 (MQ=255)
                                                                   aGACTATTCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtc       <  1:276648‑M1/90‑81 (MQ=255)
                                                                    gACTATTCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtcc      >  1:26341‑M1/1‑9 (MQ=255)
                                                                       tATTCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctc   >  1:148182‑M1/1‑6 (MQ=255)
                                                                        aTTCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcc  >  1:293630‑M1/1‑5 (MQ=255)
                                                                            |                                                                                     
TTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAACCTCATAAGATATTGCGAGAGTCAGACTGAAAAT  >  NZ_CP009273/3512373‑3512534

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