New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 403626	28 (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.