New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 28207133 (0.970)8 (0.260) 8/238 2.2 17.5% intergenic (‑88/‑207) yagA/yagE CP4‑6 prophage; putative DNA‑binding transcriptional regulator/2‑keto‑3‑deoxy gluconate (KDG) aldolase; CP4‑6 prophage
?NC_000913 = 282106 45 (1.450)intergenic (‑123/‑172) yagA/yagE CP4‑6 prophage; putative DNA‑binding transcriptional regulator/2‑keto‑3‑deoxy gluconate (KDG) aldolase; CP4‑6 prophage
Rejected: Frequency below/above cutoff threshold.

AGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAAC  >  NC_000913/281965‑282253
                                                                                                                                             |                                                                                                                                                   
aGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGccccttt                                                                                                                                              <  7:5222‑M2/149‑8 (MQ=255)
         gCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAgg                                                                                                                                     <  2:37016/149‑1 (MQ=255)
                cATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCgggg                                                                                                                              >  5:132296/1‑149 (MQ=255)
                      cGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTGAggg                                                                                                                        <  1:58209/149‑1 (MQ=255)
                             ttATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAatat                                                                                                                 >  2:170237/1‑149 (MQ=255)
                                 gCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAat                                                                                                                   >  1:154966/1‑143 (MQ=255)
                                 gCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAat                                                                                                                   <  2:154966/143‑1 (MQ=255)
                                  ccTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGa                                                                                                             <  4:26493/148‑1 (MQ=255)
                                                     ccATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCa                                                                                         >  4:78469/1‑149 (MQ=255)
                                                           ccccGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCct                                                                                   >  6:142835/1‑149 (MQ=255)
                                                            cccGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCaa                                                                                                                    >  8:94354/1‑115 (MQ=255)
                                                            cccGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCaa                                                                                                                    <  7:94354/115‑1 (MQ=255)
                                                            cccGAACAAGTGTTCACTATGCCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGcccctttggggtgtacggtccggggacatagtgaacacttgttcggggacatggtagacacttacaa                                                                                  >  7:57889‑M2/1‑82 (MQ=255)
                                                             ccGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGcccctttggggtgtacggtccggggacatagtgaacacttgttcggggacatggtagacacttacca                                                                                  <  8:57889‑M2/148‑68 (MQ=255)
                                                                      tgtTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGcccctt                                                                                                                                               <  2:17945‑M2/78‑7 (MQ=255)
                                                                      tgtTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGcccctt                                                                                                                                               >  1:17945‑M2/1‑72 (MQ=255)
                                                                      tgtTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGcc                                                                                                                                  <  1:139565/91‑1 (MQ=255)
                                                                      tgtTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGcc                                                                                                                                  >  2:139565/1‑91 (MQ=255)
                                                                       gtTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGCGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTCACTTTCCCGCACTCCCACTTCCCTTCCGa                                                                       >  6:91438/1‑149 (MQ=255)
                                                                       gtTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGa                                                                       >  4:131531/1‑149 (MQ=255)
                                                                         tCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACCGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATAtg                                                                                                                >  8:120903/1‑106 (MQ=255)
                                                                         tCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACCGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATAtg                                                                                                                <  7:120903/106‑1 (MQ=255)
                                                                           aCTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGcccctttg                                                                                                                                             >  8:106425‑M2/1‑67 (MQ=255)
                                                                           aCTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGcccctttg                                                                                                                                             <  7:106425‑M2/75‑9 (MQ=255)
                                                                            cTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTAc                                                                  <  6:132296/149‑1 (MQ=255)
                                                                                   cccGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTccc                                                                                            >  3:74105/1‑116 (MQ=255)
                                                                                   cccGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTccc                                                                                            <  4:74105/116‑1 (MQ=255)
                                                                                    ccGGCCCGTACACCCCACAGGGGAGCGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCtt                                                          <  5:91438/149‑1 (MQ=255)
                                                                                    ccGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCtt                                                          <  3:131531/149‑1 (MQ=255)
                                                                                        aCCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTc                                                      <  3:78469/149‑1 (MQ=255)
                                                                                        aCCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTc                                                      <  5:142835/149‑1 (MQ=255)
                                                                                          cGTACACACCAAAGGGGAGAGGGGCCTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGTGCCAGCTTAAATTTCACGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTCTTCTg                                                    >  7:5739/1‑149 (MQ=255)
                                                                                                 cccAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTCAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATa                                             >  4:140064/1‑149 (MQ=255)
                                                                                                  ccAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGcccctttggggtgtacggtccggggacatagtgaacacttgttcggggacatggtagacacttacaactaaggcataagaacccgtttatggagtcgcttatgcc                                            <  6:35869‑M2/149‑106 (MQ=255)
                                                                                                    aaaGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGcccctttggggtgtacggtccggggacatagtgaacacttgttcggggacatggt                                                                                              >  4:52704‑M2/1‑42 (MQ=255)
                                                                                                    aaaGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGcccctttggggtgtacggtccggggacatagtgaacacttgttcggggacatggt                                                                                              <  3:52704‑M2/97‑56 (MQ=255)
                                                                                                        gggAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGcccctttggggtgtacggtccggggacatagtgaacacttgttcggggacatggtagacacttacaactaaggcataagaacccgtttat                                                           >  8:100004‑M2/1‑38 (MQ=255)
                                                                                                        gggAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGcccctttggggtgtacggtccggggacatagtgaacacttgttcggggacatggtagacacttacaactaaggcataagaacccgtttat                                                           <  7:100004‑M2/128‑91 (MQ=255)
                                                                                                         ggAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATc                                     <  1:170237/149‑1 (MQ=255)
                                                                                                               gggACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTtaaa                               >  3:74861/1‑149 (MQ=255)
                                                                                                                  aCTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTtaaata                             >  8:73014/1‑148 (MQ=255)
                                                                                                                    tGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAaac                          >  8:1828/1‑149 (MQ=255)
                                                                                                                      cACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAaacga                        <  7:1828/149‑1 (MQ=255)
                                                                                                                       aCCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTcttcccttcc                                                                         <  6:93861/99‑1 (MQ=255)
                                                                                                                       aCCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTcttcccttcc                                                                         >  5:93861/1‑99 (MQ=255)
                                                                                                                              cATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATa                                             >  1:180844/1‑120 (MQ=255)
                                                                                                                              cATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATa                                             <  2:180844/120‑1 (MQ=255)
                                                                                                                               aTCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATa                                             <  6:115811/119‑1 (MQ=255)
                                                                                                                                    ttCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCta             <  3:140064/146‑1 (MQ=255)
                                                                                                                                       ccccTCGCCCCTTTGGGGAGAGGGCCGGTGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCTATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTtaa                                <  5:87053/124‑1 (MQ=255)
                                                                                                                                       ccccTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTtaa                                >  6:87053/1‑124 (MQ=255)
                                                                                                                                       ccccTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAAAGAACGCTGTTCTATAATg        >  1:50270/1‑148 (MQ=255)
                                                                                                                                        cccTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATg        >  6:52173/1‑147 (MQ=255)
                                                                                                                                        cccTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATg        <  5:52173/147‑1 (MQ=255)
                                                                                                                                         ccTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATg        >  6:81412/1‑146 (MQ=255)
                                                                                                                                         ccTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTACCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATg        >  7:62101/1‑146 (MQ=255)
                                                                                                                                          cTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGa    <  4:74861/149‑1 (MQ=255)
                                                                                                                                             gcccctttggggtgtacggtccggggacatagtgaacacttgttcggggacatggtagacacttacaactaaggcataagaacccgtttatggagtcgcttatgccctgggatgcgagagataccatgtcattacgtactgagtttgt  >  5:35869‑M2/1‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
AGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAAC  >  NC_000913/281965‑282253

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