New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_tpiA 1758437 =16 (0.570)5 (0.180) 5/272 3.0 23.6% coding (289/1005 nt) ldtE murein L,D‑transpeptidase
?NC_000913_3_pae_tpiA 1758456 = 17 (0.620)coding (270/1005 nt) ldtE murein L,D‑transpeptidase
Rejected: Coverage evenness skew score above cutoff.

CACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGTTTTGGCACCGGGGCGATAGTGTTATTGGCTTCAAGGATCAACATTGCCGCAGTATCAAAACGTCGGGCAATAGCCTGAAGGTTTTTATCCCCTTCTTGCAC  >  NC_000913_3_pae_tpiA/1758314‑1758596
                                                                                                                                              |                                                                                                                                            
cacGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAAGAAc                                                                                                                                        >  1:358719/1‑149 (MQ=255)
    cgTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATAACTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGtt                                                                                                                                    <  2:310072/149‑1 (MQ=255)
     gTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTtg                                                                                                                                   <  2:277382/149‑1 (MQ=255)
         ccaccaCCGGCGTTTCCAGCCCCTGCAATCCAATACCTTTTGGATAAACCTGCACAATATTTTCGCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAACTAATCGGTTATCTCCTGACGCGGTGCATCAGGTAATAACAGTTGTGaa                                                                                                                               <  2:457916/149‑1 (MQ=255)
                            ccccTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCCGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTAc                                                                                                            >  2:259451/1‑149 (MQ=255)
                                    actcAAGACCTATTGGATAAACCTGCACAATATTTTCCGCCGGCGGATAATAATAAAGGCGCAGCGCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATACCAGTTGTGAAGGAATAGTTATCGTCGTACCAGGtttt                                                                                                    <  2:293714/146‑1 (MQ=255)
                                            ccTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGTTTTGGCACCgg                                                                                            <  1:236684/149‑1 (MQ=255)
                                              tATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGCTTTGGCACCgggg                                                                                          >  2:240459/1‑149 (MQ=255)
                                                     tAAACCTGCACAATATTTGCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGTTTTGGCACCGGGGCGATAgt                                                                                   <  1:80712/149‑1 (MQ=255)
                                                         ccTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGTTTTGGCACCGGGGCGATAGTGttt                                                                               >  2:204346/1‑148 (MQ=255)
                                                                        ctcCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGTTTTGGCACCGGGGCGATAGTGTTATTGGCTTCAAGGATc                                                                <  1:156592/149‑1 (MQ=255)
                                                                                            aGGCGCAGCGCTGCAAGGTTAACGATACTCCCCTGACGCGGTGCATCAGGTAATCACAGTTGTGAAGGAATAGTTATCGGCGTACCAGGTTTTGGCACCGGGGCGATAGTGTTATTGGCTTCAAGGATCAACATTGCCGCAGTATCaaa                                            <  1:430051/149‑1 (MQ=255)
                                                                                                                aaCGATAATCCCCTGACGCGGTGCATCAGGTAATAATAGTTGTGAAGGAATAGTTATCGTCGTACCAGGTTTTGGCACCGGGGCGATAGTGTTATTGGCTTCAAGGATCAACATTGCCGCAGTATCAAAACGTCGGGCAATAGCCTGaa                        <  1:226649/149‑1 (MQ=255)
                                                                                                                   gATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGTTTTGGCACCGGGGCGATAGTGTTATTGGCTTCAAGGATCAACATTGCCGCAGTATCAAAACGTCTGGCAATAGCCTGAAGGt                     >  2:120393/1‑149 (MQ=255)
                                                                                                                           ccTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGTTTTGGCACCGGGGCGATAGTGTTATTGGCTTCAAGGATCAACATTGCCGCAGTATCAAAACGTCGGGCAATAGCCTGAAGCTTTTTATcc             >  2:310826/1‑149 (MQ=255)
                                                                                                                                 gcgGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGTTTTGGCACCGGGGCGATAGTGTTATTGGCTTCAAGGATCAACATTGCCGCAGTATCAAAACGTCGGGCAATAGCCTGAAGGTTTTTATCCCcttct       >  1:145868/1‑149 (MQ=255)
                                                                                                                                      gCATCAGGTAATAGCAGTTGTGAAGGAATAGTTATCGTCGTACCAGGTTTTGGCACCGGGGCGATAGTGTTATTGGCTTCAAGGATCAACATTGCCGCAGTATCAAAACGTCGGGCAATAGCCTGAAGGTTTTTATCCCCTTCTTGCAc  <  1:392098/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                            
CACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGTTTTGGCACCGGGGCGATAGTGTTATTGGCTTCAAGGATCAACATTGCCGCAGTATCAAAACGTCGGGCAATAGCCTGAAGGTTTTTATCCCCTTCTTGCAC  >  NC_000913_3_pae_tpiA/1758314‑1758596

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