New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_vch_pgi = 3426803NA (NA)5 (0.320)
+TTTCGATTCATCATCGTGTTGCGAAAATTTGAGAGACTCAC
3/208 2.2 17.0% intergenic (‑20/+155) rrlD/alaU 23S ribosomal RNA of rrnD operon/tRNA‑Ala
?NC_000913_3_vch_pgi 4168562 = 34 (1.550)intergenic (+115/‑79) gltT/rrlB tRNA‑Glu/23S ribosomal RNA of rrnB operon
Rejected: Frequency below/above cutoff threshold.

CTTACCGACGCTTATCGCAGATTAGCACGTCCTTCATCGCCTCTGACTGCCAGGGCATCCACCGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTAAAGAGCAAAACTTCGCAGTGAACCTTTGCAGGTACACTCTGAAGTATTTTTTATTTAATCAC  >  NC_000913_3_vch_pgi/3426699‑3426949
                                                                                                        |                                                                                                                                                  
cTTACCGACGCTTATCGCAGATTAGCACGTCCTTCATCGCCTCTGACTGCCAGGGCATCCACCGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCg                                                                                                        >  1:40097/1‑149 (MQ=11)
   aCCGACGCTTATCGCAGATTAGCACGTCCTTCATCGCCTCTGACTGCCAGGGCATCCACCGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCgttg                                                                                                     <  2:253919/149‑1 (MQ=32)
         gCTTATCGCAGATTAGCACGTCCTTCATCGCCTCTGACTGCCAGGGCATCCACCGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACATCAGAGTTGCGAAACTTTGAGAGACTCACGAACAACTTTCGTTGTTCAgt                                                                                               >  2:241977/1‑149 (MQ=17)
             aTCGCAGATTAGCACGTCCTTCATCGCCTCTGACTGCCAGGGCATCCACCGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGttt                                                                                           <  1:248015/149‑1 (MQ=32)
                gCAGATTAGCACGTCCTTCATCGCCTCTGACTGCCAGGGCATACACCGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAATTTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCaa                                                                                        >  1:207394/1‑149 (MQ=17)
                                   aTCGCCTCTGACTGCCAGGGCATCCACCGTGTACGCTTAGTCGCTTAACCTCACAACCAGAAGATGTTTCACTTCAGGGTTGCGAAAATTTGAGAGAATCCCGAACAACTTTCGTGGTTCCGTGTTACAATTTTCCGCTTCATGCCGAt                                                                     >  1:177162/1‑149 (MQ=11)
                                        ctctGACTGCCAGGGCATCCACCGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTa                                                                >  1:221743/1‑149 (MQ=32)
                                                      gCATCCACCGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTAAAGAGCAAAACTTc                                                  >  1:5228/1‑149 (MQ=37)
                                                       cATCCACCGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTAAAGAGCAAAACTTcc                                                 >  1:201419/1‑148 (MQ=37)
                                                             ccGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTAAAGAGCAAAACTTCGCAGTGa                                           >  1:370272/1‑149 (MQ=255)
                                                             ccGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGAAACCGCTAGAGAGGGTCACCAGTAACGTGCGTTGTTCAGTGATTCAATGTTCAGCTTGATCCAGATTTTTAAAGAGCAAAACTTCGCAGTGa                                           >  1:187899/1‑149 (MQ=14)
                                                               gtgtACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTAAGCTTGATCCAGATTTTTAAAGAGCAAAACTTCGCGGTGAAc                                         >  1:376284/1‑149 (MQ=255)
                                                                       gtAGTCGGTTAAGCGCAAAACGCGAAGATGTTTCATTTGAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTAAAGAGCAAAACTTCGCAGTGAACCTTTGCAg                                 <  2:363645/148‑1 (MQ=14)
                                                                           tCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTAAAGAGCAAAACTTCGCAGTGAACCTTTGCAGGTac                             >  1:39783/1‑149 (MQ=255)
                                                                                   ccTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTAAAGAGCAAAACTTCGCAGTGAACCTTTGCAGGTACACTCTGaa                     <  1:173711/149‑1 (MQ=255)
                                                                                   ccTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTAAAGAGCAAAACTTCGCAGTGAACCTTTGCAGGTACACTCTGaa                     <  1:190592/149‑1 (MQ=255)
                                                                                                   tGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTAAAGAGCAAAACTTCGCAGTGAACCTTTGCAGGTACACTCTGAAGTATTTTTTATTTAAt     <  2:165431/149‑1 (MQ=40)
                                                                                                      ttCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTAAAGAGCAAAACTTCGCAGTGAACCTTTGCAGGTACACTCTGAAGTATTTTTTATTTAATCAc  <  1:279277/149‑1 (MQ=38)
                                                                                                        |                                                                                                                                                  
CTTACCGACGCTTATCGCAGATTAGCACGTCCTTCATCGCCTCTGACTGCCAGGGCATCCACCGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTAAAGAGCAAAACTTCGCAGTGAACCTTTGCAGGTACACTCTGAAGTATTTTTTATTTAATCAC  >  NC_000913_3_vch_pgi/3426699‑3426949

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