New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_pgi = 332818119 (1.100)4 (0.250) 4/252 1.8 17.2% intergenic (+98/+58) yhbY/greA RNA binding protein associated with pre‑50S ribosomal subunits/transcript cleavage factor
?NC_000913_3_pae_pgi = 3328184 21 (1.310)intergenic (+101/+55) yhbY/greA RNA binding protein associated with pre‑50S ribosomal subunits/transcript cleavage factor
Rejected: Frequency below/above cutoff threshold.

TTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTT  >  NC_000913_3_pae_pgi/3328041‑3328318
                                                                                                                                               |                                                                                                                                      
ttATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAgcggcc                                                                                                                                  <  1:140516‑M2/150‑7 (MQ=255)
ttATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAgcggcc                                                                                                                                  <  4:82209‑M2/150‑7 (MQ=255)
        cAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCct                                                                                                                          <  2:11284/150‑1 (MQ=255)
          aCTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCcttt                                                                                                                        >  2:173672/1‑150 (MQ=255)
           cTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGGGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGATCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCctttc                                                                                                                       <  2:6736/150‑1 (MQ=255)
                 aaCGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACa                                                                                                                 <  2:26769/150‑1 (MQ=255)
                                  ccACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGa                                                                                                <  1:47352/150‑1 (MQ=255)
                                           gATTATCCTGAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGt                                                                                       >  3:94788/1‑150 (MQ=255)
                                              tATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAAt                                                                                    <  3:82209/150‑1 (MQ=255)
                                                  cTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCGCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTGCTTTACAATACATCAACATCTTTAGTATTGGGTAATTCtt                                                                                >  1:87020/1‑150 (MQ=255)
                                                   tAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCctttct                                                                                                                      >  1:116158/1‑111 (MQ=255)
                                                   tAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCctttct                                                                                                                      <  2:116158/111‑1 (MQ=255)
                                                    aaaGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTAc                                                                              >  3:72621/1‑150 (MQ=255)
                                                     aaGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACa                                                                             <  2:72690/150‑1 (MQ=255)
                                                             cacaTTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTcc                                                                     <  1:173672/150‑1 (MQ=255)
                                                                               cGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAgcggccttttatcaacgaacagagcgtggcattttg                                                                                                    >  2:4531‑M2/1‑65 (MQ=255)
                                                                               cGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAgcggccttttatcaacgaacagagcgtggcattttg                                                                                                    <  1:4531‑M2/101‑37 (MQ=255)
                                                                                gAGGGGTTTTCCGCAGGCAGGAAAGCAAAATGCCACGCACTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAAtt                                                  >  4:23849/1‑150 (MQ=255)
                                                                                                         cAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAg                                                                                               <  3:122242/80‑1 (MQ=255)
                                                                                                         cAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAg                                                                                               >  4:122242/1‑80 (MQ=255)
                                                                                                            aaTGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAgcggccttttatcaacgaacagagcgtggcattttgctctcctgcctgcggaaaa                                                                                 <  2:162589‑M2/91‑56 (MQ=255)
                                                                                                            aaTGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAgcggccttttatcaacgaacagagcgtggcattttgctctcctgcctgcggaaaa                                                                                 >  1:162589‑M2/1‑36 (MQ=255)
                                                                                                                 cACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACAACTTGAGTATTGGGTAATTCTTACAGTTATTCCACCTTAATTAGTTCAAATACTACTTCGCCGCCCGGCGTTTTGATGAcaacaa                 >  1:62374/1‑150 (MQ=255)
                                                                                                                             gTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGGGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCAtcttctt     >  2:142941/1‑150 (MQ=255)
                                                                                                                                gATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATcttcttctt  <  3:23849/150‑1 (MQ=255)
                                                                                                                                             aTAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCAtcttcttct   >  2:169684/1‑136 (MQ=255)
                                                                                                                                             aTAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCAtcttcttct   <  1:169684/136‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                      
TTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTT  >  NC_000913_3_pae_pgi/3328041‑3328318

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.