New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 157238029 (0.820)5 (0.150) 5/230 3.9 15.2% intergenic (‑335/+27) gadB/pqqL glutamate decarboxylase B, PLP‑dependent/putative periplasmic M16 family zinc metalloendopeptidase
?NC_000913 = 1572395 28 (0.830)intergenic (‑350/+12) gadB/pqqL glutamate decarboxylase B, PLP‑dependent/putative periplasmic M16 family zinc metalloendopeptidase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTAAAATATCACTATTAGTACTTGATTATTATTTTGAACGCATTTATAAAATTATTACATAAAAATAGCGAATATTGCTAAAATCCCCGCCAACGATGTGTTGACGGGGCTGTTATTATTTTGGCAATAATACTCCGGTATAAGTATTTACCGGATGAGAAAGATATTGTTTAACGGCAGTGTTAACATTCTCTACCGTC  >  NC_000913/1572251‑1572491
                                                                                                                                 |                                                                                                               
cgcATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTAAAATATCACTATTAGTACTTGATTATTATTTTGAACGCATTTATAAAATTATTACATAAAAATAGCGAATATTGCTAAAATCCCCGCCAACGATGTGTTGACgg                                                                                               <  3:73778/148‑1 (MQ=255)
         tgtCTCCTGAATTATCACGTAAAAATCAGACCTTAAAATATCACTATTAGTACTTGATTATTATTTTGAACGCATTTATAAAATTATTACATAAAAATAGCGAATATTGCTAAAATCCCCGCCAACGATGTGTTGACGGGGCTGt                                                                                         <  5:78118/145‑1 (MQ=255)
         tgtCTCCTGAATTATCACGTAAAAATCAGACCTTAAAATATCACTATTAGTACTTGATTATTATTTTGAACGCATTTATAAAATTATTACATAAAAATAGCGAATATTGCTAAAATCCCCGCCAACGATGTGTTGACGGGGCTGt                                                                                         >  6:78118/1‑145 (MQ=255)
           tctcCTGAATTATCACGTAAAAATCAGACCTTAAAATATCACTATTAGTACTTGATTATTATTTTGAACGCATTTATAAAATTATTACATAAAAATAGCGAATATTGCTAAAATCCCCGCCAACGATGTGTTGACGGGGCTGTtattat                                                                                   <  5:125242/149‑1 (MQ=255)
             tcCTGAATTATCACGTAAAAATCAGACCTTAAAATATCACTATTAGTACTTGATTATTATTTTGAACGCATTTATAAAATTATTACATAAAAATAGCGAATATTGCTAAAATCCCCGCCAACGATGTGTTGACGGGGCTGTTATTAttt                                                                                 >  8:54221/1‑149 (MQ=255)
                        cACGTAAAAATCAGACCTTAAAATATCACTATTAGTACTTGATTATTATTTTGAACGCATTTATAAAATTATTACATAAAAATAGCGAATATTGCTAAAATCCCCGCCAACGATGTGTTGACGGGGCTGt                                                                                         <  5:125350/130‑1 (MQ=255)
                        cACGTAAAAATCAGACCTTAAAATATCACTATTAGTACTTGATTATTATTTTGAACGCATTTATAAAATTATTACATAAAAATAGCGAATATTGCTAAAATCCCCGCCAACGATGTGTTGACGGGGCTGt                                                                                         >  6:125350/1‑130 (MQ=255)
                         aCGTAAAAATCAGACCTTAAAATATCACTATTAGTACTTGATTATTATTTTGAACGCATTTATAAAATTATTACATAAAAATAGCGAATATTGCTAAAATCCCCGCCAACGATGTGTTGACGGGGCTGTTATTATTTTGGCAATAATAc                                                                     >  2:118784/1‑149 (MQ=255)
                           gTAAAAATCAGACCTTAAAATATCACTATTAGTACTTGATTATTATTTTGAACGCATTTATAAAATTATTACATAAAAATAGCGAATATTGCTAAAATCCCCGCCAACGATGTGTTGAAGGGGCTGTTATTATTTTGGCAATAATACTc                                                                   >  7:73017/1‑149 (MQ=255)
                                    aGACCTTAAAATATCACTATTAGTACTTGATTATTATTTTGAACGCATTTATAAAATTATTACATAAAAATAGCGAATATTGCTAAAATCCCCGCCAACGATGTGTTGACGGGGCTGTTATTATTTTg                                                                               <  2:78816/128‑1 (MQ=255)
                                    aGACCTTAAAATATCACTATTAGTACTTGATTATTATTTTGAACGCATTTATAAAATTATTACATAAAAATAGCGAATATTGCTAAAATCCCCGCCAACGATGTGTTGACGGGGCTGTTATTATTTTg                                                                               >  1:78816/1‑128 (MQ=255)
                                         ttAAAATATCACTATTAGTACTTGATTATTATTTTGAACGCATTTATAAAATTATTACATAAAAATAGCGAATATTGCTAAAATCCCCGCCAACGATGTGTTGACGGGGCTGTTATTATTTTGGCAATAATACTCCGGTATAAGTAt                                                       >  2:38457/1‑147 (MQ=255)
                                         ttAAAATATCACTATTAGTACTTGATTATTATTTTGAACGCATTTATAAAATTATTACATAAAAATAGCGAATATTGCTAAAATCCCCGCCAACGATGTGTTGACGGGGCTGTTATTATTTTGGCAATAATACTCCGGTATAAGTAt                                                       <  1:38457/147‑1 (MQ=255)
                                                                gATTATTATTTTGAACGCATTTATAAAATTATTACATAAAAATAGCGAATATTGCTAAAATCCCCGCCAACGATGTGTTGACGGGGCTGTTATTATTTTGGCAATAATACTCCGGTATAAGTATTTACCGGATGAGAAAGATATTGtt                               <  1:118784/148‑1 (MQ=255)
                                                                      tattTTGAACGCATTTATAAAATTATTACATAAAAATAGCGAATATTGCTAAAATCCCCGCCAACGATGTGTTGACGGGGCTGTTATTATTTTGGCAATAATACTCCGGTATAAGTATTTACCGGATGAGAAAGATATTGt                                >  1:93795/1‑141 (MQ=255)
                                                                      tattTTGAACGCATTTATAAAATTATTACATAAAAATAGCGAATATTGCTAAAATCCCCGCCAACGATGTGTTGACGGGGCTGTTATTATTTTGGCAATAATACTCCGGTATAAGTATTTACCGGATGAGAAAGATATTGt                                <  2:93795/141‑1 (MQ=255)
                                                                       attTTGAACGCATTTATAAAATTATTACATAAAAATAGCGAATATTGCTAAAATCCCCGCCAACGATGTGTTGACGGGGCTGTTATTATTTTGGCAATAATACTCCGGTATAAg                                                          >  6:138693/1‑114 (MQ=255)
                                                                       attTTGAACGCATTTATAAAATTATTACATAAAAATAGCGAATATTGCTAAAATCCCCGCCAACGATGTGTTGACGGGGCTGTTATTATTTTGGCAATAATACTCCGGTATAAg                                                          <  5:138693/114‑1 (MQ=255)
                                                                       attTTGAACGCATTTATAAAATTATTACATAAAAATAGCGAATATTGCTAAAATCCCCGCCAACGATGTGTTGACGGGGCTGTTATTATTTTGGCAATAATACTCCGGTATAAGTATTTACCGGATGAGAAAGATATTGTTTAACGGc                        <  8:73017/148‑1 (MQ=255)
                                                                        ttttGAACGCATTTATAAAATTATTACATAAAAATAGCGAATATTGCTAAAATCCCCGCCAACGATGTGTTGCCGGGGCTGTTATTATTTTGGCAATAATACTCCGGTATAAGTATTTACCGGATGAGAAAGATATTGTTTAACGGCAg                      >  6:124673/1‑149 (MQ=255)
                                                                                  aTTTATAAAATTATTACATAAAAATAGCGAATATTGCTAAAATCCCCGCCAACGATGTGTTGACGGGGCTGTTATTATTTTGGCAATAATACTCCGGTATAAGTATTTACCgg                                                <  4:31523/113‑1 (MQ=255)
                                                                                  aTTTATAAAATTATTACATAAAAATAGCGAATATTGCTAAAATCCCCGCCAACGATGTGTTGACGGGGCTGTTATTATTTTGGCAATAATACTCCGGTATAAGTATTTACCgg                                                >  3:31523/1‑113 (MQ=255)
                                                                                  aTTTATAAAATTATTACATAAAAATAGCGAATATTGCTAAAATCCCCGCCAACGATGTGTTGACGGGGCTGTTATTATTTTGGCAATAATACTCCGGTATAAGTATTTACCGGATGAGAAAGATATTGTTTAACGGCAGTGTTAACAtt            >  6:27328/1‑149 (MQ=255)
                                                                                      ataAAATTATTACATAAAAATAGCGAATATTGCTAAAATCCCCGCCAACGATGTGTTGACGGGGCTGTTATTATTTTGGCAATAATACTCCGGTATAAGTATTTACCGGATGAGAAAGATATTGTTTAACGGCAGTGTTAACATtctc         >  6:137841/1‑148 (MQ=255)
                                                                                      ataAAATTATTACATAAAAATAGCGAATATTGCTAAAATCCCCGCCAACGATGTGTTGACGGGGCTGTTATTATTTTGGCAATAATACTCCGGTATAAGTATTTACCGGATGAGAAAGATATTGTTTAACGGCAGTGTTAACATtctc         <  5:137841/148‑1 (MQ=255)
                                                                                         aaaTTATTACATAAAAATAGCGAATATTGCTAAAATCCCCGCCAACGATGTGTTGACGGGGCTGTTATTATTTTGGCAATAATACTCCGGTATAAGTATTTACCGGATGAGAAAGATATTGTTTAACGGCAGTGTTAAc               <  5:63764/139‑1 (MQ=255)
                                                                                         aaaTTATTACATAAAAATAGCGAATATTGCTAAAATCCCCGCCAACGATGTGTTGACGGGGCTGTTATTATTTTGGCAATAATACTCCGGTATAAGTATTTACCGGATGAGAAAGATATTGTTTAACGGCAGTGTTAAc               >  6:63764/1‑139 (MQ=255)
                                                                                         aaaTTATTACATAAAAATAGCGAATATTGCTAAAATCCCCGCCAACGATGTGTTGACGGGGCTGTTATTATTTTGGCAATAATACTCCGGTATAAGTATTTACCGGATGAGAAAGATATTGTTTAACGGCAGTGTTAACATTCTCTAcc     <  7:54221/149‑1 (MQ=255)
                                                                                            ttattaCATAAAAATAGCGAATATTGCTAAAATCCCCGCCAACGATGTGTTGACGGGGCTGTTATTATTTTGGCAATAATACTCCGGTATAAGTATTTACCGGATGAGAAAGATATTGTTTAACGGCAGTGTTAACATTCTCTACCGTc  <  5:27328/149‑1 (MQ=255)
                                                                                                                                 |                                                                                                               
CGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTAAAATATCACTATTAGTACTTGATTATTATTTTGAACGCATTTATAAAATTATTACATAAAAATAGCGAATATTGCTAAAATCCCCGCCAACGATGTGTTGACGGGGCTGTTATTATTTTGGCAATAATACTCCGGTATAAGTATTTACCGGATGAGAAAGATATTGTTTAACGGCAGTGTTAACATTCTCTACCGTC  >  NC_000913/1572251‑1572491

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