New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_vch_pgi	= 3867983	31 (0.980)	4 (0.130)	4/282	4.2	11.2%	coding (256/333 nt)	yidQ	DUF1375 family outer membrane protein
	?	NC_000913_3_vch_pgi	= 3868014	33 (1.070)					coding (287/333 nt)	yidQ	DUF1375 family outer membrane protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CCCGCGCCAGCGCAACCATGATTGGCGACGATGAAACCAACTGGGGCACTAAGTCGCTGGCGATTCTCGATATGCCGTTTACCGCTGTGATGGACACGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAGCGATCAGTTTCGGTGATCCACAGCTGGCCGCCTTCCATCCACGCCAGCCACTGACCATCCGGCGAG  >  NC_000913_3_vch_pgi/3867840‑3868131                                                                                                                                                |                                                                                                                                                     cccGCGCCAGCGCAACCATGATTGGCGACGATGAAACCAACTGGGGCACTAAGTCGCTGGCGATTCTCGATATGCCGTTTACCGCTGTGATGGACACGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTcgcgc                                                                                                                                                 <  1:397430/149‑1 (MQ=255) cccGCGCCAGCGCAACCATGATTGGCGACGATGAAACCAACTGGGGCACTAAGTCGCTGGCGATTCTCGATATGCCGTTTACCGCTGTGATGGACACGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTcgcgc                                                                                                                                                 <  1:397405/149‑1 (MQ=255)  ccGCGCCAGCGCAACCATGATTGGCGACGATGAAACCAACTGGGGCACTAAGTCGCTGGCGATTCTCGATATGCCGTTTACCGCTGTGATGGACACGCTTTTACTGCCGTTGGATGTGTTCCGCAAAGACAGCTCGGTGCCCTcgcgcg                                                                                                                                                >  2:267926/1‑149 (MQ=255)   cgcgCCAGCGCAACCATGATTGGCGACGATGAAACCAACTGGGGCACTAAGTCGCTGGCGATTCTCGATATGCCGTTTACCGCTGTGATGGACACGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGt                                                                                                                                               >  1:78601/1‑149 (MQ=255)           cgcAACCATGATTGGCGACGATGAAACCAACTGGGGCACTAAGTCGCTGGCGATTCTCGATATGCCGTTTACCGCTGTGATGGACACGCTTTTACTGCCGTGGGATGTGTTCCTCAAAGACCGCTCGGTGCGCTCGCGCGTTGaaaaaa                                                                                                                                       <  2:283873/149‑1 (MQ=255)                        ggCGACGATGAAACCAACTGGGGCACTAAGTCGCTGGCGATTCTCGATATGCCGTTTACCGCTGTGATGGACACGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAAcgc                                                                                                                          <  1:342126/149‑1 (MQ=255)                        ggCGACGAAGAAACCCCCTGGGGCACTAAGCCGCCGGCGATACTCGATATGCCGTTTACCGCGGTGATGGACACGCTTTTACTGCCGTAGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAAcgc                                                                                                                          <  1:140439/149‑1 (MQ=255)                           gacgaTGAAACCAACTGGGGCACTAAGTCGCTGGCGATTCTCGATATGCCGTTTACCGCTGTGATGGACACGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGGGAAGCCAACGCGCa                                                                                                                       >  1:463419/1‑149 (MQ=255)                                    aCCAACTGGGGCACTAAGTCGCTGGCGATTCTCGATATGCCGTTTACCGCTGTGGTGGACACGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCaa                                                                                                              >  1:419127/1‑149 (MQ=255)                                     ccAACTGGGGCACTAAGTCGCTGGCGATTCTCGATATGCCGTTTACCGCTGTGATGGACACGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAAt                                                                                                             >  1:528940/1‑149 (MQ=255)                                         cTGGGGAACTAAGTCGGTGCCGATTCTCGATATGCCGTTTACCGCTGTGATGGACACGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGTGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCgg                                                                                                         <  2:370491/149‑1 (MQ=255)                                            gggCACTAAGTCGCTGGCGATTCTCGATATGCCGTTTACCGCTGTGATGGACACGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCCCCAATGCGGttc                                                                                                      >  2:287212/1‑148 (MQ=255)                                            gggCACTAAGTCGCTGGCGATTCTCGATATGCCGTTTACCGCTGTGATGGACACGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTa                                                                                                      >  2:350034/1‑149 (MQ=255)                                              gCACTAAGTCGCTGGCGATTCTCGATATGCCGTTTACCGCTGTGATGGACACGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTAtt                                                                                                    >  1:315540/1‑149 (MQ=255)                                                  tAAGTCGCTGGCGATTCTCGATATGCCGTTTACCGCTGTGATGGACACGCTTTTACTGCCGTGGGGTGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGGTATTCCAc                                                                                                >  1:309449/1‑149 (MQ=255)                                                         cTGGCGATTCTCGATATGCCGTTTACCGCTGTGATGGACACGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAg                                                                                         >  1:292068/1‑149 (MQ=255)                                                                  ctcGATATGCCGTTTACCGCTGTGATGGACACGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTACGCTCGCGCGATGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCAAGAATGCCTGa                                                                                >  1:232991/1‑149 (MQ=255)                                                                                    gCTGTGATGGACACGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAGCGATCa                                                              >  2:390450/1‑149 (MQ=255)                                                                                                  gCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAGCGATCAGTTTCGGTGATCca                                                >  1:225339/1‑149 (MQ=255)                                                                                                   cTTTTACTGCCGTGGGATGTGTTCCGCAAAGAAAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAGCGATCAGTTTCGGTGATCcac                                               >  2:535339/1‑149 (MQ=255)                                                                                                             cGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAGCGATCAGTTTCGGTGATCCACAGCTGgccgc                                     >  1:391031/1‑149 (MQ=255)                                                                                                               tGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAGCGATCAGTTTCGGTGATCCACAGCTGGCCGCCt                                   <  1:353465/149‑1 (MQ=255)                                                                                                                       gtTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAGCGATCAGTTTCGGTGATCCACAGCTGGCCGCCTtccatcca                           >  2:519904/1‑149 (MQ=255)                                                                                                                            gCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCCCTGCGGTTATAACGCCAACCAGAATGCCTGATAACTGATTTAGCGATCAGTTTCGGGGATCCCCAGCTGGCCGCCTTCCATCCACgccc                      >  1:216116/1‑148 (MQ=255)                                                                                                                            gCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAGCGATCAGTTTCGGTGATCCACAGCTGGCCGCCTTCCATCCACgcca                      <  1:20532/149‑1 (MQ=255)                                                                                                                            gCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAGCGATCAGTTTCGGTGATCCACAGCTGGCCGCCTTCCATCCACgcca                      >  2:501019/1‑149 (MQ=255)                                                                                                                             cAAAGACAGCTGGGTGAGCTCGTGCGTTGAAGAAAGCGAAGCCAACGCGCAAGGCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAGCGATCAGTTTCGGTGATCCACAGCTGGCCGCCTTCCATCCACgccag                     <  2:227850/149‑1 (MQ=255)                                                                                                                                    aGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAGCGATCAGTTTCGGTGATCCACAGCTGGCCGCCTTCCATCCACGCCAGCCACTGa              >  2:456867/1‑149 (MQ=255)                                                                                                                                    aGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAGCGATCAGTTTCGGTGATCCACAGCTGGCCGCCTTCCATCCACGCCAGCCACTGa              <  1:519904/149‑1 (MQ=255)                                                                                                                                    aGCGCGGTGCGCTCGTGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAGCGATCAGTTTCGGTGATCCACAGCTGGCCGCCTTCCATCCACGCCAGCCACTGa              <  1:514758/149‑1 (MQ=255)                                                                                                                                                tCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAGCGATCAGTTTCGGTGATCCACAGCTGGCCGCCTTCCATCCACGCCAGCCACTGACCATCCGGCgag  >  1:542994/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                     CCCGCGCCAGCGCAACCATGATTGGCGACGATGAAACCAACTGGGGCACTAAGTCGCTGGCGATTCTCGATATGCCGTTTACCGCTGTGATGGACACGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAGCGATCAGTTTCGGTGATCCACAGCTGGCCGCCTTCCATCCACGCCAGCCACTGACCATCCGGCGAG  >  NC_000913_3_vch_pgi/3867840‑3868131

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