New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_tpiA 2729630 =NA (NA)11 (0.290) 7/274 4.0 NA noncoding (1528/1542 nt) rrsG 16S ribosomal RNA of rrnG operon
?NC_000913_3_pae_tpiA 2729652 = NA (NA)noncoding (1506/1542 nt) rrsG 16S ribosomal RNA of rrnG operon
Rejected: Coverage evenness skew score above cutoff.

AACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTG‑CTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAA  >  NC_000913_3_pae_tpiA/2729512‑2729714
                                                                                                                                             |                                                              
aaCGGCCTTATTCTCTTCAGCCTCACTCCCAAAGCGTCACAGCCATGTCCACT‑ATGATTCGTCAGTGAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTa                                                        <  2:665482/149‑1 (MQ=11)
 aCGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTG‑CTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTAc                                                       <  1:352429/149‑1 (MQ=40)
       ttATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTG‑CTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGGTCCCCTACGGTTACCTTGTTACGACTTc                                                 >  2:371796/1‑149 (MQ=255)
               ttCAGCCTCACTCCCAACGCGTAAACGCCTTG‑CTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGt                                         >  1:295984/1‑149 (MQ=255)
                  aGCCTCACTCCCAACGCGTAAACGCCTTG‑CTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCAt                                      <  1:544938/149‑1 (MQ=40)
                  aGCCTCACTCCCAACGCGTAAACGCCTTG‑CTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTAAt                                      >  1:122487/1‑149 (MQ=255)
                      tCACTCCCAACGCGTAAACGCCTTG‑CTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAAt                                  >  2:519700/1‑149 (MQ=255)
                       cACTCCCAACGCGTAAACGACTTG‑CTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCCGGTTCCCCTACGGTTACCTTGTTACGTCTTCACCCCAGGCAAGacgt                                 >  2:83467/1‑146 (MQ=14)
                            ccAACGCGTAAACGCCTTG‑CTATTCGCTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACaaa                            <  1:535962/149‑1 (MQ=255)
                            ccAACGCGTAAACGCCTTG‑CTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACaaa                            <  1:519700/149‑1 (MQ=40)
                                 gcgTAAACGCCTTG‑CTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGt                       >  1:377952/1‑149 (MQ=39)
                                    tAAACGCCTTG‑CTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAg                    >  2:271161/1‑149 (MQ=255)
                                           cTTG‑CTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGCTAGGCGCCCTc             <  1:426990/149‑1 (MQ=255)
                                                 tATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGaa        <  1:632761/149‑1 (MQ=35)
                                                    tCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGCAGGt     >  2:312496/1‑149 (MQ=17)
                                                     cACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCGTCTTTAATGTAAGGAGGGGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGtt    <  2:247853/149‑1 (MQ=12)
                                                       cTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTaa  <  2:553267/149‑1 (MQ=11)
                                                                                                                                             |                                                              
AACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTG‑CTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAA  >  NC_000913_3_pae_tpiA/2729512‑2729714

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.