New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_pgi	1786712 =	11 (0.630)	3 (0.180)	3/264	2.2	21.7%	coding (401/2379 nt)	ppsA	phosphoenolpyruvate synthase
	?	NC_000913_3_pae_pgi	1786747 =	11 (0.650)					coding (366/2379 nt)	ppsA	phosphoenolpyruvate synthase

TACATGTTTCACTGCCACGAGAACGGCGTCAAAACCCTGAACGTTGAGGAAGGTTTCCTGCTGACCGGCAAAAGAAGCGTCCGGCATATCTTCTGCGGTGGCGGAGGAGCGCACCGCAAAAGAGGCGTTTTCGTCATCGGCGGAAAGCTGTGCATAGGCTTCGCGGATGGCGTTTTCCAGCTCAGGCTGGAAGGGAGTGTCGATAATCCACTGGCGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTT  >  NC_000913_3_pae_pgi/1786624‑1786890                                                                                                                            |                                                                                                                                                tACATGTTTCACTGCCACGAGAACGGCGTCAAAACCCTGAACGTTGAGGAAGGTTTCCTGCTGACCGGCAAAAGAAGCGTCCGGCATATCTTCTGCGGTGGCGGAGGAGCGCACCGCAAAAGAGGCGTTTTCGTCATCGGCGGAAAGCtg                                                                                                                       <  2:32907/150‑1 (MQ=255)  tggattatcgacactcccttccagcctgagctggaaaacgccatccgcgaagcctatgcacagctttccgccgatgacgaaaacgccTCTTCTGCGGTGGCGGAGGAGCGCACCGCAAAAGAGGCGTTTTCGTCATCGGCGGAAAGCtgt                                                                                                                      >  4:64986‑M1/88‑150 (MQ=255)      gTTTCACTGCCACGAGAACGGCGTCAAAACCCTGAACGTTGAGGAAGGTTTCCTGCTGACCGGCAAAAGAAGCGTCCGGCATATCTTCTGCGGTGGCGGAGGAGCGCACCGCAAAAGAGGCGTTTTCGTCATCGGCGGAAAGCTGTGCAt                                                                                                                  <  2:183628/150‑1 (MQ=255)                         gcctgagctggaaaacgccatccgcgaagcctatgcacagctttccgccgatgacgaaaacgccTCTTCTGCGGTGGCGGAGGAGCGCACCGCAAAAGAGGCGTTTTCGt                                                                                                                                       <  1:46310‑M1/46‑1 (MQ=255)                         gcctgagctggaaaacgccatccgcgaagcctatgcacagctttccgccgatgacgaaaacgccTCTTCTGCGGTGGCGGAGGAGCGCACCGCAAAAGAGGCGTTTTCGt                                                                                                                                       >  2:46310‑M1/65‑110 (MQ=255)                                                 gAAGGTTTCCTGCTGACCGGCAAAAGAAGCGTCCGGCATATCTTCTGCGGTGGCGGAGGAGCGCACCGCAAAAGAGGCGTTTTCGTCATCGGCGGAAAGCTGTGCATAGGCTTCGCGGATGGCGTTTTCCAGCTCAGGCTGGAAGGGAgt                                                                       >  2:153848/1‑150 (MQ=255)                                                             ctgACCGGCAAAAGAAGCGTCCGGCATATCTTCTGCGGTGGCGGAGGAGCGCACCGCAAAAGAGGCGTTTTCGTCATCGGCGGAAAGCTGTGCATAGGCTTCGCGGATGGCGTTTTCCAGCTCAGGCTGGAAGGGAGTGTCGATAATCCa                                                           <  3:54024/150‑1 (MQ=255)                                                             acagctttccgccgatgacgaaaacgcctcttttgcggtgcgctcctccgccaccgcagaagaGGCGTTTTCGTCATCGGCGGAAAGCTGTGCATAGGCTTCGCGGATGGCGTTTTCCAGCTCAGGCTGGAAGGGAGTGTCGATAATCCa                                                           >  4:64986‑M2/64‑150 (MQ=255)                                                                    gCAAAAGAAGCGTCCGGCATATCTTCTGCGGTGGCGGAGGAGCGCACCGCAAAAGAGGCGTTTTCGTCATCGGCGGAAAGCTGTGCATAGGCTTCGCGGATGGCGTTTTCCAGCTCAGGCTGGAAGGGAGTGTCGATAATCCACTGGCgg                                                    >  2:176612/1‑150 (MQ=255)                                                                      aaaaGAAGCGTCCGGCATATCTTCTGCGGTGGCGGAGGAGCGCACCGCAAAAGAGGCGTTTTCGTCATCGGCGGAAAGCTGTGCATAGGCTTCGCGGATGGCGTTTTCCAGCTCAGGCTGGAAGGGAGTGTCGATAATCCACTGGCGGAt                                                  >  4:33090/1‑150 (MQ=255)                                                                              acgaaaacgcctcttttgcggtgcgctcctccgccaccgcagaagaGGCGTTTTCGTCATCGGCGGAAAGCTGTGCATAGGCTTCGCGGATGGCGTTTTCCAGCTCAGGc                                                                                  >  2:46310‑M2/47‑110 (MQ=255)                                                                              acgaaaacgcctcttttgcggtgcgctcctccgccaccgcagaagaGGCGTTTTCGTCATCGGCGGAAAGCTGTGCATAGGCTTCGCGGATGGCGTTTTCCAGCTCAGGc                                                                                  <  1:46310‑M2/64‑1 (MQ=255)                                                                                  cGGCATATCTTCTGCGGTGGCGGAGGAGCGCACCGCAAAAGAGGCGTTTTCGTCATCGGCGGAAAGCTGTGCATCGGCTTCGCGGATGGCGTTTTCCAGCTCAGGCTGGAAGGGAGTGTCGATAATCCACTGGCGGATTTGCGCGCCCGc                                      >  4:27705/1‑150 (MQ=255)                                                                                     cATATCTTCTGCGGTGGCGGAGGAGCGCACCGCAAAAGAGGCGTTTTCGTCATCGGCGGAAAGCTGTGCATAGGCTTCGCGGATGGCGTTTTCCAGCTCAGGCTGGAAGGGAGTGTCGATAATCCACTGGCgg                                                    >  3:33835/1‑133 (MQ=255)                                                                                     cATATCTTCTGCGGTGGCGGAGGAGCGCACCGCAAAAGAGGCGTTTTCGTCATCGGCGGAAAGCTGTGCATAGGCTTCGCGGATGGCGTTTTCCAGCTCAGGCTGGAAGGGAGTGTCGATAATCCACTGGCgg                                                    <  4:33835/133‑1 (MQ=255)                                                                                                             gcgcACCGCAAAAGAGGCGTTTTCGTCATCGGCGGAAAGCTGTGCATAGGCTTCGCGGATGGCGTTTTCCAGCTCAGGCTGGAAGGGAGTGTCGATAATCCACTGGCGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCatt           >  4:4553/1‑148 (MQ=255)                                                                                                                      aaaaGAGGCGTTTTCGTCATCGGCGGAAAGCTGTGCATAGGCTTCGCGGATGGCGTTTTCCAGCTCAGGCTGGAAGGGAGTGTCGATAATCCACTGGCGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGtttt  <  1:153848/150‑1 (MQ=255)                                                                                                                            |                                                                                                                                                TACATGTTTCACTGCCACGAGAACGGCGTCAAAACCCTGAACGTTGAGGAAGGTTTCCTGCTGACCGGCAAAAGAAGCGTCCGGCATATCTTCTGCGGTGGCGGAGGAGCGCACCGCAAAAGAGGCGTTTTCGTCATCGGCGGAAAGCTGTGCATAGGCTTCGCGGATGGCGTTTTCCAGCTCAGGCTGGAAGGGAGTGTCGATAATCCACTGGCGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTT  >  NC_000913_3_pae_pgi/1786624‑1786890

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.