New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =2 (0.050)38 (1.010) 31/248 0.2 92.9% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 4 (0.110)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

ATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAA  >  NC_000913/1209482‑1209752
                                                                                                                                         |                                                                                                                                     
tgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTa                                                                                                                                  <  2:224008‑M2/6‑1 (MQ=255)
         ggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTCATGAAAATga                                                                                                                         >  2:529356‑M2/129‑143 (MQ=255)
              tccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccgtctcccaaaccaaCGTTTATGAAAATGAAGaaa                                                                                                                    <  1:719126‑M2/20‑1 (MQ=255)
                 ccccccactgaatgattactgattgctcaatgccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATaa                                                                                                                 <  2:533790‑M2/23‑1 (MQ=255)
                   gggTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACa                                                                                                               <  2:701710/143‑1 (MQ=255)
                      ccgcgaatgattactgattgctcaacaccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGc                                                                                                            <  1:601447‑M2/28‑1 (MQ=255)
                         tgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatatagccgttgatgaattgagttttgctccttctcccaaaccaaAGTTTATGAAATTGAAGAAATAACAAGCaaa                                                                                                         >  2:196600‑M2/113‑143 (MQ=255)
                           aatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcatcgagttttcctccttctcccaaaccaaAGTTTATGAAAATGAAGAAATAACAAGCAAATg                                                                                                       <  2:408670‑M2/33‑1 (MQ=255)
                           aatgatgactggttgctcaatacccacaatcagagccggtatagcgatgaacccggtgcttcctaatgtagccgttgctgcattgagttttgcaccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATg                                                                                                       <  2:243065‑M2/33‑1 (MQ=255)
                              gattactgattgctcaataccaccaatcgtaaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTCATGAACATGAAGAACTAACAAGCAACTCGca                                                                                                    >  2:642748‑M2/108‑143 (MQ=255)
                               agtactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat                                                                                                   <  1:233315‑M2/37‑1 (MQ=255)
                                 tactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                                                                                                 >  1:724380‑M2/105‑143 (MQ=255)
                                      attgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                                                                                            <  1:373783‑M2/44‑1 (MQ=255)
                                       ttgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAGATAACAAGCAAATGGCATCATTCCTg                                                                                           <  1:760182‑M2/45‑1 (MQ=255)
                                       ctgcccaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg                                                                                           <  1:298614‑M2/45‑1 (MQ=255)
                                        tgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAGAATGAAGAAATAACAAGCAAATGGCATCATTCCtgg                                                                                          >  1:636260‑M2/98‑142 (MQ=255)
                                              ataccaccaatcatatccgggatatctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAc                                                                                    <  2:39317‑M2/52‑1 (MQ=255)
                                                         cctcaccgttatagctatgaaaccggtgcgtcctaatgtagccgttgccgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAATTAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt                                                                         <  1:67581‑M2/63‑1 (MQ=255)
                                                         cataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATACCAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt                                                                         >  1:165664‑M2/81‑143 (MQ=255)
                                                          ataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCCTtaaccaggggcgataa                                                                        >  2:206886‑M2/80‑137 (MQ=255)
                                                           gcccgccGGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCCACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa                                                                       <  1:34901/143‑1 (MQ=255)
                                                               cggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctcattctcccaaacgaaCGTTTATGACAATGAAGAAATAACAAGCAAATGGCACCATTCGTGCTTGTACCAGGGGGatgtaaactg                                                                   >  2:663907‑M2/75‑139 (MQ=255)
                                                               cggtatagctatgaaaccggcgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATg                                                                   >  1:414551‑M2/75‑143 (MQ=255)
                                                                   ttagctatgataccggtggtccctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAGCAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTa                                                               <  2:505511‑M2/73‑1 (MQ=255)
                                                                   atagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGATTa                                                               <  2:318255‑M2/73‑1 (MQ=255)
                                                                     agctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGCAGAAATAACAAGCAAATGGCATCATTCCTACTTTTACCAGTGGGATTTAACACGcttctt                                                             >  1:654778‑M2/69‑143 (MQ=255)
                                                                                 cgtgcgttctaatgtggcccttgctgcattgagtttttctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAATAAGCAAATGGCATCATTCCTGCTTTTACAAGCGAAAATTAACATGCTTATTGGCTATGTAcgc                                                 <  2:475724‑M2/87‑1 (MQ=255)
                                                                                     ggtcACGGAGGCGGGACACCAAGAGGTAGTGCAGAGCTTTGTCCGAAGCGAACGTGTATGAAAGTGAAGAAATAGCAAGCAAATGGCATCATTCCTGCGTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTAt                                             <  2:19351/140‑1 (MQ=255)
                                                                                       tcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATca                                           >  1:523446‑M2/51‑143 (MQ=255)
                                                                                       ccACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTACCAAACCAACGTATATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCAGTTACCAGGGGGCTGTAACATGCCTATTGGCTATGCACGCGTAGca                                           >  2:742848/1‑143 (MQ=255)
                                                                                           aatgtagccgttgcttcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaacaa                                       <  1:763440‑M2/97‑1 (MQ=255)
                                                                                             tgtagccgttgctgcattgagttttgctccttctcccaatccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAAt                                     <  1:93936‑M2/99‑1 (MQ=255)
                                                                                             tgtagccgttgctgcattgaggtttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAAt                                     <  2:231330‑M2/99‑1 (MQ=255)
                                                                                                agccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGAc                                  <  2:500591‑M2/102‑1 (MQ=255)
                                                                                                agccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGAc                                  <  1:441619‑M2/102‑1 (MQ=255)
                                                                                                  ctgttgcttcattgagtttttctccttcccccaaaccagCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCa                                <  2:695523‑M2/104‑1 (MQ=255)
                                                                                                     ttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCCTCATTCCTGCTTTTACCAGGGGGATTTAAAATGGCTACATGCGATGTACGCGTATCAACAAATGACCAGaa                             >  2:341776‑M2/37‑143 (MQ=255)
                                                                                                       gctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAGTGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCGATGTACGCGTATCAACACATGACCAGAaca                           >  2:75750‑M2/35‑143 (MQ=255)
                                                                                                                  ttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACaa                >  2:774707‑M2/24‑143 (MQ=255)
                                                                                                                   tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAAc               <  1:774707‑M2/121‑1 (MQ=255)
                                                                                                                   tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAAc               <  2:356172‑M2/121‑1 (MQ=255)
                                                                                                                    ttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTAGGTACGCGTATCAACAAATGACCAGAACACAGAACTACAACg              >  1:407363‑M2/22‑143 (MQ=255)
                                                                                                                        tccTTCTCCCAAACCAACGTTTCTGCAAATGAAGAAATAACACGCAACTGCCATGATGCGTGCCCTCACGCCGATTTTTTAACCTGCTTATTGGCTCGGTACGCGTGTGAACAAATGACAAGAACACAGATCTACCACGTAAt          >  2:681405/2‑143 (MQ=255)
                                                                                                                               cccaaaccaaCGATTAAGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATATACAACGTAATGCGCAGa   >  2:765122‑M2/11‑143 (MQ=255)
                                                                                                                                ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGaa  >  2:724138‑M2/10‑143 (MQ=255)
                                                                                                                                         |                                                                                                                                     
ATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAA  >  NC_000913/1209482‑1209752

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