New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4467477 =149 (1.160)6 (0.050) 5/162 NT 4.2% coding (47/54 nt) mgtL regulatory leader peptide for mgtA
?NC_000913 4467511 = 138 (1.190)intergenic (+27/‑114) mgtL/mgtA regulatory leader peptide for mgtA/magnesium transporter
Rejected: Frequency below cutoff threshold.

TTCAGGCGTGGCGAAGCCTTACCTGAAGGTAATAAAGCTGAGCAGGGACGTTAAAACGTCACGGCCTTACGCATCCGGTAACGCCGTGAAAAGACGGGCTTACCGGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4467583‑4467477
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttaccggaTGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGA  >  NC_000913/4467511‑4467605
                                                                                                                                                                                                          
TTCAGGCGTGGCGAAGCCTTACCTGAAGGTAATAAAGCTGAGCAGGGACGTTAAAACGTCACGGCCTTACGCATCCGGTAACGCCGTGAAAAGACGGGCT                                                                                                        <  1:912074/100‑1
 TCAGGCGTGGCGAAGCCTTACCTGAAGGTAATAAAGCTGAGCAGGGACGTTAAAACGTCACGGCCTTACGCATCCGGTAACGCCGTGAAAAGACGGGCTT                                                                                                       >  2:1737993/1‑100
 TCAGGCGTGGCGAAGCCTTACCTGAAGGTAATAAAGCTGAGCAGGGACGTTAAAACGTCACGGCCTTACGCATCCGGTAACGCCGTGAAAAGACGGGCTT                                                                                                       >  2:1831962/1‑100
   AGGCGTGGCGAAGCCTTACCTGAAGGTAATAAAGCTGAGCAGGGACGTTAAAACGTCACGGCCTTACGCATCCGGTAACGCCGTGAAAAGACGGGCTTAC                                                                                                     <  2:1050739/100‑1
    GGCGTGGCGAAGCCTTACCTGAAGGTAATAAAGCTGAGCAGGGACGTTAAAACGTCACGGCCTTACGCATCCGGTAACGCCGTGAAAAGACGGGCTTACC                                                                                                    >  2:1076991/1‑100
    GGCGTGGCGAAGCCTTACCTGAAGGTAATAAAGCTGAGCAGGGACGTTAAAACGTCACGGCCTTACGCATCCGGTAACGCCGTGAAAAGACGGGCTTACC                                                                                                    <  2:2535815/100‑1
    GGCGTGGCGAAGCCTTACCTGAAGGTAATAAAGCTGAGCAGGGACGTTAAAACGTCACGGCCTTACGCAGCCGGTAACGCCGGGAAAAGACGGGCTTACC                                                                                                    >  2:1888133/1‑100
                                                  TTAAAACGTCACGGCCTTACGCATCCGGTAACGCCGTCAAAAGACGGGCTTACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATT                                                      >  2:3389156/1‑100
                                                                       CATCCGGTAACGCCGTGAAAAGACGGGCTTACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCC                                 >  2:611224/1‑100
                                                                         TCCGGTAACGCCGTGAAAAGACGGGCTTACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCAC                               >  1:2086254/1‑100
                                                                           CGGTAACGCCGTGAAAAGACGGGCTTACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGC                             >  2:747850/1‑100
                                                                           CGGTAACGCCGTGAAAAGACGGGCTTACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGC                             >  2:3331699/1‑100
                                                                            GGTAACGCCGTGAAAAGACGGGCTTACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCC                            >  2:437710/1‑100
                                                                                                   TTACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTC     <  1:483992/100‑1
                                                                                                   TTACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTC     <  1:2511228/100‑1
                                                                                                    TACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCA    >  2:2634956/1‑100
                                                                                                    TACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCA    <  1:717865/100‑1
                                                                                                    TACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCA    <  2:744606/100‑1
                                                                                                     ACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAG   <  2:1374904/100‑1
                                                                                                     ACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAG   <  2:2814680/100‑1
                                                                                                     ACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAG   <  2:393193/100‑1
                                                                                                      CCGGATGCGTAGGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTCCA    <  1:3389131/98‑1
                                                                                                      CCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGA  >  2:1856819/1‑100
                                                                                                      CCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGA  >  2:3176269/1‑100
                                                                                                                                                                                                          
TTCAGGCGTGGCGAAGCCTTACCTGAAGGTAATAAAGCTGAGCAGGGACGTTAAAACGTCACGGCCTTACGCATCCGGTAACGCCGTGAAAAGACGGGCTTACCGGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4467583‑4467477
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttaccggaTGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGA  >  NC_000913/4467511‑4467605

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