New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_tpiA = 3944978NA (NA)8 (0.150) 7/274 6.4 33.0% noncoding (1275/2904 nt) rrlC 23S ribosomal RNA of rrnC operon
?NC_000913_3_pae_tpiA = 4038811 17 (0.300)noncoding (1293/2905 nt) rrlA 23S ribosomal RNA of rrnA operon
Rejected: Coverage evenness skew score above cutoff.

ATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCCCTAA  >  NC_000913_3_pae_tpiA/4038665‑4038871
                                                                                                                                                  |                                                            
aTGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCg                                                            >  2:294927/1‑149 (MQ=255)
      cGAAGCCGCGGCAGGGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGaa                                                      <  2:412656/149‑1 (MQ=255)
                 cAGCGACCCTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGGGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGtt                                           <  1:263689/149‑1 (MQ=255)
                     gACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCCGTGAAGGTGAGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGACCGCCCGCCGGAAGACCCAGCGTTCCTg                                       >  2:642718/1‑149 (MQ=11)
                      acacTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGCAAAGACCGCTCGCCGGAAGACCAAGGGTTCCTGt                                      >  2:361510/1‑149 (MQ=255)
                        acTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGCAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTcc                                    >  1:155137/1‑149 (MQ=25)
                            tgtgTTGTTGGGTAGGGGAGTGTTCTGTAATCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGCTTCCTGTCCAGCg                                >  1:983621/1‑149 (MQ=255)
                                    tgGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCgg                        <  1:106733/149‑1 (MQ=32)
                                     gggTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCggg                       <  2:780290/149‑1 (MQ=32)
                                         aGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAg                   <  1:272945/149‑1 (MQ=32)
                                            ggAGCGTTGTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATACAGCGGGGGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGt                <  1:450124/149‑1 (MQ=12)
                                                 gTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTAACAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTc           <  1:838111/149‑1 (MQ=17)
                                                  ttCTGTAAGCCTGTGATGGCTTTCTGTGAGCCGTGATGCGCCTATTAGAATTGCGAATCATGACATAAGGAACGATAAAGCGCGGGAAAAGCCCGCTAGCCGGAAGACCAAGGGTTCCTGTAAAACGTGACTCGGGGCAGGGTGAGTCg          >  2:975790/1‑149 (MQ=2)
                                                  ttCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGCAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCCAGGGTTCCTGTCCAAAGTTAATCGGGGCAGGGTGAGTCg          >  2:827512/1‑149 (MQ=12)
                                                    cTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCGTACATAAGTAACGATAAAGCTCGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGAc        <  1:14133/149‑1 (MQ=255)
                                                     tGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGACTGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCACGGGTTCCTGTCCAAGGTTAATCGGGGCAGGGGGGGTGGAcc       >  2:847745/1‑149 (MQ=11)
                                                          gCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGCCATAAGTAACGAAAAAGCGGGGGAAAAGCCCGCTCGCCGGAAGAACAAGGGTTCCTGTCCAACGATAATCGGGGCAAGGTAAGGCGACCCCTaa  >  2:401596/1‑149 (MQ=11)
                                                                                                                                                  |                                                            
ATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCCCTAA  >  NC_000913_3_pae_tpiA/4038665‑4038871

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