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

AACAGACCTCGACTTCATGGCCCGCCTGGCCGCCGAGGAAGGCATTGTCTACCGCTTCGACCACACGCCGAAACAGCACTTGCTAATCATCACTGACCGGCTG‑‑CTGGCGCTGGGCCAGATCAGCCGGGGGGCCATCAAGACCGATGATGAGGATGAAGGGTTTCACAATGACGACGTGGGCCCCGATCAGGTGCTGTACCGATCCGTTGGAGGCGGCGATCAGGCCAGGCCTTGCCTGCGGCGCTTGCGTTACAGCG  >  NC_002947/3833831‑3834087                                                                                                                                  |                                                                                                                                  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)                        tcggtacagcacctgatcggggcccacgtcgtcattgtgaaacccttcatcctcatcatcggtcttgatggccccccggctgctggcgctgggccagatcagccggGGGGCCATCAAGACCGAtgatga                                                                                                             >  1:135710‑M2/107‑129 (MQ=255)                        tcggtacagcacctgatcggggcccacgtcgtcattgtgaaacccttcatcctcatcatcggtcttgatggccccccggctgctggcgctgggccagatcagccggGGGGCCATCAAGACCGAtgatga                                                                                                             <  2:135710‑M2/23‑1 (MQ=255)                           tGGCCGCCGAGGAAGGCATTGTCTACCGCTTCGACCACACGCCGAAACAGCACTCGCTAATCATCACTGACCGGCTG‑‑CTGGCGCTGGGCCAGATCAGCCGGGGGGCCATCAAGACCGATGATGAGGATGAAGGGTTTCACAATGacgacg                                                                                   >  2:262572/1‑150 (MQ=255)                             gccgccGAGGAAGGCATTGTCTACCGCTTCGACCACACGCCGAAACAGCACTTGCTAATCATCACTGACCGGCTG‑‑CTGGCGCTGGGCCAGATCAGCCGGGGGGCCATCAAGACCGATGATGAGGATGAAGGGTTTCACAATGACGACGTg                                                                                 <  2:253739/150‑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)                                                                                                     ggctgatctggcccagcgccagcagccggGGGGCCATCAAGACCGATGATGAGGATGAAGGGTTTCACAATGACGACGTGGGCCCCGATCAg                                                                     >  1:16649‑M2/30‑92 (MQ=255)                                                                                                     ggctgatctggcccagcgccagcagccggGGGGCCATCAAGACCGATGATGAGGATGAAGGGTTTCACAATGACGACGTGGGCCCCGATCAg                                                                     <  2:16649‑M2/63‑1 (MQ=255)                                                                                                       G‑‑CTGGCGCTGGGCCAGATCAGCCGGGGGGCCATCAAGACCGATGATGAGGATGAAGGGTTTCACAATGACGACGTGGGCCCCGATCAGGTGCTGTACCGATCCGTTGGAGGCGGCGATCAGGCCAGGCCTTGCCTGCGGCGCTTGCGTTa       >  1:268273‑M1/1‑150 (MQ=255)                                                                                                          ctgGCGCTGGGCCAGATCAGCCGGGGGGCCATCAAGACCGATGATGAGGATGAAGGGTTTCACAATGACGACGTGGGCCCCGATCAGGTGCTGTACCGATCCGTTGGAGGCGGCGATCAGGCCAGGCCTTGCCTGCGGCGCTTGCGTTAc      <  1:442910/150‑1 (MQ=255)                                                                                                              cgcTGGGCCAGATCAGCCGGGGGGCCATCAAGACCGATGATGAGGATGAAGGGTTTCACAATGACGACGTGGGCCCCGATCAGGTGCTGTACCGATCCGTTGGAGGCGGCGATCAGGCCAGGCCTTGCCTGCGGCGCTTGCGTTACAGCg  <  2:418302/150‑1 (MQ=255)                                                                                                                                  |                                                                                                                                  AACAGACCTCGACTTCATGGCCCGCCTGGCCGCCGAGGAAGGCATTGTCTACCGCTTCGACCACACGCCGAAACAGCACTTGCTAATCATCACTGACCGGCTG‑‑CTGGCGCTGGGCCAGATCAGCCGGGGGGCCATCAAGACCGATGATGAGGATGAAGGGTTTCACAATGACGACGTGGGCCCCGATCAGGTGCTGTACCGATCCGTTGGAGGCGGCGATCAGGCCAGGCCTTGCCTGCGGCGCTTGCGTTACAGCG  >  NC_002947/3833831‑3834087

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.