New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 40362628 (0.300)5 (0.060) 5/182 NT 17.4% intergenic (+25/‑77) psiF/yaiC PsiF family protein/diguanylate cyclase, cellulose regualtor
?NC_000913 = 403636 22 (0.260)intergenic (+35/‑67) psiF/yaiC PsiF family protein/diguanylate cyclase, cellulose regualtor

TTGACGCCACAGCAGCAAAAGATGCGCGAATGCAATAATCAAGCAACACAACAATCTCTGAAAGGTGATGATCGTAATAAGTTTATGAGTGCCTGCCTCAAGAAAGCCGCCTGATACCTGATAGTGCTAACGGGTGAGCTACGAAAATGGCTCACCCGAAATATCATACTTCTGCCTTTAGCTCCGTCTCTATAATTTGGGAAAATTGTTTCTGAATGTTCCCAAAAATAATGAATGATGAAAACTTTTTCAAAAAAGCGGCGGCGCACGGGGAGGAACCTC  >  NC_000913/403488‑403769
                                                                                                                                          |                                                                                                                                               
ttGACGCCACAGCAGCAAAAGATGCGCGAATGCAATAATCAAGCAACACAACAATCTCTGAAAGGTGATGATCGTAATAAGTTTATGAGTGCCTGCCTCAAGAAAGCCGCCTGATACCTGATAGTGCTAACGGGTGAGCTACGAAAAt                                                                                                                                        <  3:220835/148‑1 (MQ=255)
ttGACGCCACAGCAGCAAAAGATGCGCGAATGCAATAATCAAGCAACACAACAATCTCTGAAAGGTGATGATCGTAATAAGTTTATGAGTGCCTGCCTCAAGAAAGCCGCCTGATACCTGATAGTGCTAACGGGTGAGCTACGAAAAt                                                                                                                                        >  4:220835/1‑148 (MQ=255)
      ccACAGCAGCAAAAGATGCGCGAATGCAATAATCAAGCAACACAACAATCTCTGAAAGGTGATGATCGTAATAAGTTTATGAGTGCCTGCCTCAAGAAAGCCGCCTGATACCTGATAGTGCTAACGGGTGAGCcattttcgtagctca                                                                                                                                  <  4:857646‑M1/148‑16 (MQ=255)
                aaaaGATGCGCGAATGCAATAGTCAAGCAACACAACAATCTCTGAAAGGTGATGATCGTAATAAGTTTATGAGTGCCTGCCTCAAGAAAGCCGCCTGATACCTGATAGTGCTAACGGGTGAGCTACGAAAATGGCTCACCCGAAatat                                                                                                                        >  2:1015265/1‑148 (MQ=255)
                             aTGCAATAATCAAGCAACACAACAATCTCTGAAAGGTGATGATCGTAATAAGTTTATGAGTGCCTGCCTCAAGAAAGCCGCCTGATACCTGATAGTGCTAACGGGTGAGCTACGAAAATg                                                                                                                                       >  3:810276/1‑120 (MQ=255)
                             aTGCAATAATCAAGCAACACAACAATCTCTGAAAGGTGATGATCGTAATAAGTTTATGAGTGCCTGCCTCAAGAAAGCCGCCTGATACCTGATAGTGCTAACGGGTGAGCTACGAAAATg                                                                                                                                       <  4:810276/120‑1 (MQ=255)
                                   tagtCAAGCAACACAACAATCTCTGAAAGGTGATGATCGTAATAAGTTTATGAGTGCCTGCCTCAAGAAAGCCGCCTGATACCTGATAGTGCTAACGGGTGAGCTACGAAAATGGCTCACCCGAAATATCATACTTCTGCCTTTAGCt                                                                                                     <  1:1015265/145‑1 (MQ=255)
                                      tCAAGCAACACAACAATCTCTGAAAGGTGATGATCGTAATAAGTTTATGAGTGCCTGCCTCAAGAAAGCCGCCTGATACCTGATAGTGCTAACGGGTGAGCTACGAAAAt                                                                                                                                        <  1:148388/110‑1 (MQ=255)
                                      tCAAGCAACACAACAATCTCTGAAAGGTGATGATCGTAATAAGTTTATGAGTGCCTGCCTCAAGAAAGCCGCCTGATACCTGATAGTGCTAACGGGTGAGCTACGAAAAt                                                                                                                                        >  2:148388/1‑110 (MQ=255)
                                               acaacaATCTCTGAAAGGTGATGATCGTAATAAGTTTATGAGTGCCTGCCTCAAGAAAGCCGCCTGATACCTGATAGTGCTAACGGGTGAGCTACGAAAATGgctc                                                                                                                                   >  2:90328‑M2/1‑102 (MQ=255)
                                               acaacaATCTCTGAAAGGTGATGATCGTAATAAGTTTATGAGTGCCTGCCTCAAGAAAGCCGCCTGATACCTGATAGTGCTAACGGGTGAGCTACGAAAATGgctc                                                                                                                                   <  1:90328‑M2/106‑5 (MQ=255)
                                                    aaTCTCTGAAAGGTGATGATCGTAATAAGTTTATGAGTGCCTGCCTCAAGAAAGCCGCCTGATACCTGATAGTGCTAACGGGTGAGCTACGAAAAt                                                                                                                                        <  3:189674/96‑1 (MQ=255)
                                                    aaTCTCTGAAAGGTGATGATCGTAATAAGTTTATGAGTGCCTGCCTCAAGAAAGCCGCCTGATACCTGATAGTGCTAACGGGTGAGCTACGAAAAt                                                                                                                                        >  4:189674/1‑96 (MQ=255)
                                                     aTCTCTGAAAGGTGATGATCGTAATAAGTTTATGAGTGCCTGCCTCAAGAAAGCCGCCTGATACCTGATAGTGCTAACGGGTGAGCTa                                                                                                                                               <  2:1095442/88‑1 (MQ=255)
                                                     aTCTCTGAAAGGTGATGATCGTAATAAGTTTATGAGTGCCTGCCTCAAGAAAGCCGCCTGATACCTGATAGTGCTAACGGGTGAGCTa                                                                                                                                               >  1:1095442/1‑88 (MQ=255)
                                                                      atCGTAATAAGTTTATGAGTGCCTGCCTCAAGAAAGCCGCCTGATACCTGATAGTGCTAACGGGTGAGCTACGAAAATGGCTCACCCGAAATATCATACTTCTGCCTTTAGCTCCGTCTCta                                                                                            >  2:1097634/1‑122 (MQ=255)
                                                                      atCGTAATAAGTTTATGAGTGCCTGCCTCAAGAAAGCCGCCTGATACCTGATAGTGCTAACGGGTGAGCTACGAAAATGGCTCACCCGAAATATCATACTTCTGCCTTTAGCTCCGTCTCta                                                                                            <  1:1097634/122‑1 (MQ=255)
                                                                                gTTTATGAGTGCCTGCCTCAAGAAAGCCGCCTGATACCTGATAGTGCTAACGGGTGAGCTACGAAAATGGCTCACCCGAAATATCATACTTCTg                                                                                                              >  4:907771/1‑94 (MQ=255)
                                                                                gTTTATGAGTGCCTGCCTCAAGAAAGCCGCCTGATACCTGATAGTGCTAACGAGTGAGCTACGAAAATGGCTCACCCGAAATATCATACTTCTg                                                                                                              <  3:907771/94‑1 (MQ=255)
                                                                                         tgcctgccTCAAGAAAGCCGCCTGATACCTGATAGTGCTAACGGGTGAGCTACGAAAATGGCTCACCCGAAATATCATACTTCTGCCTTTAGCTc                                                                                                    >  2:902469/1‑95 (MQ=255)
                                                                                         tgcctgccTCAAGAAAGCCGCCTGATACCTGATAGTGCTAACGGGTGAGCTACGAAAATGGCTCACCCGAAATATCATACTTCTGCCTTTAGCTc                                                                                                    <  1:902469/95‑1 (MQ=255)
                                                                                                 tCAAGAAAGCCGCCTGATACCTGATAGTGCTAACGGGTGAGCcattttcgtagctc                                                                                                                                   >  1:964571‑M1/1‑42 (MQ=255)
                                                                                                 tCAAGAAAGCCGCCTGATACCTGATAGTGCTAACGGGTGAGCcattttcgtagctc                                                                                                                                   <  2:964571‑M1/56‑15 (MQ=255)
                                                                                                  cAAGAAAGCCGCCTGATACCTGATAGTGCTAACGGGTGAGCcattttcgtagctcacc                                                                                                                                <  4:1057909‑M1/58‑18 (MQ=255)
                                                                                                  cAAGAAAGCCGCCTGATACCTGATAGTGCTAACGGGTGAGCcattttcgtagctcacc                                                                                                                                >  3:1057909‑M1/1‑41 (MQ=255)
                                                                                                  cAAGAAAGCCGCCTGATACCTGATAGTGCTAACGGGTGAGCTACGAAAATGGCTCACCCGAAATATCATACTTCTGCCTTTAGCTCCg                                                                                                  >  4:962106/1‑88 (MQ=255)
                                                                                                  cAAGAAAGCCGCCTGATACCTGATACTGCTAACGGGTGAGCTACGAAAATGGCTCACCCGAAATATCATACTTCTGCCTTTAGCTCCg                                                                                                  <  3:962106/88‑1 (MQ=255)
                                                                                                      aaaGCCGCCTGATACCTGATAGTGCTAACGGGTGAGCTACGAAAATGGCTCACCCGaa                                                                                                                            <  3:263145/58‑1 (MQ=255)
                                                                                                      aaaGCCGCCTGATACCTGATAGTGCTAACGGGTGAGCTACGAAAATGGCTCACCCGaa                                                                                                                            >  4:263145/1‑58 (MQ=255)
                                                                                                         gccgccTGATACCTGATAGTGCTAACGGGTGAGCTACGAAAATGGCTCACCCGAAATATCATACTTCTGCCTTTAGCTCCGTCTCTATAATTTGGGAAAAt                                                                              >  3:387574/1‑101 (MQ=255)
                                                                                                         gccgccTGATACCTGATAGTGCTAACGGGTGAGCTACGAAAATGGCTCACCCGAAATATCATACTTCTGCCTTTAGCTCCGTCTCTATAATTTGGGAAAAt                                                                              <  4:387574/101‑1 (MQ=255)
                                                                                                          ccgccTGATACCTGATAGTGCTAACGGGTGAGCTACGAAAATGGCTCACCCGAAATATCAt                                                                                                                     <  3:465447/61‑1 (MQ=255)
                                                                                                          ccgccTGATACCTGATAGTGCTAACGGGTGAGCTACGAAAATGGCTCACCCGAAATATCAt                                                                                                                     >  4:465447/1‑61 (MQ=255)
                                                                                                                      tgataGTGCTAACGGGTGAGCTACGAAAATGGCTCACCCGAAATATCATACTTCTGCCTTTAGCTCCGTCTCTATAAtt                                                                                       >  1:429075/1‑79 (MQ=255)
                                                                                                                      tgataGTGCTAACGGGTGAGCTACGAAAATGGCTCACCCGAAATATCATACTTCTGCCTTTAGCTCCGTCTCTATAAtt                                                                                       <  2:429075/79‑1 (MQ=255)
                                                                                                                               tgaCGGGTGAGCTACGAAAATGGCTCACCCGAAATATCATACTTCTGCCTTTAGCTCCGTCTCTATAATTTGGGAAAAt                                                                              <  2:381125/77‑1 (MQ=255)
                                                                                                                               tAACGGGTGAGCTACGAAAATGGCTCACCCGAAATATCATACTTCTGCCTTTAGCTCCGTCTCTATAATTTGGGAAAAt                                                                              >  1:381125/1‑79 (MQ=255)
                                                                                                                                    ggTGAGCTACGAAAATGgctcacccgttagcactatcaggtatcaggcggctttcttg                                                                                              <  4:1057909‑M2/58‑42 (MQ=255)
                                                                                                                                    ggTGAGCTACGAAAATGgctcacccgttagcactatcaggtatcaggcggctttcttg                                                                                              >  3:1057909‑M2/1‑17 (MQ=255)
                                                                                                                                      tGAGCTACGAAAATGgctcacccgttagcactatcaggtatcaggcggctttcttgaggcaggcactcataaacttattacgatcatcacctttcagagattgttgtgttgcttgattattgcattcgcgcatcttttgctgctgtgg  <  4:857646‑M2/148‑134 (MQ=255)
                                                                                                                                       gAGCcattttcgtagctcacccgttagcactatcaggtatcaggcggctttcttgaggcaggcactcataaacttattacgatcatcacctttcagagattgttgt                                           >  2:90328‑M1/1‑4 (MQ=255)
                                                                                                                                       gAGCcattttcgtagctcacccgttagcactatcaggtatcaggcggctttcttgaggcaggcactcataaacttattacgatcatcacctttcagagattgttgt                                           <  1:90328‑M1/106‑103 (MQ=255)
                                                                                                                                       gAGCTACGAAAATGgctcacccgttagcactatcaggtatcaggcggctttcttga                                                                                             >  1:964571‑M2/1‑14 (MQ=255)
                                                                                                                                       gAGCTACGAAAATGgctcacccgttagcactatcaggtatcaggcggctttcttga                                                                                             <  2:964571‑M2/56‑43 (MQ=255)
                                                                                                                                          cTACGAAAATGGCTCACCCGAAATATCATACTTCTGCCTTTAGCTCCGTCTCTATAATTTGGGAAAATTGTTTCTGAATGTTCCCAAAaataa                                                     >  2:643417/1‑93 (MQ=255)
                                                                                                                                          cTACGAAAATGGCTCACCCGAAATATCATACTTCTGCCTTTAGCTCCGTCTCTATAATTTGGGAAAATTGTTTCTGAATGTTCCCAAAaataa                                                     <  1:643417/93‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                               
TTGACGCCACAGCAGCAAAAGATGCGCGAATGCAATAATCAAGCAACACAACAATCTCTGAAAGGTGATGATCGTAATAAGTTTATGAGTGCCTGCCTCAAGAAAGCCGCCTGATACCTGATAGTGCTAACGGGTGAGCTACGAAAATGGCTCACCCGAAATATCATACTTCTGCCTTTAGCTCCGTCTCTATAATTTGGGAAAATTGTTTCTGAATGTTCCCAAAAATAATGAATGATGAAAACTTTTTCAAAAAAGCGGCGGCGCACGGGGAGGAACCTC  >  NC_000913/403488‑403769

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