New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3683801 =30 (0.730)5 (0.130) 5/272 4.8 15.2% coding (1817/1989 nt) yhjK cyclic‑di‑GMP phosphodiesterase
?NC_000913 3683817 = 27 (0.690)coding (1801/1989 nt) yhjK cyclic‑di‑GMP phosphodiesterase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGCAGCCCTGCGTAGCCCATGCCGAAATCATCCAGCGCCACCCGAACTCCGGCATTGCGCAGC  >  NC_000913/3683691‑3683959
                                                                                                                              |                                                                                                                                              
cgCCCTCCGCAATCCCACCGCCCGCGTTTGCCCGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCCGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCAGCAAGCCt                                                                                                                                <  1:386048/142‑1 (MQ=255)
    cTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTcaa                                                                                                                            <  2:118084/143‑1 (MQ=255)
          aaTACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCATCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCaacaaaca                                                                                                                      >  2:307998/1‑143 (MQ=255)
                        cTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCCGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTc                                                                                                        <  2:12263/143‑1 (MQ=255)
                         ttttGCCAGACAGTCATCTTGTGCTTCAGTCTCAATGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGTAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCa                                                                                                       <  1:399874/143‑1 (MQ=255)
                               cagccagTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTCTCTCCCGGCCAGCCTTCGACAACCATTTTGTCGACTTCCAGTACGt                                                                                                 >  2:335191/1‑143 (MQ=255)
                                    attggcgTTTGGCTGCAGTCTCCACGCCCTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGAtt                                                                                            <  1:268499/139‑1 (MQ=255)
                                      tCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTCGTCGACTTTCAGTACGTCGACTgg                                                                                          >  1:556782/1‑143 (MQ=255)
                                        gcgTTGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCa                                                                                        >  2:677029/1‑143 (MQ=255)
                                              tgCTTCAGTCTCCACGCCTTCGGCAATCATGTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCAACGAtt                                                                                  <  2:604950/143‑1 (MQ=255)
                                                 ttCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGCTTTTCACTACGGCGATTGGCAACCATCTCa                                                                               >  1:140211/1‑143 (MQ=255)
                                                           tcGCCTTCTGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGTAACGATTTCATATGCTGCAg                                                                     <  1:143203/142‑1 (MQ=255)
                                                            ctccagctgCAATCATGCTTGAGTTCATGCTCTGCGCCAGCATCATCATTGCTGCAATCATGTTGCTTTCTCCCGGCATGCCTTCAACGATCATTTTGTCGATTTTCAGTACGTTGATTGGTAACGATTTCATATGCTGCAGc                                                                    <  1:231171/135‑1 (MQ=255)
                                                                  cGGCAATCATGTGTAAGTGCAGGCCCTGCGCCAGCCTGACAATCCAACGGGGGGGGCGTTTGTCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTGTCAGTACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGc                                                              <  2:133606/143‑1 (MQ=255)
                                                                          aTTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGCAGCCCTGc                                                      <  1:307998/143‑1 (MQ=255)
                                                                           tttGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCAACGATTTCATCTGCTGCCGCTGCCGCAGCCCTGCg                                                     >  2:530920/1‑143 (MQ=255)
                                                                             tGTAAGTTCAGGCTCTGCGCGAGCATGATAATTGTAGCAATCATGCTGCTATCTCCCGGCAAGACTTCAACAAACATTTGGTCGATTTTCAGTACGTCGATTGGCAACGATTTAATATGCTGCAGCTGACGCAGCCCTGCGTa                                                   <  1:174551/143‑1 (MQ=255)
                                                                              gTACGTTCAGACGCTTCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTCCGTCGATTGGCAACGCTTTCATATGCTGCAGCTGACGCCGCCCTCCGCAg                                                  >  1:258749/1‑143 (MQ=255)
                                                                                aatttCAGGCTCCCCGCCTCCTTGAGAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGCAGCCCTGCGTAGcc                                                <  2:188909/140‑1 (MQ=255)
                                                                                aaGTTCAGGCTCCGCGCCAGCATGATAATTGCAGCAATCATGCTGCCATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGCAGCCCTGCGTAGcc                                                <  1:333243/143‑1 (MQ=255)
                                                                                                   gCAGGATAATTGGAGCAAGCGTGCTGCTATCTCCCGGCAAGCCTTCAACAAACAGTTTGTCGATTTTCAGGACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGCAGCCCTGCGTAGCCCATGCCGAAATCATCCAgc                             <  2:219227/143‑1 (MQ=255)
                                                                                                        attagtGCAGCAGTCTTTCTGCTGTCTCCCGGCAAGCCTTTAACGAGCATGTTGTCGATTTTCAGTACGTCGCTTGGCTACGATTTCAGCTGCTGTAGCTCACGCAGCCCTGCGTATCCCATGCCGAAATCATCCAGCGCCAc                        <  1:192605/138‑1 (MQ=255)
                                                                                                        aTAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGCAGCCCTGCGTAGCCCATGCCGAAATCATCCAACGCCAc                        >  2:593807/1‑143 (MQ=255)
                                                                                                          aaTTGCAGCAATCATGCTGCTATCTCCCGGCTAGCCTTCAACAAACATTTTGTCGATTTTCAGGACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGCAGCCCTGCGTAGCCCATGCCGAAATCATCCAGCGCCAccc                      <  2:506409/143‑1 (MQ=255)
                                                                                                           atagcagcatgattgctgcTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGCCGCCCTGCGTAGCCCATCCCGAAATCATCCAGCGCCACCCg                     >  1:608744/14‑143 (MQ=255)
                                                                                                            tagcagcatgattgctgcTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGCAGCCCTGCCTAGCCCATGCCGAACTCATCCAGCGCCCCCCgc                    >  2:702617/13‑142 (MQ=255)
                                                                                                              gcagcaATCATGCTGCTATCTCCCGGCAAGCCTTCCACACACATTTTGTCGATTTTCAGTACGTCGATTGGCAACGATTTCCTATGCTGCAGCTGACGCAGCCCTGCGTAGCCCCTGCCGAAATCATCCCCCGCCACCCGAAc                  >  2:33198/1‑143 (MQ=255)
                                                                                                                 gcaATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGCAGCCCTGCGTAGCCCATGCCGAAATCATCCAGCGCCACCCGAACTcc               >  1:607896/1‑143 (MQ=255)
                                                                                                                   actCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGCAGCCCTGCGTAGCCCATGCCGAAATCATCCAGCGCCACCCGAACTCCgg             <  2:40034/141‑1 (MQ=255)
                                                                                                                        tgctgcTTTCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTTGGTGCGCCCACCGCCAACGATTTCATATGCTGCAGCTGACGCAGCCCTGCGTAGCCCATGCCGAAATCATCCAGCGCCACCCGAACTCCGGCATTg        <  2:710876/143‑1 (MQ=255)
                                                                                                                              tATCTCCCGGCGAGCCGTCAACGATCATTTTGTCGATTTTCGGTACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGCAGCCCTGCGTAGCCCATGCCGAAATCATCCAGCGCCACCCGAACTCCGGCATTGCGCAGc  <  1:174285/143‑1 (MQ=255)
                                                                                                                              |                                                                                                                                              
AGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGCAGCCCTGCGTAGCCCATGCCGAAATCATCCAGCGCCACCCGAACTCCGGCATTGCGCAGC  >  NC_000913/3683691‑3683959

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