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.

GGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGC  >  NC_000913_3_vch_pgi/4168414‑4168709
                                                                                                                                                    |                                                                                                                                                   
accttaccgacgcttatcgcagattagcacgtccttcatcgcctctgactgccagggcatccaccgtgtacgcttagtcgcttaacctcacaacccgaagatgtttctttcgattcatcatcgtgttgcgaaaatttgagagactcacg                                                                                                                                                     >  1:112652‑M2/149‑149 (MQ=255)
 ccctaccgacgcttatcgcagattagcacgtccttcatcgcctctgactgccagggcatccaccgtgtacgcttagtcgcttaacctcacaacccgaagatgtttctttcgattcatcatcgtgttgcgaaaatttgagagactcacGa                                                                                                                                                    <  1:6475‑M2/2‑1 (MQ=255)
  tAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGaa                                                                                                                                                   >  2:169712/1‑149 (MQ=38)
  tAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGCGTGCCGTTGGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGaa                                                                                                                                                   <  1:269373/149‑1 (MQ=255)
     cAGGGGGTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGaacaa                                                                                                                                                <  2:343838/149‑1 (MQ=37)
     accgacgcttatcgtagattagcacgtccttcatcgcctctgactgccagggcatccaccgtgtacgcttagtcgcttaacctcacaacccgaagatgtttctttcgattcatcatcgtgttgcgaaaatttgagagactcacGaacaa                                                                                                                                                <  2:156704‑M2/6‑1 (MQ=255)
       ggggTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACg                                                                                                                                              >  1:35707/1‑149 (MQ=37)
         ggTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACgag                                                                                                                                            <  1:232130/149‑1 (MQ=38)
          gTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTGATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACgagg                                                                                                                                           >  1:65356/1‑148 (MQ=38)
          gTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCAACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACgaga                                                                                                                                           <  2:373151/149‑1 (MQ=25)
                 tCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTc                                                                                                                                    >  1:297000/1‑149 (MQ=255)
                   cgttAGGGGACGCCACGTGTTGGATTGTGAGGGAAAGTCGCCGACCTGAATATCTCAAAACGCATCTTCGGGTGAGGTGTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACGCGGAACAACGAGAGTTGTGCGt                                                                                                                                  <  2:261186/146‑1 (MQ=14)
                    ccTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAATCACTGAACAACGAGAGTTGTTCGTg                                                                                                                                 >  1:382593/1‑149 (MQ=255)
                       aGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGt                                                                                                                              >  2:38722/1‑149 (MQ=255)
                              gCCAGTTGCTGGTTTGTGAGTGAAAGTCGCCGGCTTTAATATCTGAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACGCCGAGAGTTGTCCGTGAGTCCGTCaa                                                                                                                       >  2:149076/1‑149 (MQ=11)
                                cACTTGCTGGTTTGTGAGTGAAAGTCGCCGACTTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAAt                                                                                                                     >  2:15676/1‑149 (MQ=255)
                                      tcgcctctgactgccagggcatccaccgtgtacgcttagtcgcttaacctcacaacccgaagatgtttctttcgattcatcatcgtgttgcgaaaatttgagagactcacGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGc                                                                                                               >  1:262319‑M2/111‑149 (MQ=255)
                                                     agggcatccaccgtgtacgcttagtcgcttaacctcacaacccgaagatgtttctttcgattcatcatcgtgttgcgaaaatttgagagactcacGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAAt                                                                                                <  2:262319‑M2/54‑1 (MQ=255)
                                                      aGTCGCCGACCTTAATATCTCAAAACGCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATc                                                                                               <  2:306043/149‑1 (MQ=255)
                                                           ccGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTACTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCACATTTTCGCAAGCCGATGATGAATCaaaag                                                                                          >  1:158083/1‑149 (MQ=11)
                                                           ccGACCTTAATATCTCAAAACTCATCTTCGGGTGATGCTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaag                                                                                          >  2:373078/1‑149 (MQ=25)
                                                                         tcAAAACTCATCTTCGGGTTATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGCGTTGTTCGTGAGTCCCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGt                                                                            <  1:341812/149‑1 (MQ=12)
                                                                           aaaaCTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTtg                                                                          >  2:357441/1‑149 (MQ=32)
                                                                               cTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAg                                                                      <  1:335832/149‑1 (MQ=32)
                                                                               cTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAg                                                                      <  2:309149/149‑1 (MQ=32)
                                                                               cTCATCTTCGGGTGATGTTTGAGATATCTGCCCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGCGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAg                                                                      <  1:301566/149‑1 (MQ=12)
                                                                                     ttggggTGATGTTTGAGATATTCGCTCTTTAAAAATCCGGATAAATCAGAAAATGGAAACACTGAACAACGAGAGTTGTTCGTGTGTCTCTCAAATTTTCGCAACACGATGATGAAGCGAAAGAAACATCTTCGGGTTGTGAGGTTAAg                                                                <  1:72874/146‑1 (MQ=11)
                                                                                      tCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAAGAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTGGGGTTTTAAGGTTAAGc                                                               >  2:19300/1‑149 (MQ=255)
                                                                                        gggTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGa                                                             <  2:302042/149‑1 (MQ=32)
                                                                                                                tttAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACCCGATGATGAATCGGAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACAAGGTGGATGCCCt                                     >  1:332387/1‑149 (MQ=12)
                                                                                                                tttAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTATCAAATTTTAGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTTTAAGGTTAAGCGACTAAGCGTAGACGGGGGATGCCCt                                     >  1:332452/1‑149 (MQ=11)
                                                                                                                 tgaaaaaTCTGGATCAAGCGGAAAAGTGAAACACTGAACGACGAGGGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGGAAGAAACCTCGTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTg                                    <  2:73324/147‑1 (MQ=11)
                                                                                                                      aaTCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGt                               <  1:15676/149‑1 (MQ=32)
                                                                                                                      aaTCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGt                               <  2:204115/149‑1 (MQ=32)
                                                                                                                               cAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGa                      <  2:142428/149‑1 (MQ=32)
                                                                                                                               cAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTCTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGa                      <  1:236037/149‑1 (MQ=17)
                                                                                                                                   cTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGaa                  <  2:77358/149‑1 (MQ=32)
                                                                                                                                              aaCACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTaa       <  1:230349/149‑1 (MQ=32)
                                                                                                                                                   tGAACAACGAGAGTTGTTCGTGAGGCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGc  <  1:311101/149‑1 (MQ=17)
                                                                                                                                                    |                                                                                                                                                   
GGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGC  >  NC_000913_3_vch_pgi/4168414‑4168709

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.