New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	3833927 =	12 (0.740)	4 (0.260)	4/242	NT	25.5%	coding (549/2178 nt)	vgrG‑II	type VI secretion system protein
	?	NC_002947	3833958 =	12 (0.780)					coding (580/2178 nt)	vgrG‑II	type VI secretion system protein

TCAAGTGTCTCGACCCGCATCAGGTGCGCGAATTTTGCGTGCAGGCCGGTGAAACAGACCTCGACTTCATGGCCCGCCTGGCCGCCGAGGAAGGCATTGTCTACCGCTTCGACCACACGCCGAAACAGCACTTGCTAATCATCACTGACCGGCTG‑‑CTGGCGCTGGGCCAGATCAGCCGGGGGGCCATCAAGACCGATGATGAGGATGAAGGGTTTCACAATGACGACGTGGGCCCCGATCAGGTGCTGTACCGATCCGTTGGAGGC  >  NC_002947/3833779‑3834044                                                                                                                                                     |                                                                                                                        tCAAGTGTCTCGACCCGCATCAGGTGCGCGAATTTTGCGTGCAGGCCGGTGAAACAGACCTCGACTTCATGGCCCGCCTGGCCGCCGAGGAAGGCATTGTCTACCGCTTCGACCACACGCCGAAACAGCACTTGCTAATCATCACTGAcc                                                                                                                        <  1:220572/150‑1 (MQ=255)   aaGTGTCTCGACCCGCATCAGGTGCGCGAATTTTGCGTGCAGGCCGGTGAAACAGACCTCGACTTCATGGCCCGCCTGGCCGCCGAGGAAGGCATTGTCTACCGCTTCGACCACACGCCGAAACAGCACTTGCTAATCATCACTGACCgg                                                                                                                      >  1:253739/1‑150 (MQ=255)                                                     aaCAGACCTCGACTTCATGGCCCGCCTGGCCGCCGAGGAAGGCATTGTCTACCGCTTCAACCACACGCCGAAACAGCACTTGCTAATCATCACTGACCGGCTG‑‑CTGGCGCTGGGCCAGATCAGCCGGGGGGCCATCAAGACCGAtgatga                                                                  >  1:190737/1‑150 (MQ=255)                                                                          ccGCCTGGCCGCCGAGGAAGGCATTGTCTACCGCTTCGACCACACGCCGAAACAGCACTTGCTAATCATCACTGACCGGCTG‑‑CTGGCGCTGGGCCAGATCAGCCggg                                                                                        <  2:41454/107‑1 (MQ=255)                                                                          ccGCCTGGCCGCCGAGGAAGGCATTGTCTACCGCTTCGACCACACGCCGAAACAGCACTTGCTAATCATCACTGACCGGCTG‑‑CTGGCGCTGGGCCAGATCAGCCggg                                                                                        >  1:41454/1‑107 (MQ=255)                                                                          ccGCCTGGCCGCCGAGGAAGGCATTGTCTACCGCTTCGACCACACGCCGAAACAGCACTTGCTAATCATCACTGACCGGCTG‑‑CTGGCGCTGGGCCAGATCAGCCGGGGGGCCATCAAGACCGATGATGAGGATGAAGGGTTTCACAATga                                             >  2:57259/1‑150 (MQ=255)                                                                               tGGCCGCCGAGGAAGGCATTGTCTACCGCTTCGACCACACGCCGAAACAGCACTCGCTAATCATCACTGACCGGCTG‑‑CTGGCGCTGGGCCAGATCAGCCGGGGGGCCATCAAGACCGATGATGAGGATGAAGGGTTTCACAATGacgacg                                        >  2:262572/1‑150 (MQ=255)                                                                                 gccgccGAGGAAGGCATTGTCTACCGCTTCGACCACACGCCGAAACAGCACTTGCTAATCATCACTGACCGGCTG‑‑CTGGCGCTGGGCCAGATCAGCCGGGGGGCCATCAAGACCGATGATGAGGATGAAGGGTTTCACAATGACGACGTg                                      <  2:253739/150‑1 (MQ=255)                                                                                      ctgatcggggcccacgtcgtcattgtgaaacccttcatcctcatcatcggtcttgatggccccCCGGCTG‑‑CTGGCGCTGGGCCAGATCAGcc                                                                                           >  1:16649‑M1/64‑92 (MQ=255)                                                                                      ctgatcggggcccacgtcgtcattgtgaaacccttcatcctcatcatcggtcttgatggccccCCGGCTG‑‑CTGGCGCTGGGCCAGATCAGcc                                                                                           <  2:16649‑M1/29‑1 (MQ=255)                                                                                            aGGCATTGTCTACCGCTTCGACCACACGCCGAAACAGCACTTGCTAATCATCACTGACCGGCTG‑‑CTGGCGCTGGGCCAGATCAGCCGGGGGGCCATCAAGACCGATGATGAGGATGAAGGGTTTCACAATGACGACGTGGGCCCCGATCa                           <  1:13904/150‑1 (MQ=255)                                                                                                                aCCACACGCCGAAACAGCACTTGCTAATCATCACTGACCGGCTG‑‑CTGGCGCTGGGCCAGATCAGCCGGGGGGCCATCAAGACCGATGATGAGGATGAAGGGTTTCACAATGACGACGTGGGCCCCGATCAGGTGCTGTACCGATCCGTTg       <  2:145392/150‑1 (MQ=255)                                                                                                                  cacaCGCCGAAACAGCACTTGCTAATCATCACTGACCGGCTG‑‑CTGGCGCTGGGCCAGATCAGCCGGGGGGCCATCAAGACCGATGATGAGGATGAAGGGTTTCACAATGACGACGTGGGCCCCGATCAGGTGCTGTACCGATCCGTTGGa     <  1:190063/150‑1 (MQ=255)                                                                                                                     acGCCGAAACAGCACTTGCTAATCATCACTGACCGGCTG‑‑CTGGCGCTGGGCCAGATCAGCCGGGGGGCCATCAAGACCGAtgatga                                                                  <  1:354812/86‑1 (MQ=255)                                                                                                                     acGCCGAAACAGCACTTGCTAATCATCACTGACCGGCTG‑‑CTGGCGCTGGGCCAGATCAGCCGGGGGGCCATCAAGACCGAtgatga                                                                  >  2:354812/1‑86 (MQ=255)                                                                                                                     acGCCGAAACAGCACTTGCTAATCATCACTGACCGGCTG‑‑CTGGCGCTGGGCCAGATCAGCCGGGGGGCCATCAAGACCGATGATGAGGATGAAGGGTTTCACAATGACGACGTGGGCCCCGATCAGGTGCTGTACCGATCCGTTGGAggc  >  2:268273/1‑150 (MQ=255)                                                                                                                              tcatcatcggtcttgatggccccCCGGCTGATCTGGCCCAGCGCCAG‑‑CAGCCGGGGGGCCATCAAGACCGATGATGAGGATGAAGGGTTTCACAATGACGACGTGGGCCCCGATCAGGTGCTGTACCGa              <  2:135710‑M1/106‑1 (MQ=255)                                                                                                                              tcatcatcggtcttgatggccccCCGGCTGATCTGGCCCAGCGCCAG‑‑CAGCCGGGGGGCCATCAAGACCGATGATGAGGATGAAGGGTTTCACAATGACGACGTGGGCCCCGATCAGGTGCTGTACCGa              >  1:135710‑M1/24‑129 (MQ=255)                                                                                                                                                     |                                                                                                                        TCAAGTGTCTCGACCCGCATCAGGTGCGCGAATTTTGCGTGCAGGCCGGTGAAACAGACCTCGACTTCATGGCCCGCCTGGCCGCCGAGGAAGGCATTGTCTACCGCTTCGACCACACGCCGAAACAGCACTTGCTAATCATCACTGACCGGCTG‑‑CTGGCGCTGGGCCAGATCAGCCGGGGGGCCATCAAGACCGATGATGAGGATGAAGGGTTTCACAATGACGACGTGGGCCCCGATCAGGTGCTGTACCGATCCGTTGGAGGC  >  NC_002947/3833779‑3834044

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