New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? W3110S.gb 2018571 =2 (0.050)33 (0.920) 21/92 NT 94.6% coding (567/687 nt) fliH flagellar biosynthesis protein
?W3110S.gb 2519723 = NA (NA)intergenic (+33/‑54) nupC/insL nucleoside (except guanosine) transporter/predicted transposase

CTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTC  >  W3110S.gb/2519658‑2519787
                                                                 |                                                                
cTGGCTGCGTCAATCGCAGCACTGGTGCTGTAATACCATACATAAAAAAGCCGGGGATAATTCCCATAAg                                                              <  1:675474/70‑1 (MQ=255)
  gTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAgcgc                                                           >  1:628056/1‑71 (MQ=255)
    cTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCt                                                          >  1:343371/1‑70 (MQ=255)
     tGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTaa                                                        >  1:1354864/1‑71 (MQ=255)
     tGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTaa                                                        >  1:505551/1‑71 (MQ=255)
     tGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTaa                                                        >  1:2453840/1‑71 (MQ=255)
            aTCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAATCCCATAAgct                                                            <  1:2647388/60‑2 (MQ=255)
               gcagcaCTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTa                                                    >  1:2031791/1‑65 (MQ=255)
                cagcaCTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTtg                                             <  1:1916027/71‑1 (MQ=255)
                        gTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCAATAAGCGCTAACTTaa                                                   >  1:356046/1‑57 (MQ=255)
                              tAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCt                                >  1:3148073/1‑70 (MQ=255)
                                    catacataAAAAAGCCGGGGATAATTCCCATAAGCGCTAACtt                                                     <  1:3141579/43‑1 (MQ=255)
                                     atacataAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGa                         <  1:1760714/70‑1 (MQ=255)
                                     atacataAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGa                         <  1:1442533/70‑1 (MQ=255)
                                       acataAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAtt                       >  1:3495547/1‑70 (MQ=255)
                                        cataAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTa                                                    >  1:2267282/1‑40 (MQ=255)
                                          taAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGa                              >  1:3682176/1‑60 (MQ=255)
                                                   cGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATg          <  1:1941745/71‑1 (MQ=39)
                                                   cGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATg          <  1:2085932/71‑1 (MQ=39)
                                                     gggATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTaa                    <  1:1290718/59‑1 (MQ=255)
                                                     gggATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGaa        >  1:2732417/1‑71 (MQ=39)
                                                      ggATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCt                                <  1:3133714/46‑1 (MQ=255)
                                                      ggATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCt                                <  1:3143252/46‑1 (MQ=255)
                                                       gATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTa                                        >  1:2907334/1‑37 (MQ=255)
                                                         tAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACg                  >  1:1046023/1‑57 (MQ=38)
                                                           aTTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATg                          >  1:213060/1‑47 (MQ=35)
                                                           aTTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActc  >  1:1865375/1‑71 (MQ=32)
                                                                 |                                                                
CTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTC  >  W3110S.gb/2519658‑2519787

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