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.

GTACGGTGGTGATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGAT  >  NC_000913/4316942‑4317212
                                                                                                                                             |                                                                                                                                 
gTACGGTGGTGATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCt                                                                                                                            >  2:57297/1‑149 (MQ=255)
            tGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCa                                                                                                                  >  1:208735/1‑147 (MQ=255)
            tGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCa                                                                                                                  <  2:208735/147‑1 (MQ=255)
                         ccaatgggagcctcttcatgattatcaataacgataagctggtgctggaaaacgaggtggtaagcggttcgctggaggtgcagaacggcgaaatccgcgcctttgccgaaagccagAGCCGCCTTCGCCGTCCATCGCCTCCGGCAggc                                                                                                   <  5:163448‑M1/33‑1 (MQ=255)
                           tCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAggcgg                                                                                                 >  5:157257/1‑149 (MQ=255)
                                 agcctcttcatgattatcaataacgttaagctggtgctggaaaacgaggtggtaagcggttcgctggaggtgcagaacggcgaaatccgcgcctttgccgaaagccagAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTgg                                                                                           >  6:163448‑M1/109‑149 (MQ=255)
                                           gctCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCa                                                                                  >  8:133594/1‑148 (MQ=255)
                                                           gCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGcc                                                                 <  1:57297/149‑1 (MQ=255)
                                                                              gaggtggtaagcggttcgctggaggtgcagaacggcgaaatccgcgcctttgccgaaagccagAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCgg                                                                                    <  8:140401‑M1/48‑1 (MQ=255)
                                                                              gaggtggtaagcggttcgctggaggtgcagaacggcgaaatccgcgcctttgccgaaagccagAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCgg                                                                                    >  7:140401‑M1/64‑111 (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                             <  6:157257/148‑1 (MQ=255)
                                                                                                tATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGtt                             <  7:133594/148‑1 (MQ=255)
                                                                                                 aTCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGttt                            <  1:84077/148‑1 (MQ=255)
                                                                                                    cAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTc                          >  6:80660/1‑147 (MQ=255)
                                                                                                     ggATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGttt                            <  8:65810/143‑1 (MQ=255)
                                                                                                     ggATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGttt                            >  7:65810/2‑144 (MQ=255)
                                                                                                     aGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAg                       <  5:80660/149‑1 (MQ=255)
                                                                                                                          aTCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGa   >  2:37246/1‑148 (MQ=255)
                                                                                                                          aTCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGAt  <  8:60909/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                 
GTACGGTGGTGATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGAT  >  NC_000913/4316942‑4317212

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