New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_pgi = 416283918 (1.040)3 (0.180) 3/258 2.1 14.3% coding (432/1101 nt) trmA tRNA m(5)U54 methyltransferase, SAM‑dependent; tmRNA m(5)U341 methyltransferase
?NC_000913_3_pae_pgi = 4162849 19 (1.160)coding (422/1101 nt) trmA tRNA m(5)U54 methyltransferase, SAM‑dependent; tmRNA m(5)U341 methyltransferase
Rejected: Frequency below/above cutoff threshold.

TTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTATTACGCACACCC  >  NC_000913_3_pae_pgi/4162740‑4162992
                                                                                                             |                                                                                                                                               
tttCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAg                                                                                                         >  4:65714/1‑150 (MQ=255)
   cccTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGctt                                                                                                      <  2:62797/150‑1 (MQ=255)
           ccGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATcgcgcagaatctgaatgtgcatctgattggtcgggcaacgaaaaccaaaat                                                                                              >  3:50497‑M2/1‑99 (MQ=255)
              gCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGtata                                                                                           <  3:64236/150‑1 (MQ=255)
                      tCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGaa                                                                                   <  2:18491/150‑1 (MQ=255)
                        aTCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAac                                                                                 >  1:154408/1‑150 (MQ=255)
                            aTGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATCCCAGGGCAaccacc                                                                             >  4:112503/1‑150 (MQ=255)
                                   ccTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGAt                                                                      <  3:65714/150‑1 (MQ=255)
                                                              ttGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAgg                                                                                      >  1:103149/1‑107 (MQ=255)
                                                              ttGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAgg                                                                                      <  2:103149/107‑1 (MQ=255)
                                                                gCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTgg                                         >  2:148253/1‑150 (MQ=255)
                                                                  ccGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGaa                                       >  1:53752/1‑150 (MQ=255)
                                                                         aTCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATcgcgca                                                                                                                                           <  2:18307‑M2/43‑7 (MQ=255)
                                                                         aTCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATcgcgca                                                                                                                                           >  1:18307‑M2/1‑37 (MQ=255)
                                                                          tCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATcgcgcagaatctgaatgtgcatctgattggtcgggcaacgaaaaccaaaatcgagctggatcaggattacat                                                                         <  2:55704‑M2/108‑73 (MQ=255)
                                                                          tCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATcgcgcagaatctgaatgtgcatctgattggtcgggcaacgaaaaccaaaatcgagctggatcaggattacat                                                                         >  1:55704‑M2/1‑36 (MQ=255)
                                                                                  acaTTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGgcgc                          <  1:86159/147‑1 (MQ=255)
                                                                                  acaTTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGgcgc                          >  2:86159/1‑147 (MQ=255)
                                                                                      tCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACgg                   >  2:192902/1‑150 (MQ=255)
                                                                                        aGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGa                 >  2:77591/1‑150 (MQ=255)
                                                                                            tCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGattat             <  2:154408/150‑1 (MQ=255)
                                                                                               gcgcgcgcAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACt                                                                  >  1:102450/1‑94 (MQ=255)
                                                                                               gcgcgcgcAGTGCATCGCGCAGGGCCTCCGCTTCCAGACCCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACt                                                                  <  2:102450/94‑1 (MQ=255)
                                                                                                       aGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTATTACGCACAccc  >  1:45869/1‑150 (MQ=255)
                                                                                                             |                                                                                                                                               
TTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTATTACGCACACCC  >  NC_000913_3_pae_pgi/4162740‑4162992

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