New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_vch_pgi 1429011 =30 (1.350)3 (0.140) 3/268 2.9 9.5% intergenic (+27/‑38) lomR/stfR pseudogene, Rac prophage lom homolog;Phage or Prophage Related; interrupted by IS5 and N‑ter deletion/Rac prophage; putative tail fiber protein
?NC_000913_3_vch_pgi 1429033 = 29 (1.380)intergenic (+49/‑16) lomR/stfR pseudogene, Rac prophage lom homolog;Phage or Prophage Related; interrupted by IS5 and N‑ter deletion/Rac prophage; putative tail fiber protein
Rejected: Frequency below/above cutoff threshold.

GGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGG  >  NC_000913_3_vch_pgi/1428864‑1429151
                                                                                                                                                   |                                                                                                                                            
tccggattttctgaggccagcgtgttcaccaccaccgtggtgctgttacgttttgctttcagctggattgtgcagttctgtaccggttttcctgtgccgtctttcagtacacctgaaatctttacttccatattcaccccacaaaaaAg                                                                                                                                             <  2:69432‑M1/2‑1 (MQ=255)
tccggattttctgaggccagcgtgttcaccaccaccgtggtgctgttacgttttgctttcagctggattgtgcagttctgtaccggttttcctgtgccgtctttcagtacacctgaaatctttactgccatattcaccccacaaaaaAg                                                                                                                                             <  2:265697‑M1/2‑1 (MQ=255)
   ggattttctgaggccagcgtgttcaccaccaccgtggtgctgttacgttttgctttcagctggattgtgcagttctgtaccggttttcctgtgccgtctttcagtacacctgaaatctttactgccatattcaccccacaaaaaAGccc                                                                                                                                          <  1:295024‑M1/5‑1 (MQ=255)
           tGCAGTTTAACCCGATCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACATCCCGCCGGTTc                                                                                                                                  <  1:247185/149‑1 (MQ=255)
                        ttcaccaccaccgtggtgctgttacgttttgctttcagctggattgtgcagttctgtaccggttttcctgtgccgtctttcagtacacctgaaatctttactgccatattcaccccacaaaaaAGCCCGCCGGTTCAGGCGGGCttttt                                                                                                                     >  1:69432‑M1/124‑149 (MQ=255)
                           cgaATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGGGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGCTTCAGGCGGGCTTTTTtgt                                                                                                                  <  2:331487/149‑1 (MQ=255)
                            gaATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTtgtg                                                                                                                 >  2:374464/1‑149 (MQ=255)
                              aTCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTggg                                                                                                               <  1:142857/149‑1 (MQ=255)
                               tCCGTGGCCATTGATATTGCTTATGCAGGCCCCGGCAGTGGCGTCTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTATCACAGTGTTATGACAGCCCGCCGGTTCACGCGGGCTTTTTTGTgggg                                                                                                              >  1:298394/1‑149 (MQ=255)
                                    ggCCCTTGATCTCGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCCCTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAat                                                                                                         <  1:11068/149‑1 (MQ=255)
                                        aTTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGGGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGCATCTgg                                                                                                     >  2:369073/1‑149 (MQ=255)
                                          tGATATTGCTTAGGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCa                                                                                                   <  1:144491/149‑1 (MQ=255)
                                                tGCTTATGAAGGCCCTGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAg                                                                                             <  1:41811/149‑1 (MQ=255)
                                                        aaGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTCGCCAGGTAACACAGTGTTATGACCGCCCGACGGTCCAGGCGGGCTTTTTTGTGGGGTGAATATGGAAGTAACTATTTCAgc                                                                                     >  1:80588/1‑148 (MQ=255)
                                                           gCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGgtgt                                                                                  >  1:204897/1‑149 (MQ=255)
                                                               cGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGATAATATGGCAGTAAAGATTTCAGCTGTACTg                                                                              >  2:384659/1‑149 (MQ=255)
                                                                  cAGCGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGaaa                                                                           <  1:225547/149‑1 (MQ=255)
                                                                                    ggttttcctgtgccgtctttcagtacacctgaaatctttactgccatattcaccccacaaaaaAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCa                                                                                                   >  2:318263‑M1/64‑107 (MQ=255)
                                                                                    ggttttcctgtgccgtctttcagtacacctgaaatctttactgccatattcaccccacaaaaaAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCa                                                                                                   <  1:318263‑M1/44‑1 (MQ=255)
                                                                                     gACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGTAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCgg                                                        <  1:374464/149‑1 (MQ=255)
                                                                                     gACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCa                                                                                                   <  1:320288/106‑1 (MQ=255)
                                                                                     gACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCa                                                                                                   >  2:320288/1‑106 (MQ=255)
                                                                                     gACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCgg                                                        <  1:190719/149‑1 (MQ=255)
                                                                                     gACGGTTTCATCGTGGGTGGCGGTTATAAGTTCTGATTAGCCAGGTAACCCAGAGTGGTTTCAACCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCgg                                                        <  1:190169/149‑1 (MQ=255)
                                                                                      aCGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGa                                                                                                           >  2:69521/1‑97 (MQ=255)
                                                                                      aCGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGa                                                                                                           <  1:69521/97‑1 (MQ=255)
                                                                                          tttCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAg                                                   >  2:176492/1‑149 (MQ=255)
                                                                                            gtagcGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACGGAAAGACGGCACAGGAAAACCGGTACAGaa                                                 <  1:312336/145‑1 (MQ=255)
                                                                                              aTCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACGGGAAAACCGGTACAGAACt                                               >  1:268203/1‑149 (MQ=255)
                                                                                              aTCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACt                                               >  2:239596/1‑149 (MQ=255)
                                                                                               tCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTg                                              <  1:176492/149‑1 (MQ=255)
                                                                                                        tCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCa                                     <  2:204897/149‑1 (MQ=255)
                                                                                                         cGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAg                                    <  1:15679/149‑1 (MQ=255)
                                                                                                                   ttCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCaa                          >  2:295024/1‑149 (MQ=255)
                                                                                                                          tAGCCAGGTAACACAGTGGTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTaa                   <  2:309215/149‑1 (MQ=255)
                                                                                                                                           gtTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTTTTTTGTTTTTATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACggtggtgg  >  2:39214/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
GGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGG  >  NC_000913_3_vch_pgi/1428864‑1429151

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