New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? W3110S.gb = 20185832 (0.110)12 (0.670) 10/96 NT 86.2% coding (579/687 nt) fliH flagellar biosynthesis protein
?W3110S.gb = 2521059 NA (NA)intergenic (+170/+30) insL/yfeA predicted transposase/predicted diguanylate cyclase

GTGTGGATGATATGCTCGGCGCTACCTTAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCCTCCATCCTGGCGGCTGTAAAGTCTCCGCCGATGAAGGCGATCTCGACGCCAGTGTCGCCACTCGCTGG  >  W3110S.gb/2018519‑2018652
                                                                |                                                                     
gcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatGGGCGATCCCACCCTCCAt                                                                 >  1:1648527‑M1/53‑71 (MQ=255)
 cacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatGGGCGATCCCACCCTCCAt                                                                 <  1:780110‑M1/19‑1 (MQ=255)
  acgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatGGGCGATCCCACCCTCCATcc                                                               >  1:315582‑M1/51‑71 (MQ=255)
       aaatcatatcaggcgtaataccacaacccttaagttagcgcttatGGGCGATCCCACCCTCCATCCTggcg                                                          <  1:642323‑M1/26‑1 (MQ=255)
        gatATGCTCGGCGCTACCTTAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaact                                                          <  1:207660‑M1/70‑14 (MQ=255)
             tatcaggcgtaataccacaacccttaagttagcgcttatGGGCGATCCCACCCTCCATCCTGGCGGCTGTa                                                    >  1:955251‑M1/40‑71 (MQ=255)
             tatcaggcgtaataccacaacccttaagttagcgcttatGGGCGATCCCACCCTCCATCCTGGCGGCTGTa                                                    >  1:835108‑M1/40‑71 (MQ=255)
              cTCGGCGCTACCTTAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagc                                                                 >  1:1428831‑M1/1‑51 (MQ=255)
              cTCGGCGCTACCTTAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagc                                                                 >  1:596892‑M1/1‑51 (MQ=255)
                   cgcTACCTTAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaac                                              >  1:1763179‑M1/1‑46 (MQ=255)
                      taataccacaacccttaagttagcgcttatGGGCGATCCCACCCTCCATCCTGGCGGCTGTAAAGTCTccg                                           >  1:925167‑M1/31‑71 (MQ=255)
                       aCCTTAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaag                                                      >  1:1441736‑M1/1‑42 (MQ=255)
                       aCCTTAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataa                                                                   >  1:774261‑M1/1‑42 (MQ=255)
                        ataccacaacccttaagttagcgcttatGGGCGATCCCACCCTCCATCCTGGCGGCTGTAAAGTCTccgc                                          >  1:1397331‑M1/29‑70 (MQ=255)
                        ataccacaacccttaagttagcgcttatGGGCGATCCCACCCTCCATCCTGGCGGCTGTAAAGTCTccgc                                          <  1:102838‑M1/42‑1 (MQ=255)
                         cTTAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatctg                                        <  1:458490‑M1/71‑32 (MQ=255)
                          ttAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatctg                                        >  1:363174‑M1/1‑39 (MQ=255)
                          ttAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgcta                                                             <  1:1739900‑M1/49‑11 (MQ=255)
                          ttAAGTTTGAATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgcta                                                             <  1:732563‑M1/49‑11 (MQ=255)
                           ccacaacccttaagttagcgcttatGGGCGATCCCACCCTCCATCCTGGCGGCTGTAAAGTCTCCGCCGAt                                      >  1:1179381‑M1/26‑71 (MQ=255)
                               aacccttaagttagcgcttatGGGCGATCACACCCTcc                                                                   <  1:1190830‑M1/17‑1 (MQ=255)
                                  gCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagtgttgaaccatctgaag                                     <  1:1542521‑M1/65‑35 (MQ=255)
                                   cttaagttagcgcttatGGGCGATCCCACCATCCATCCTGGCGGCTGTAAAGTCTCCGCCGATGAAGGCGa                              <  1:1763894‑M1/54‑1 (MQ=255)
                                    aTGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaac                                              >  1:447864‑M1/1‑29 (MQ=255)
                                     tggctggcGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatctgaagaatgcga                              >  1:1686541‑M1/1‑28 (MQ=255)
                                     taagttagcgcttatGGGCGATCCCACCCTCCATCCTGGCGGCTGTAAAGTCTCCGCCGATGAAGGCGAtc                            <  1:1055132‑M1/56‑1 (MQ=255)
                                      aagttagcgcttatGGGCGATCCCACCCTCCATCCTGGCGGCTGTAAAGTCTCCGCCGATGAAGGCGAtct                           <  1:944673‑M1/57‑1 (MQ=255)
                                           agcgcttatGGGCGATCCCACCCTCCATCCTGGCGGCTGTAAAGtctc                                             <  1:864957‑M1/39‑1 (MQ=255)
                                            gcgcttatGGGCGATCCCACCCTCCATCCTGGCGGCTGTAAAGTCTCCGCCGATGaa                                   <  1:297003‑M1/49‑1 (MQ=255)
                                               cttatGGGCGATCCCACCCTCCATCCTGGCGGCTGTAAAGTCTCCGCCGAt                                      >  1:427574‑M1/6‑51 (MQ=255)
                                                  atGGGCGATCCCACCCTCCATCCTGGCGGCTGTAAAGTCTCCGCCGa                                       <  1:997915‑M1/45‑1 (MQ=255)
                                                           cccACCCTCCATCCTGGCGGCTGTAAAGTCTccgccg                                        >  1:902774/1‑37 (MQ=255)
                                                           cccACCCTCCATCCTGGCGGCTGTAAAGTCTCCGCCGATGAAGGCGATCTCGACGCCAGTGTCGCCACTCg      >  1:565643/1‑71 (MQ=255)
                                                               ccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcc  >  1:881493‑M1/1‑2 (MQ=255)
                                                               ccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcc  >  1:714250‑M1/1‑2 (MQ=255)
                                                               ccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcc  <  1:1639826‑M1/71‑70 (MQ=255)
                                                               ccataagcgctaacttaagggttgaaccatctgaagaatgc                                >  1:924404‑M1/1‑2 (MQ=255)
                                                                |                                                                     
GTGTGGATGATATGCTCGGCGCTACCTTAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCCTCCATCCTGGCGGCTGTAAAGTCTCCGCCGATGAAGGCGATCTCGACGCCAGTGTCGCCACTCGCTGG  >  W3110S.gb/2018519‑2018652

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.