New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 3607987	18 (0.930)	4 (0.220)	4/244	NT	17.0%	coding (1173/2265 nt)	PP_3182	GGDEF domain‑containing protein
	?	NC_002947	= 3608006	22 (1.200)					coding (1154/2265 nt)	PP_3182	GGDEF domain‑containing protein

GTGGCCGACAGCTTCTGCACCAGTTCGTGGACCCAGGCGTTGATCTGCTCCTCGCTGCGCTCACCGGCCAGCAACACCACGAATTCGTCGCCGCCAAACCGCGAAGCCTCGTCGCCGGGCTCCAGCAGGCGCTGGATACGCCCGGCCACAGCCTGTAGCAGCAGGTCACCAATCTTGTGGCCGAGCGAGTCGTTGATCGACTTGAAACGGTCCATGTCGATGAACAGGATCGCCATCAGTCGGCGCTTGCCACGTTGCCGAGT  >  NC_002947/3607868‑3608130                                                                                                                        |                                                                                                                                                gTGGCCGACAGCTTCTGCACCAGTTCGTGGACCCAGGCGTTGATCTGCTCCTCGCTGCGCTCACCGGCCAGCAACACCACGAATTCGTCGCCGCCAAACCGCGAAGCCTCGTCGCCGGGCTCCAGCAGGCGCTGGATACGCCCGGCcaca                                                                                                                   <  2:454144/150‑1 (MQ=255)      cGACAGCTTCTGCACCAGTTCGTGGACCCAGGCGTTGATCTGCTCCTCGCTGCGCTCACCGGCCAGCAACACCGCGAATTCGTCGCCGCCAAACCGCGAAGCCTCGTCGCCGGGCTCCAGCAGGCGCTGGATACGCCCGGCCACAGCCTg                                                                                                              >  2:425064/1‑150 (MQ=255)      cGACAGCTTCTGCACCAGTTCGTGGACCCAGGCGTTGATCTGCTCCTCGCTGCGCTCACCGGCCAGCAACACCACGAATTCGTCGCCGCCAAACCGCGAAGCCTCGTCGCCGGGCgtatccagcgcctgctggagcccggcgac                                                                                                                    >  1:202494‑M1/1‑115 (MQ=255)      cGACAGCTTCTGCACCAGTTCGTGGACCCAGGCGTTGATCTGCTCCTCGCTGCGCTCACCGGCCAGCAACACCACGAATTCGTCGCCGCCAAACCGCGAAGCCTCGTCGCCGGGCgtatccagcgcctgctggagcccggcgac                                                                                                                    <  2:202494‑M1/144‑30 (MQ=255)           gCTTCTGCACCAGTTCGTGGACCCAGGTGTTGATCTGCTCCTCGCTGCGCTCACCGGCCAGCAACACCACGAATTCGTCGCCGCCAAACCGCGAAGCCTCGTCGCCGGGCTCCAGCAGGCGCTGGATACGCCCGGCCACAGCCTGTagca                                                                                                         <  2:51552/150‑1 (MQ=255)               cTGCACCAGTTCGTGGACCCAGGCGTTGATCTGCTCCTCGCTGCGCTCACCGGCCAGCAACACCACGAATTCGTCGCCGCCAAACCGCGAAGCCTCGTCGCCGGGCTCCAGCAGGCGCTGGATACGCCCGGCCACAGCCTGTAGcagcag                                                                                                     <  2:208978/150‑1 (MQ=255)                tGCACCAGTTCGTGGACCCAGGCGTTGATCTGCTCCTCGCTGCGCTCACCGGCCAGCAACACCACGAATTCGTCGCCGCCAAACCGCGAAGCCTCGTCGCCGGGCTCCAGCAGGCGCTGGATACGCCCGGCCACAGCCTGTAGCAGCAgg                                                                                                    <  1:425064/150‑1 (MQ=255)                        ttCGTGGACCCAGGCGTTGAGCTGCTCCTCGCTGCGCTCACCGGCCAGCAACACCACGAATTCGTCGCCGCCAAACCGCGAAGCCTCGTCGCCGGGCTCCAGCAGGCGCTGGATACGCCCGGCCACAGCCTGTAGCAGCAGGTCACCAAt                                                                                            >  2:15878/1‑150 (MQ=255)                                         tGATCTGCTCCTCGCTGCGCTCACCGGCCAGCAACACCACGAATTCGTCGCCGCCAAACCGCGAAGCCTCGTCGCCGGGCgtatccagcgcctgctggagcccggc                                                                                                                       <  1:121530‑M1/106‑27 (MQ=255)                                         tGATCTGCTCCTCGCTGCGCTCACCGGCCAGCAACACCACGAATTCGTCGCCGCCAAACCGCGAAGCCTCGTCGCCGGGCgtatccagcgcctgctggagcccggc                                                                                                                       >  2:121530‑M1/1‑80 (MQ=255)                                          gATCTGCTCCTCGCTGCGCTCACCGGCCAGCAACACCACGAATTCGTCGCCGCCAAACCGCGAAGCCTCGTCGCCGGGCTCCAGCAGGCGCTGGATACGCCCGGCCACAGCCTGTAGCAGCAGGTCACCAATCTTGTGGCCGAGCGAGTc                                                                          >  1:114011/1‑150 (MQ=255)                                              tGCTCCTCGCTGCGCTCACCGGCCAGCAACACCACGAATTCGTCGCCGCCAAACCGCGAAGCCTCGTCGCCGGGCTCCAGCAGGCGCTGGATACGCCCGGCCACAGCCTGTAGCAGCAGGTCACCAATCTTGTGGCCGAGCGAGTCGTTg                                                                      <  2:114011/150‑1 (MQ=255)                                                ctcctcGCTGCGCTCACCGGCCAGCAACACAACGAATTCGTCGCCGCCAAACCGCGAAGCCTCGTCGCCAGGCTCCAGCAGGCGCTGGATACGCCCGGCCACAGCCGCGGGCAGCAAGTCACCAATCTTGTGGCCGAGCGAGGCGTTGAt                                                                    >  2:64672/1‑150 (MQ=255)                                                      gCTGCGCTCACCGGCCAGCAACACCACGAATTCGTCGCCGCCAAACCGCGAAGCCTCGTCGCCGGGCTCCAGCAGGCGCTGGATACGCCCGGCCACAGCCTGTAGCAGCAGGTCACCAATCTTGTGGCCGAGCGAGTCGTTGATCGACtt                                                              <  2:260122/150‑1 (MQ=255)                                                           gcTCACCGGCCAGCAACACCACGAATTCGTCGCCGCCAAACCGCGAAGCCTCGTCGCCGGGCTCCAGCAGGCGCTGGATACGCCCGGCCACAGCCTGTAGCAGCAGGTCACCAATCTTGTGGCCGAGCGAGTCGTTGATCGACTTGAAAc                                                         >  2:124960/1‑150 (MQ=255)                                                              cACCGGCCAGCAACACCACGAATTCGTCGCCGCCAAACCGCGAAGCCTCGTCGCCGGGCTCCAGCAGGCGCTGGATACGCCCGGCCACAGCCGGTAGCAGCAGGTCACCAATCTTGTGCCCGAGCGAGTCGTTGATCGACTTGAAACGGt                                                      <  1:64672/150‑1 (MQ=255)                                                                 cGGCCAGCAACACCACGAATTCGTCGCCGCCAAACCGCGAAGCCTCGTCGCCGGGCTCCAGCAGGCGCTGGATACGCCCGGCCACAGCCTGTAGCAGCAGGTCACCAATCTTGTGGCCGAGCGAGTCGTTGATCGACTTGAAACGGTCCa                                                   >  2:30651/1‑150 (MQ=255)                                                                                                  aCCGCGAAGCCTCGTCGCCGGGCTCCAGCAGGCGCTGGATACGCCCGGCCACAGCCTGTAGCAGCAGGTCACCAATCTTGTGGCCGAGCGAGTCGTTGATCGACTTGAAACGGTCCATGTCGATGAACAGGATCGCCATCAGTCGGCGCt                  >  2:412450/1‑150 (MQ=255)                                                                                                    cgcgAAGCCTCGTCGCCGGGCTCCAGCAGGCGCTGGATACGCCCGGCCACAGCCTGTAGCAGCAGGTCACCAATCTTGTGGCCGAGCGAGTCGTTGATCGACTTGAAACGGTCCATGTCGATGAACAGGATCGCCATCAGTCGGCGCTTg                <  1:124960/150‑1 (MQ=255)                                                                                                      cgAAGCCTCGTCGCCGGGCTCCAGCAGGCGCTGGATACGCCCGGCCACAGCCTGTAGCAGCAGGTCACCAATCTTGTGGCCGAGCGAGTCGTTGATCGACTTGAAACGGTCCATGTCGATGAACAGGATCGCCATCAGTCGGCGCTTGcc              >  1:376024/1‑150 (MQ=255)                                                                                                           ccTCGTCGCCGGGCTCCAGCAGGCGCTGGATACGCCCGGCCACAGCCTGTAGCAGCAGGTCACCAATCTTGTGGCCGAGCGAGTCGTTGATCGACTTGAAACGGTCCATGTCGATGAACAGGATCGCCATCAGTCGGCGCTTGCCACGtt         >  2:403340/1‑150 (MQ=255)                                                                                                               gtcgCCGGGCTCCAGCAGGCGCTGGATACgcccggcgacgaggcttcgcggtttggcggcgacgaattcgtggtgttgctggccggtgagcgcagcgaggagcagatcaacgcctgggtccacgaactggtgcagaagctgtcg           <  2:202494‑M2/144‑116 (MQ=255)                                                                                                               gtcgCCGGGCTCCAGCAGGCGCTGGATACgcccggcgacgaggcttcgcggtttggcggcgacgaattcgtggtgttgctggccggtgagcgcagcgaggagcagatcaacgcctgggtccacgaactggtgcagaagctgtcg           >  1:202494‑M2/1‑29 (MQ=255)                                                                                                                  gCCGGGCTCCAGCAGGCGCTGGATACgcccggcgacgaggcttcgcggtttggcggcgacgaattcgtggtgttgctggccggtgagcgcagcgaggagcagatca                                              >  2:121530‑M2/1‑26 (MQ=255)                                                                                                                  gCCGGGCTCCAGCAGGCGCTGGATACgcccggcgacgaggcttcgcggtttggcggcgacgaattcgtggtgttgctggccggtgagcgcagcgaggagcagatca                                              <  1:121530‑M2/106‑81 (MQ=255)                                                                                                                  gCCGGGCTCCAGCAGGCGCTGGATACGCCCGGCCACAGCCTGTAGCAGCAGGTCACCAATCTTGTGGCCGAGCGAGTCGTTGATCGACTTGAAACGGTCCATGTCGATGAACAGGATCGCCATCAGTCGGCGCTTGCCACGTTGCCGAGt  <  1:412450/150‑1 (MQ=255)                                                                                                                       gCTCCAGCAGGCGCTGGATACGCCCGGCCACAGCCTGTAGCAGCAGGTCACCAATCTTGTGGCCGAGCGAGTCGTTGATCGACTTGAAACGGTCCATGTCGATGAACAGGa                                    <  2:197124/111‑1 (MQ=255)                                                                                                                       gCTCCAGCAGGCGCTGGATACGCCCGGCCACAGCCTGTAGCAGCAGGTCACCAATCTTGTGGCCGAGCGAGACGTTGATCGACTTGAAACGGTCCATGTCGATGAACAGGa                                    >  1:197124/1‑111 (MQ=255)                                                                                                                        |                                                                                                                                                GTGGCCGACAGCTTCTGCACCAGTTCGTGGACCCAGGCGTTGATCTGCTCCTCGCTGCGCTCACCGGCCAGCAACACCACGAATTCGTCGCCGCCAAACCGCGAAGCCTCGTCGCCGGGCTCCAGCAGGCGCTGGATACGCCCGGCCACAGCCTGTAGCAGCAGGTCACCAATCTTGTGGCCGAGCGAGTCGTTGATCGACTTGAAACGGTCCATGTCGATGAACAGGATCGCCATCAGTCGGCGCTTGCCACGTTGCCGAGT  >  NC_002947/3607868‑3608130

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