New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 591569 =20 (1.780)3 (0.280) 3/250 NT 13.9% coding (1610/2238 nt) nfrB bacteriophage N4 receptor, inner membrane subunit
?NC_000913 591595 = 18 (1.680)coding (1584/2238 nt) nfrB bacteriophage N4 receptor, inner membrane subunit

TGCCGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTCGGCAATCAGCGAGGAAGGGATCTGCCAGGCATCGATGGATTCCCACGCCACGCCGTTTTGCTCTGCCAGCGCCTGTGCCAGCTGCTCGGCGCTAATCAGCCCCTGCATCAGCATTGAACCGCCCAGGCGTAGACCTTCGACGCGATTACGCAGTGCTGTATCGAGTTGTTCTTCAGTGATGACCTGATTTTCCAGCAGAATTTGACCTAAC  >  NC_000913/591426‑591712
                                                                                                                                               |                                                                                                                                               
tGCCGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTCGGCAATCAGCGAGGAAGGGATCTGCCAGGCATCGATGGATTCCCACGCCACGCCGTTTTGCTCTGCCAGCGcc                                                                                                                                            >  3:98647/1‑149 (MQ=255)
        aTTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTCGGCAATCAGCGAGGAAGGGATCTGCCAGGCATCGATGGATTCCCACGCCACGCCGTTTTGCTCTGCCAGCGCCTGTGCCAg                                                                                                                                    >  2:103129/1‑149 (MQ=255)
             cTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTCGGCAATCAGCGAGGAAGGGATCTGCCAGGCATCGATGGATTCCCACGCCACGCCGTTTTGCTCTGCCAGCGCCTGTGCCAgctgct                                                                                                                               <  3:42311/149‑1 (MQ=255)
                     gcccgttagttcaaattctgctggaaaatcaggtcatcactgaagaacaactcgatacagcactgcgtaatcgcgtcgaaggtctacgcctgggcggttcaatgctgatgcaggggctgattAGCGCCTGTGCCAGCTGCTCGGCGCTa                                                                                                                       >  2:32287‑M1/123‑149 (MQ=255)
                           gACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTCGGCAATCAGCGAGGAAGGGATCTGCCAGGCATCGATGGATTCCCACGCCACGCCGTTTTGCTCTGCCAGCGCCTGTGCCAGCTGCTCGGCGCTAATCAGc                                                                                                                 <  2:83248/149‑1 (MQ=255)
                                aGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTCGGCAATCAGCGAGGAAGGGATCTGCCAGGCATCGATGGATTCCCACGCCACGCCGTTTTGCTCTGCCAGCGCCTGTGCCAGCTGCTCGGCGCTAATCAGCCCCTg                                                                                                            >  4:21071/1‑149 (MQ=255)
                                              aTAATGCAGCGCCACGGAGGCCGGCATTTCGGCAATCAGCGAGGAAGGGATCTGCCAGGCATCGATGGATTCCCACGCCACGCCGTTTTGCTCTGCCAGCGCCt                                                                                                                                           >  1:10881/1‑104 (MQ=255)
                                              aTAATGCAGCGCCACGGAGGCCGGCATTTCGGCAATCAGCGAGGAAGGGATCTGCCAGGCATCGATGGATTCCCACGCCACGCCGTTTTGCTCTGCCAGCGCCt                                                                                                                                           <  2:10881/104‑1 (MQ=255)
                                                   gCAGCGCCACGGAGGCCGGCATTTCGGCAATCAGCGAGGAAGGGATCTGCCAGGCATCGATGGATTCCCACGCCACGCCGTTTTGCTCTGCCAGCGCCTGTGCCAGCTGCTCGGCGCTAATCAGCCCCTGCATCAGCATTGAACCGccc                                                                                         >  2:37971/1‑149 (MQ=255)
                                                           aCGGAGGCCGGCATTTCGGCAATCAGCGAGGAAGGGATCTGCCAGGCATCGATGGATTCCCACGCCACGCCGTTTTGCTCTGCCAGCGCCTGTGCCAGCTGCTCGGCGCTAATCAGCCCCTGCATCAGCATTGAACCGCCCAGGCGTAg                                                                                 >  1:61236/1‑149 (MQ=255)
                                                           aCGGAGGCCGGCATTTCGGCAATCAGCGAGGAAGGGATCTGCCAGGCATCGATGGATTCCCACGCCACGCCGTTTTGCTCTGCCAGCGCCTGTGCCAGCTGCTCGGCGCTAATCAGCCCCTGCATCAGCATTGAACCGCCCAGGCGTAg                                                                                 >  2:16237/1‑149 (MQ=255)
                                                            cGGAGGCCGGCATTTCGGCAATCAGCGAGGAAGGGATCTGCCAGGCATCGATGGATTCCCACGCCACGCCGTTTTGCTCTGCCAGCGCCTGTGCCAGCTGCTCGGCGCTAATCAGCCCCTGCATCAGCATTGAACCGCCCAGGCGTAGa                                                                                <  4:98647/149‑1 (MQ=255)
                                                                    ggCATTTCGGCAATCAGCGAGGAAGGGATCTGCCAGGCATCGATGGATTCCCACGCCACGCCGTTTTGCTCTGCCAGCGCCTGTGCCAGCTGCTCGGCGCTAATCAGCCCCTGCATCAGCATTGAACCGCCCAGGCGTAGACCTTCGAc                                                                        >  3:44279/1‑149 (MQ=255)
                                                                                       aGGAAGGGATCTGCCAGGCATCGATGGATTCCCACGCCACGCCGTTTTGCTCTGCCAGCGCCTGTGCCAGCTGCTCGGCGCTAATCAGCCCCTGCATCAGCATTGAACCGCCCAGGCGTAGACCTTCGACGCGATTACGCAGTGCTGTa                                                     >  1:84498/1‑149 (MQ=255)
                                                                                        ggAAGGGATCTGCCAGGCATCGATGGATTCCCACGCCACGCCGTTTTGCTCTGCCAGCGCCTGTGCCAGCTGCTCGGCGCTAATCAGCCCCTGCATCAGCATTGAACCGCCCAGGCGTAGACCTTCGACGCGATTACGCAGTGCTGTAt                                                    <  1:16237/149‑1 (MQ=255)
                                                                                                              aTGGATTCCCACGCCACGCCGTTTTGCTCTGCCAGCGCCTGTGCCAGCTGCTCGGCGCTAATCAGCCCCTGCATCAGCATTGAACCGCCCAGGCGTAGACCTTCGACGCGATTACGCAGTGCTGTATCGAGTTGTTCTTCAGTGATGAc                              <  2:84498/149‑1 (MQ=255)
                                                                                                                  aTTCCCACGCCACGCCGTTTTGCTCTGCCAGCGCCTGTGCCAGCTGCTCGGCGCTAATCAGCCCCTGCATCAGCATTGAACCGCCCAGGCGTAGACCTTCGACGCGATTACGCAGTGCTGTATCGAGTTGTTCTTCAGTGATGACCTGa                          <  2:103447/149‑1 (MQ=255)
                                                                                                                      ttcaatgctgatgcaggggctgattAGCGCCTGTGCCAGCTGCTCGGCGCTAATCa                                                                                                                   <  3:92316‑M1/31‑1 (MQ=255)
                                                                                                                      ttcaatgctgatgcaggggctgattAGCGCCTGTGCCAGCTGCTCGGCGCTAATCa                                                                                                                   >  4:92316‑M1/26‑56 (MQ=255)
                                                                                                                      ccacgccacgCCGTTTTGCTCTGCCAGCGCCTGTGCCAGCTGCTCGGCGCTAATCAGCCCCTGCATCAGCATTGAACCGCCCAGGCGTAGACCTTCGACGCGATTACGCAGTGCTGTATCGAGTTGTTCTTCAGTGATGACCTGAtttt                      >  1:51601/1‑149 (MQ=255)
                                                                                                                       cacgccacgcCGTTTTGCTCTGCCAGCGCCTGTGCCAGCTGCTCGGCGCTAATCAGCCCCTGCATCAGCATTGAACCGCCCAGGCGTAGACCTTCGACGCGATTACGCAGTGCTGTATCGAGTTGTTCTTCAGTGATGACCTGATTTTc                     <  4:22058/149‑1 (MQ=255)
                                                                                                                                    tttGCTCTGCCAGCGCCTGTGCCAGCTGCTCGGCGCTAATCAGCCCCTGCATCAGCATTGAACCGCCCAGGCGTAGACCTTCGACGCGATTACGCAGTGCTGTATCGAGTTGTTCTTCAGTGATGACCTGATTTTCCAGCAGAATTTGa        <  2:56925/149‑1 (MQ=255)
                                                                                                                                     ttGCTCTGCCAGCGCCTGTGCCAGCTGCTCGGCGCTAATCAGCCCCTGCATCAGCATTGAACCGCCCAGGCGTAGACCTTCGACGCGATTACGCAGTGCTGTATCGAGTTGTTCTTCAGTGATGACCTGATTTTCCAGCAGAATTTGAc       <  3:21071/149‑1 (MQ=255)
                                                                                                                                          ctGCCAGCGCCTGTGCCAGCTGCTCGGCGCTAATCAGCCCCTGCATCAGCATTGAACCGCCCAGGCGTAGACCTTCGACGCGATTACGCAGTGCTGTATCGAGTTGTTCTTCAGTGATGACCTGATTTTCCAGCAGAATTTGACCTAAc  <  1:87003/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                               
TGCCGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTCGGCAATCAGCGAGGAAGGGATCTGCCAGGCATCGATGGATTCCCACGCCACGCCGTTTTGCTCTGCCAGCGCCTGTGCCAGCTGCTCGGCGCTAATCAGCCCCTGCATCAGCATTGAACCGCCCAGGCGTAGACCTTCGACGCGATTACGCAGTGCTGTATCGAGTTGTTCTTCAGTGATGACCTGATTTTCCAGCAGAATTTGACCTAAC  >  NC_000913/591426‑591712

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