New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4317083 =	19 (0.520)	4 (0.120)	4/224	2.8	17.1%	coding (136/1137 nt)	phnM	ribophosphonate triphosphate hydrolase
	?	NC_000913	4317120 =	21 (0.620)					coding (99/1137 nt)	phnM	ribophosphonate triphosphate hydrolase
Rejected: Frequency below/above cutoff threshold.

CGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTC  >  NC_000913/4316975‑4317256                                                                                                                                                  |                                                                                                                                         agcctcttcatgattatcaataacgttaagctggtgctggaaaacgaggtggtaagcggttcgctggaggtgcagaacggcgaaatccgcgcctttgccgaaagccagAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTgg                                                                                                                                       >  6:163448‑M1/109‑149 (MQ=255)           gctCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCa                                                                                                                              >  8:133594/1‑148 (MQ=255)                           gCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGcc                                                                                                             <  1:57297/149‑1 (MQ=255)                                              gaggtggtaagcggttcgctggaggtgcagaacggcgaaatccgcgcctttgccgaaagccagAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCgg                                                                                                                                >  7:140401‑M1/64‑111 (MQ=255)                                              gaggtggtaagcggttcgctggaggtgcagaacggcgaaatccgcgcctttgccgaaagccagAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCgg                                                                                                                                <  8:140401‑M1/48‑1 (MQ=255)                                                        gaagaaTTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTacca                                                                                >  1:128265/1‑149 (MQ=255)                                                              tttATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGcc                                                                                                             <  3:69110/114‑1 (MQ=255)                                                              tttATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGcc                                                                                                             >  4:69110/1‑114 (MQ=255)                                                              tttATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGt                                                                          >  7:60909/1‑149 (MQ=255)                                                                tATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGtt                                                                         <  7:133594/148‑1 (MQ=255)                                                                tATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGtt                                                                         <  6:157257/148‑1 (MQ=255)                                                                 aTCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGttt                                                                        <  1:84077/148‑1 (MQ=255)                                                                    cAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTc                                                                      >  6:80660/1‑147 (MQ=255)                                                                     ggATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGttt                                                                        >  7:65810/2‑144 (MQ=255)                                                                     ggATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGttt                                                                        <  8:65810/143‑1 (MQ=255)                                                                     aGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAg                                                                   <  5:80660/149‑1 (MQ=255)                                                                                          aTCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGa                                               >  2:37246/1‑148 (MQ=255)                                                                                          aTCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGAt                                              <  8:60909/149‑1 (MQ=255)                                                                                                  ccgaaagccagagccgcctgccggaggcgatggacggcgaaggcggctCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTc                                                                            <  8:140401‑M2/63‑1 (MQ=255)                                                                                                  ccgaaagccagagccgcctgccggaggcgatggacggcgaaggcggctCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTc                                                                            >  7:140401‑M2/49‑111 (MQ=255)                                                                                                         ccagagccgcctgccggaggcgatggacggcgaaggcggctCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCt                               >  6:163448‑M2/42‑149 (MQ=255)                                                                                                                 gcctgccggaggcgatggacggcgaaggcggctCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTATCGTTATTGATAATCATGAAGAGGCTCCCATTgg                       <  5:163448‑M2/116‑1 (MQ=255)                                                                                                                 GCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTgg                       >  2:128265‑M1/1‑149 (MQ=255)                                                                                                                          gTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCg              >  3:170153/1‑149 (MQ=255)                                                                                                                               tCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAg                 <  4:100219/141‑1 (MQ=255)                                                                                                                               tCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAg                 >  3:100219/1‑141 (MQ=255)                                                                                                                                   cTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTgg                       >  8:26499/1‑131 (MQ=255)                                                                                                                                        gCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTc  >  7:8382/1‑147 (MQ=255)                                                                                                                                        gCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTc  <  8:8382/147‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                         CGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTC  >  NC_000913/4316975‑4317256

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