New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_vch_pgi	1487376 =	23 (1.040)	4 (0.190)	4/272	2.5	16.0%	coding (142/801 nt)	ydcF	DUF218 superfamily protein, SAM‑binding
	?	NC_000913_3_vch_pgi	1487398 =	20 (0.940)					coding (164/801 nt)	ydcF	DUF218 superfamily protein, SAM‑binding
Rejected: Frequency below/above cutoff threshold.

GGAACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGTGGCTGGCGCAGGATGATTTCTCCGGTGAGGTGCCGTATCAGGCCGATTGCGTGATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTG  >  NC_000913_3_vch_pgi/1487229‑1487517                                                                                                                                                    |                                                                                                                                              ggaACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGTGGCTGGCGCAGGATGATTTCTCCGGTGAGGTGCCGTATCAGGCCGATTGCGTGATCCTTGCAGGCAATGCGGTTAt                                                                                                                                              >  1:239544/1‑149 (MQ=255)   gccaggatggtcgcttctgctctgccagtggtgcggatagtgttgtagtgcggatgctgtgcgatggcgctatacaaaaaagttgtcgagtgaccgataccaccactaatcagtaaaggaatttgctgatcgcgggcaatcttacATGc                                                                                                                                            >  1:266508‑M1/146‑149 (MQ=255)                 cctcGTTTGCGACGGTTTTCCCTGCAACTATATATGCCATAAATTTTACCGGAGAGTGGCTGGCGCAGGATGATTTCTCCGGTGAGGGGCCGTATCAGGCCGATTGCGTGATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCg                                                                                                                              <  2:116686/146‑1 (MQ=255)                  cccGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGTGGCTGGCGCAGGATGATTTCTCCGGTGAGGTGCCGTATCAGGCCGATTGCGTGATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCgg                                                                                                                             >  1:368766/1‑149 (MQ=255)                                          aaCTATAGATGGCATAATTGTTATCGGACAGTGGCTGGCGCAGGATGATTTCTCCGGTGAGGTGCCGTATCAGGCCGATTGCGTCATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGc                                                                                                     <  2:232216/149‑1 (MQ=255)                                            cTATAGATGCCATAAATGTTATCGGACAGTGGCTGGCGCAGGATGATTTCTCCGGGGAGGTGCCGTATCAGGCCGATTGCGTGATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCaa                                                                                                   <  2:62586/149‑1 (MQ=255)                                                  aTGCCATAAATGTTATCGGACAGTGGCTGGCGCAGGATGATTTCTCCGGTGAGGTGCCGTATCAGGCCGATTGCGTGATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCt                                                                                             >  2:263336/1‑149 (MQ=255)                                                    gCCATAAATGTTATCGGACAGTGGCTGGCGCAGGATGATTTCTCCGGTGAGGTGCCGTATCAGGCCGATTGCATGATCCTTGCAGGCAATGCGGTTATGCCGACTATAGATGAGGCATGTAAGATTGGCCGCGATCAGCAAGTTACttt                                                                                           >  1:160379/1‑149 (MQ=255)                                                       aTAAATGTTATCGGACAGTGGCTGGCGCAGGATGATTTCTCCGGTGAGGTGCCGTATCAGGCCGATTGCGTGATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACt                                                                                        >  1:363050/1‑149 (MQ=255)                                                              ttATCGGACAGTGGCTGGCGCAGGATGATTTCTCCGGTGAGGTGCCGTATCAGGCCGATTGCGTGATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAgt                                                                                 >  1:391713/1‑149 (MQ=255)                                                                    gACAGTGGCTGGCGCAGGATGATTTCTCCGGTGAGGTGCCGTATCAGGCCGATTGCGTGATCCTTGCAGGCAATGCGGTTATGCCGACTGCACAGTTTGCATGTAAGATTGCGCGCGATCAGCAAATTCCTTTACTGATTAGTggtggt                                                                           >  1:189845/1‑149 (MQ=255)                                                                               aaaaagttgtcgagtgaccgataccaccactaatcagtaaaggaatttgctgatcgcgggcaatcttacATGCCGACTATCGATGCGGCATGTaa                                                                                                                      <  1:273540‑M1/26‑1 (MQ=255)                                                                               aaaaagttgtcgagtgaccgataccaccactaatcagtaaaggaatttgctgatcgcgggcaatcttacATGCCGACTATCGATGCGGCATGTaa                                                                                                                      >  2:273540‑M1/70‑95 (MQ=255)                                                                                 gcAGGATGATTTCTCCGGTGAGGTGCCGTATCAGGCCGATTGCGTGATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGTAAATTCCTTTACTGCTTAGTGGTGGTATCGGTCACTCGa                                                              >  2:218071/1‑149 (MQ=255)                                                                                     gatgatTTCTCCGGTGAGGTGCCGTATCAGGCCGATTGCGTGATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAAc                                                          >  1:106851/1‑149 (MQ=255)                                                                                          tttCTCCGGTGAGGTGCCGTATCAGGCCGATTGCGTGATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACttttt                                                     <  1:263336/149‑1 (MQ=255)                                                                                           ttatcCGGTGAGGTGTCGTCTCTGGGCGATTGCGTGATCCTTTGAGGTAATGCGGTTATGCCGACTATGGATGCGGAATGTAAGATTGTCCGCGATCAGCAAATTCCTTTCCTGATTAGTGGTGGTATCGGTCACTCGACAACtttttt                                                    <  2:156819/146‑1 (MQ=255)                                                                                                  gataccaccactaatcagtaaaggaatttgctgatcgcgggcaatcttacATGCCGACTATCGATGCGGCATGTAAGATTg                                                                                                                 <  1:341571‑M1/31‑1 (MQ=255)                                                                                                  gataccaccactaatcagtaaaggaatttgctgatcgcgggcaatcttacATGCCGACTATCGATGCGGCATGTAAGATTg                                                                                                                 >  2:341571‑M1/51‑81 (MQ=255)                                                                                                  gTGAGGTGCCGTATCAGGCCGATTGCGTGATCCTTGCAGGCAATGCGGTTCTCCCGACTATCGATGCGGCATGCAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACCCGACAACTTTTTTGTATAgc                                             >  2:358557/1‑149 (MQ=255)                                                                                                               tCAGGCCGATTGCGTGATCCTTGCAGGCAATTCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGc                                <  2:106851/149‑1 (MQ=255)                                                                                                                    ccGATTGCGTGATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGACATCGCGCAGCATCCg                           >  1:247304/1‑149 (MQ=255)                                                                                                                      gATTGCGTGATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCa                                                                                                       >  2:122311/1‑71 (MQ=255)                                                                                                                      gATTGCGTGATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCa                                                                                                       <  1:122311/71‑1 (MQ=255)                                                                                                                          gCGTGATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCACGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCTGCACTac                     <  2:239544/149‑1 (MQ=255)                                                                                                                           cGTGATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTaca                    >  1:374116/1‑149 (MQ=255)                                                                                                                                  cTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTAt             <  2:374116/149‑1 (MQ=255)                                                                                                                                             tGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTg  >  1:202537/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                              GGAACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGTGGCTGGCGCAGGATGATTTCTCCGGTGAGGTGCCGTATCAGGCCGATTGCGTGATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTG  >  NC_000913_3_vch_pgi/1487229‑1487517

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.