New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_vch_tpiA 1207790 =13 (0.440)19 (0.730) 18/254 0.5 71.0% noncoding (11571/15193 nt) cryptic prophage e14 cryptic prophage e14
?NC_000913_3_vch_tpiA 1209619 = 4 (0.150)noncoding (13400/15193 nt) cryptic prophage e14 cryptic prophage e14

GTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATG  >  NC_000913_3_vch_tpiA/1209477‑1209764
                                                                                                                                              |                                                                                                                                                 
gTCGTATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGAAGAGCCTTCTCCCAAACCAACGCTTAt                                                                                                                                             >  2:90363/1‑149 (MQ=255)
   gatgcttttgcggcaatccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaa                                                                                                                                          <  2:93468‑M2/10‑1 (MQ=255)
        ttttgcggcactccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcctcctaatgtagccgttgctgcagtgagttttgctccttctcccaaaccaaCGTTTATGAAAATga                                                                                                                                     <  2:159812‑M2/15‑1 (MQ=255)
         tttgcgggattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaa                                                                                                                                    <  2:498890‑M2/16‑1 (MQ=255)
          ttgcggcattccacccccacggaatgattactgatggctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaag                                                                                                                                   <  1:51184‑M2/17‑1 (MQ=255)
               gcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaa                                                                                                                                         >  2:484119‑M2/128‑138 (MQ=255)
               gcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaa                                                                                                                                         <  1:484119‑M2/11‑1 (MQ=255)
                    ccacccccactgaatgattactgatcgctcaataccaccaatcataaccggtatagctatgatgccgctgcgtcctaatgtagccgctgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAg                                                                                                                         <  1:505594‑M2/27‑1 (MQ=255)
                            actgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                                                                                                                 >  1:213079‑M2/115‑149 (MQ=255)
                                                tcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAg                                                                                             >  2:103289‑M2/95‑149 (MQ=255)
                                                    gagcaccaatcagaaccggtatggctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg                                                                                         <  2:450700‑M2/59‑1 (MQ=255)
                                                     aTTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa                                                                                        <  1:310131/149‑1 (MQ=255)
                                                         ccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATGTa                                                                                    >  1:320252‑M2/86‑149 (MQ=255)
                                                                         tagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTa                                                                    <  2:320252‑M2/80‑1 (MQ=255)
                                                                            gCGGAGGGGCAGAAGGCCACGGAACCGGGACACCAAAAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGt                                                                 <  2:139435/149‑1 (MQ=255)
                                                                                  aaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGt                                                           >  2:380508‑M2/61‑149 (MQ=255)
                                                                                        tgcgtcgtaatgtagccgttgcttcgttgagttttgctccttctcccaaaccaaCGTTTATGAATATGAAGAAATAACAAGCCAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaac                                                     <  1:464423‑M2/95‑1 (MQ=255)
                                                                                        tgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaac                                                     <  1:103289‑M2/95‑1 (MQ=255)
                                                                                                   gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGa                                          >  1:468782‑M2/44‑149 (MQ=255)
                                                                                                      gccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAaca                                       >  1:383477‑M2/41‑149 (MQ=255)
                                                                                                                 attgagtgttgctccttctcccaatccaaCGTTTATGAAACTGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCCGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACaa                            <  1:416033‑M2/120‑1 (MQ=255)
                                                                                                                  ttgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAAc                           >  1:242716‑M2/29‑149 (MQ=255)
                                                                                                                   gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACACCg                          >  1:364382/1‑149 (MQ=255)
                                                                                                                               cttctcccaaaccaaCGTTTATGAAAATGAAGAAATAGCAAGCAAATGGCATCATTCCTGCTTTTACCAGTGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGaa              <  1:474526‑M2/134‑1 (MQ=255)
                                                                                                                               cttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGaa              <  2:22914‑M2/134‑1 (MQ=255)
                                                                                                                                           caaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGACTGTGTACGCTTTTCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATg  >  1:504347‑M2/4‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
GTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATG  >  NC_000913_3_vch_tpiA/1209477‑1209764

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.