New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 2873992155 (1.090)6 (0.040) 3/264 NT 3.9% coding (1/606 nt)
coding (1428/1428 nt)
cysC
cysN
adenosine 5'‑phosphosulfate kinase
sulfate adenylyltransferase, subunit 1
?NC_000913 = 2874083 147 (1.090)coding (1337/1428 nt) cysN sulfate adenylyltransferase, subunit 1
Rejected: Frequency below/above cutoff threshold.

AACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2873838‑2873992
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgccatCTGAATTCAGTGCATTCGAACTGGAATTGAATGCTCTGGTTCGTCGCCACTTTCCGCACTGGGGCGCGCGCGATTTGCTGGGGGATAAATAATGGCGCTGCATGACGAAAACGTCGTCTGGCATAGCCATCCGGTCACTGTGCAACAAC  <  NC_000913/2874083‑2873935
                                                                                                                                                                                                                                                                                                                
AACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGC                                                                                                                                                            >  2:1912075/1‑150
AACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGC                                                                                                                                                            >  1:489414/1‑150
     CTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCAT                                                                                                                                                       >  2:2001806/1‑150
     CTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCAT                                                                                                                                                       >  1:2094836/1‑150
     CTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCAT                                                                                                                                                       >  1:1839911/1‑150
                                      GGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATCTGAATTCAGTGCATTCGAACTGGAATTGAATG                                                                                                                      >  2:516466/1‑150
                                      GGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATCTGAATTCAGTGCATTCGAACTGGAATTGAATG                                                                                                                      >  1:213178/1‑150
                                                                                                          CAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATCTGAATTCAGTGCATTCGAACTGGAATTGAATGCTCTGGTTCGTCGCCACTTTCCGCACTGGGGCGCGCGCGATTTGCTGGGGGATAAATAATGGCGCTGC                                                  <  1:516466/150‑1
                                                                                                          CAGTGACCGGATGGCTATGCCAGACGACGATTTCGTCATGCAGCGCCATCTGAATTCAGTGCATTCGAACTGGAATTGAATGCTCTGGTTCGTCGCCACTTTCCGCACTGGGGCGCGCGCGATTTGCTGGGGGATAAATAATGGCGCTGC                                                  <  2:213178/150‑1
                                                                                                               ACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCTGCGCCATCTGAATTCAGTTCATTAGAACTGGAATTGAATGCTCTGGTTCGTCGCCACTTTCCGCACTGGTGCGCGCTCGATTTGCTGGGGGATAAATAATGGCGCTGCATGAC                                             <  2:2201048/150‑1
                                                                                                               ACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCTGCGCCATCTGAATTCAGTGCATTCGAACTGGAATTGAATGCTCTGGTTCGTCGCCACTTTCCGCACTGGGGCGCGCGCGATTTGCTGGGGGATAAATAATGGCGCTGCATGAC                                             <  2:2201049/150‑1
                                                                                                                                                    GCGCCATCAGAATTCAGTGCATTCGAACTGGAAATGAATGCTCTGGTTCGTCGCCACTTTCCGCACTGGGGCGCGCGCGATTTGCTGGGGGATAAATAATGGCGCTGCATGACGAAAACGTCGTCTGGCATAGCCATCCGGTCACTGTGC        <  2:1925689/150‑1
                                                                                                                                                     CGCCATCTGAATTCAGTGCATTCGAACTGGAATTGAATGCTCTGGTTCGTCGCCACTTTCCGCACTGGGGCGCGCGCGATTTGCTGGGGGATAAATAATGGCGCTGCATGACGAAAACGTCGTCTGGCATAGCCATCCGGTCACTGTGCA       >  2:279941/1‑150
                                                                                                                                                     CGCCATCTGAATTCAGTGCATTCGAACTGGAATTGAATGCTCTGGTTCGTCGCCACTTTCCGCACTGGGGCGCGCGCGATTTGCTGGGGGATAAATAATGGCGCTGCATGACGAAAACGTCGTCTGGCATAGCCATCCGGTCACTGTGCA       >  1:961187/1‑150
                                                                                                                                                     CGCCATCTGAATTCAGTGCATTCGAACTGGAATTGAATGCTCTGGTTCGTCGCCACTTTCCGCACTGGGGCGCGCGCGATTTGCTGGGGGATAAATAATGGCGCTGCATGACGAAAACGTCGTCTGGCATAGCCATCCGGACACTGTGCA       >  2:409546/1‑150
                                                                                                                                                     CGCCATCTGAATTCAGTGCATTCGAACTGGAATTGAATGCTCTGGTTCGTCGCCACTTTCCGCACGGGGGCGCGCGCGATTTGCTGGGGGATAAATAATGGCGCTGCATGACGAAAACGTCGTCTGGCATAGCCATCCGGTCACTGTGCA       >  2:494434/1‑150
                                                                                                                                                      GCCATCTGAATTCAGTGCATTCGAACTGGAATTGAATGCTCTGGTTCGTCGCCACTTTCCGCACTGGGGCGCGCGCGATTTGCTGGGGGATAAATAATGGCGCTGCATGACGAAAACGTCGTCTGGCATAGCCATCCGGTCACTGTGCAc      >  1:2352868/1‑149
                                                                                                                                                      GCCATCTGAATTCAGTGCATTCGAACTGGAATTGAATGCTCTGGTTCGTCGCCACTTTCCGCACTGGGGCGCGCGCGATTTGCTGGGGGATAAATAATGGCGCTGCATGACGAAAACGTCGTCTGGCATAGCCATCCGGTCACTGTGCAA      >  1:1163273/1‑150
                                                                                                                                                       CCATCTGAATTCAGTGCATTCGAACTGGAATTGAATGCTCTGGTTCGTCGCCACTTTCCGCACTGGGGCGCGCGCGATTTGCTGGGGGATAAATAATGGCGCTGCATGACGAAAACGTCGTCTGGCATAGCCATCCGGTCACTGTGCAAC     >  1:1061218/1‑150
                                                                                                                                                       CCATCTGAATTCAGTGCATTCGAACTGGAATTGAATGCTCTGGTTCGTCGCCACTTTCCGCACTGGGGCGCGCGCGATTTGCTGGGGGATAAATAATGGCGCTGCATGACGAAAACGTCGTCTGGCATAGCCATCCGGTCACTGTGCAAC     <  1:1071294/150‑1
                                                                                                                                                       CCATCTGAATTCAGTGCATTCGAACTGGAATTGAATGCTCTGGTTCGTCGCCACTTTCCGCACTGGGGCGCGCGCGATTTGCTGGGGGATAAATAATGGCGCTGCATGACGAAAACGTCGTCTGGCATAGCCATCCGGTCACTGTGCAAC     <  1:469599/150‑1
                                                                                                                                                       CCATCTGAATTCAGTGCATTCGAACTGGAATTGAATGCTCTGGTTCGTCGCCACTTTCCGCACTGGGGCGCGCGCGATTTGCTGGGGGATAAATAATGGCGCTGCATGACGAAAACGTCGTCTGGCATAGCCATCCGGTCACTGTGCAAC     <  1:469585/150‑1
                                                                                                                                                       CCATCTGAATTCAGTGCATTCGAACTGGAATTGAATGCTCTGGTTCGTCGCCACTTTCCGCACTGGGGCGCGCGCGATTTGCTGGGGGATAAATAATGGCGCTGCATGACGAAAACGTCGTCTGGCATAGCCATCCGGTCACTGTGCAAC     <  1:2230244/150‑1
                                                                                                                                                       CCATCTGAATTCAGTGCATTCGAACTGGAATTGAATGCTCTGGTTCGTCGCCACTTTCCGCACTGGGGCGCGCGCGATTTGCTGGGGGATAAATAATGGCGCTGCATGACGAAAACGTCGTCTGGCATAGCCATCCGGTCACTGTGCAAC     <  1:1542173/150‑1
                                                                                                                                                        CATCTGAATTCAGTGCATTCGAACTGGAATTGAATGCTCTGGTTCGTCGCCACTTTCCGCACTGGGGCGCGCGCGATTTGCTGGGGGATAAATAATGGCGCTGCATGACGAAAACGTCGTCTGGCATAGCCATCCGGTCACTGTGCAACA    >  1:1426033/1‑150
                                                                                                                                                        CATCTGAATTCAGTGCATTCGAACTGGAATTGAATGCTCTGGTTCGTCGCCACTTTCCGCACTGGGGCGCGCGCGATTTGCTGGGGGATAAATAATGGCGCTGCATGACGAAAACGTCGTCTGGCATAGCCATCCGGTCACTGTGCAACA    >  1:1652631/1‑150
                                                                                                                                                          TCTGAATTCAGTGCATTCGAACTGGAATTGAATGCTCTGGTTCGTCGCCACTTTCCGCACTGGGGCGCGCGCGATTTGCTGGGGGATAAATAATGGCGCTGCATGACGAAAACGTCGTCTGGCATAGCCATCCGGTCACTGTGCAACAAC  >  1:694747/1‑150
                                                                                                                                                                                                                                                                                                                
AACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2873838‑2873992
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgccatCTGAATTCAGTGCATTCGAACTGGAATTGAATGCTCTGGTTCGTCGCCACTTTCCGCACTGGGGCGCGCGCGATTTGCTGGGGGATAAATAATGGCGCTGCATGACGAAAACGTCGTCTGGCATAGCCATCCGGTCACTGTGCAACAAC  <  NC_000913/2874083‑2873935

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