New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_pgi 3312700 =17 (0.980)4 (0.250) 3/254 2.1 21.7% coding (78/945 nt) truB tRNA pseudouridine synthase B: tRNA pseudouridine(55) synthase and putative tmRNA pseudouridine(342) synthase
?NC_000913_3_pae_pgi 3312734 = 13 (0.800)coding (44/945 nt) truB tRNA pseudouridine synthase B: tRNA pseudouridine(55) synthase and putative tmRNA pseudouridine(342) synthase
Rejected: Coverage evenness skew score above cutoff.

TCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTCCCCGAGGCAAATCGGCAACATGCCGGTCGCCAGCGGGTCCAGCGCACCGGTATGCCCGGCACGGTTGGCGTTATATATACGTTTCACTTTTTGCAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACG  >  NC_000913_3_pae_pgi/3312560‑3312834
                                                                                                                                            |                                                                                                                                      
tCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTCCCCGAGGCAAATCGGCAACATGCCGGTCGCCAGCGGGTCCAGCGCACCGGTATGCCCGGCACGGTTGGCGTTATATATACGTTTCACTTTTTGCAGCGCATCGTTGCTGGAc                                                                                                                               >  1:48305/1‑150 (MQ=255)
tCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTCCCCGAGGCAAATCGGCAACATGCCGGTCGCCAGCGGGTCCAGCGCACCGGTATGCCCGGCACGGTTGGCGTTATATATACGTTTCACTTTTTGCAGCGCATCGTTGCTGGAc                                                                                                                               <  2:149641/150‑1 (MQ=255)
                     gAAAACTTCGTCGCTTCCCCGAGGCAAATCGGCAACATGCCGGTCGCCAGCGGGTCCAGCGCACCGGTATGCCCGGCACGGTTGGCGTTATATATACGTTTCACTTTTTGCAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGc                                                                                                          <  2:192999/150‑1 (MQ=255)
                           ttCGTCGCTTCCCCGAGGCAAATCGGCAACATGCCGGTCGCCAGCGGGTCCAGCGCACCGGTATGCCCGGCACGGTTGGCGTTATATATACGTTTCACTTTTTGCAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACaaa                                                                                                    >  3:111930/1‑150 (MQ=255)
                                                               gTCGCCAGCGGGTCCAGCGCACCGGTATGCCCGGCACGGTTGGCGTTATATATACGTTTCACTTTTTGCAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGgacga                                                                >  4:102087/1‑150 (MQ=255)
                                                                tCGCCAGCGGGTCCAGCGCACCGGTATGCCCGGCACGGTTGGCGTTATATATACGTTTCACTTTTTGCAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGgacgac                                                               <  3:102087/150‑1 (MQ=255)
                                                                 cGCCAGCGGGTCCAGCGCACCGGTATGCCCGGCACGGTTGGCGTTATATATACGTTTCACTTTTTGCAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACt                                                              >  3:36717/1‑150 (MQ=255)
                                                                  gCCAGCGGGTCCAGCGCACCGGTATGCCCGGCACGGTTGGCGTTATATATACGTTTCACTTTTTGCAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTc                                                             <  4:19799/150‑1 (MQ=255)
                                                                    cAGCGGGTCCAGCGCACCGGTATGCCCGGCACGGTTGGCGTTATATATACGTTTCACTTTTTGCAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCAt                                                           >  2:110568/1‑150 (MQ=255)
                                                                      gCGGGTCCAGCGCACCGGTATGCCCGGCACGGTTGGCGTTATATATACGTTTCACTTTTTGCAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTa                                                         <  1:109963/150‑1 (MQ=255)
                                                                           tCCAGCGCACCGGTATGCCCGGCACGGTTGGCGTTATATATACGTTTCACTTTTTGCAGCGCATCGTTGCTGGa                                                                                                                                <  1:115356/74‑1 (MQ=255)
                                                                           tCCAGCGCACCGGTATGCCCGGCACGGTTGGCGTTATATATACGTTTCACTTTTTGCAGCGCATCGTTGCTGGa                                                                                                                                >  2:115356/1‑74 (MQ=255)
                                                                                             ccGGCACGGTTGGCGTTATATATACGTTTCACTTTTTGCAGCGCATCGTTGCTGGACGTAGCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCGCCTTGCTGTCGTCCGGGt                                  <  2:48305/150‑1 (MQ=255)
                                                                                                         gCGTTATATATACGTTTCACTTTTTGCAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAACACGACGtt                      <  2:103961/150‑1 (MQ=255)
                                                                                                                    aCGTTTCACTTTTTGCAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGcaac                                                                                                       <  3:73942/58‑1 (MQ=255)
                                                                                                                    aCGTTTCACTTTTTGCAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGcaac                                                                                                       >  4:73942/1‑58 (MQ=255)
                                                                                                                             tttttGCAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACg  >  4:82257/1‑150 (MQ=255)
                                                                                                                                            |                                                                                                                                      
TCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTCCCCGAGGCAAATCGGCAACATGCCGGTCGCCAGCGGGTCCAGCGCACCGGTATGCCCGGCACGGTTGGCGTTATATATACGTTTCACTTTTTGCAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACG  >  NC_000913_3_pae_pgi/3312560‑3312834

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