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.

TTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTAT  >  NC_000913_3_pae_pgi/4162740‑4162981
                                                                                                   |                                                                                                                                              
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)
                            aTGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGatgcactgcgcgcgcagaatctgaatgtgcatctga                                                                                                            >  1:55704‑M1/1‑72 (MQ=255)
                            aTGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGatgcactgcgcgcgcagaatctgaatgtgcatctga                                                                                                            <  2:55704‑M1/108‑37 (MQ=255)
                            aTGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATCCCAGGGCAaccacc                                                                  >  4:112503/1‑150 (MQ=255)
                                   ccTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGAt                                                           <  3:65714/150‑1 (MQ=255)
                                                 aTTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGatgcactgcgcgcgcagaatctgaatgtgcatctgattggtcgggcaacgaaaaccaaaatcgagctggatcaggattacatcgatgaacgtctgccgg                                             >  3:50497‑M1/1‑51 (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               >  2:86159/1‑147 (MQ=255)
                                                                                  acaTTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGgcgc               <  1:86159/147‑1 (MQ=255)
                                                                                      tCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACgg        >  2:192902/1‑150 (MQ=255)
                                                                                        aGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGa      >  2:77591/1‑150 (MQ=255)
                                                                                            tCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGattat  <  2:154408/150‑1 (MQ=255)
                                                                                              tGCGCGatgcactgcgcgcgcagaatctgaatgtgcatctgat                                                                                                           <  2:18307‑M1/43‑38 (MQ=255)
                                                                                              tGCGCGatgcactgcgcgcgcagaatctgaatgtgcatctgat                                                                                                           >  1:18307‑M1/1‑6 (MQ=255)
                                                                                               gcgcgcgcAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACt                                                       >  1:102450/1‑94 (MQ=255)
                                                                                               gcgcgcgcAGTGCATCGCGCAGGGCCTCCGCTTCCAGACCCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACt                                                       <  2:102450/94‑1 (MQ=255)
                                                                                                   |                                                                                                                                              
TTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTAT  >  NC_000913_3_pae_pgi/4162740‑4162981

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.