New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 11450531 (0.910)6 (0.200) 6/220 3.0 16.9% intergenic (+18/+17) guaC/hofC GMP reductase/assembly protein in type IV pilin biogenesis, transmembrane protein
?NC_000913 = 114511 31 (1.010)intergenic (+24/+11) guaC/hofC GMP reductase/assembly protein in type IV pilin biogenesis, transmembrane protein
Rejected: Frequency below/above cutoff threshold.

CGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATCCGC  >  NC_000913/114373‑114650
                                                                                                                                    |                                                                                                                                                 
cggcCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAAcacgccagcg                                                                                                                                   >  7:84911‑M2/1‑139 (MQ=255)
              cTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCAc                                                                                                                     >  4:45844/1‑149 (MQ=255)
                     aTACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGTTGCTCcacgccagcgttggga                                                                                                                             <  8:59178‑M2/134‑17 (MQ=255)
                     aTACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGTTGCTCcacgccagcgttggga                                                                                                                             >  7:59178‑M2/1‑118 (MQ=255)
                              ggCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTa                                                                                                      >  2:51033/1‑148 (MQ=255)
                                   cACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATTTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATg                                                                                                  <  2:71063/147‑1 (MQ=255)
                                   cACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATgg                                                                                                 <  1:51033/148‑1 (MQ=255)
                                      gCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTG                                                                                                                                                   >  8:9757‑M1/1‑95 (MQ=255)
                                         tGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGTTGCTCcacgccagcgttgggagattacaggttgttgaagatgcggttttcttgttc                                                                                          <  6:146985‑M2/149‑52 (MQ=255)
                                           aaaGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGTTGCTCcacgccagcgttgggagattacaggttgttgaagatgcggttttcttgttcct                                                                                        >  5:146985‑M2/1‑96 (MQ=255)
                                              gagCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTaaa                                                                                                    <  4:107354/134‑1 (MQ=255)
                                              gagCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTaaa                                                                                                    >  3:107354/1‑134 (MQ=255)
                                                gCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATa                                                                                   <  3:45844/149‑1 (MQ=255)
                                                    aCCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACa                                                                                 <  3:84020/147‑1 (MQ=255)
                                                                acGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGgagcaacacgccagcgttgggagattacaggttgttgaagatgcggtt                                                                                                   <  1:131278‑M1/117‑49 (MQ=255)
                                                                acGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGgagcaacacgccagcggtgggagattacaggttgttgaagatgcggtt                                                                                                   >  2:131278‑M1/1‑69 (MQ=255)
                                                                     tATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTAc                                                               >  8:109858/1‑148 (MQ=255)
                                                                          gtgtGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCaataa                                                         <  7:109858/149‑1 (MQ=255)
                                                                              gCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATtc                                                      >  2:71141/1‑148 (MQ=255)
                                                                                aGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGgagcaacacgccagcgttgggagattacaggttgttgaagatgcggttttcttgttcctgcacacgaataaacgtggtgcgcttggtcagctcttt                                                   >  5:146985‑M1/1‑53 (MQ=255)
                                                                                aGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAAt                                                                                                   <  5:57069/101‑1 (MQ=255)
                                                                                aGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAAt                                                                                                   >  6:57069/1‑101 (MQ=255)
                                                                                  gAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGgagcaacacgccagcgttgggagattacaggttgttgaagatgcggttttcttgttcctgcacacgaataaacgtggtgcgcttggtcagctctttca                                                 <  6:146985‑M1/149‑99 (MQ=255)
                                                                                     cAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGaa                                                                                               >  7:105541/1‑100 (MQ=255)
                                                                                     cAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGaa                                                                                               <  8:105541/100‑1 (MQ=255)
                                                                                          aaaCCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGAt                                         <  1:71141/149‑1 (MQ=255)
                                                                                           aaCCGCATCTTCAACAACCTGTAATCTCCCACCGCTGGCGTGTTGCTCcacgccagcgttgggagattacaggttgttgaagatgcggttttcttgttcctgcacacgaataaacgt                                                                        >  2:131278‑M2/1‑48 (MQ=255)
                                                                                           aaCCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGTTGCTCcacgccagcgttgggagattacaggttgttgaagatgcggttttcttgttcctgcacacgaataaacgt                                                                        <  1:131278‑M2/117‑70 (MQ=255)
                                                                                               gCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGccaccacca                                                                     >  3:57536/1‑116 (MQ=255)
                                                                                               gCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGccaccacca                                                                     <  4:57536/116‑1 (MQ=255)
                                                                                                   cTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAAtt                                                       >  1:78170/1‑126 (MQ=255)
                                                                                                   cTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAAtt                                                       <  2:78170/126‑1 (MQ=255)
                                                                                                    ttCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGcc                               >  1:189662/1‑149 (MQ=255)
                                                                                                           acCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTc                          >  3:63596/1‑147 (MQ=255)
                                                                                                                     tcCCAACGCTGGCGTGgagcaacacgccagcgttgggagattacaggttgttgaagatgcggttttcttgttcctgcacacgaataaacgtggtgcgcttggtcagctctttcaggcgtgaagccccaacgtat                             >  7:59178‑M1/1‑16 (MQ=255)
                                                                                                                     tcCCAACGCTGGCGTGgagcaacacgccagcgttgggagattacaggttgttgaagatgcggttttcttgttcctgcacacgaataaacgtggtgcgcttggtcagctctttcaggcgtgaagccccaacgtat                             <  8:59178‑M1/134‑119 (MQ=255)
                                                                                                                      cccAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGc                                                                                        <  4:129558/74‑1 (MQ=255)
                                                                                                                      cccAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGc                                                                                        >  3:129558/1‑74 (MQ=255)
                                                                                                                         aaCGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCaat                                                           >  4:11981/1‑100 (MQ=255)
                                                                                                                         aaCGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCaat                                                           <  3:11981/100‑1 (MQ=255)
                                                                                                                           cGCTGGCGTGttgctccacgccagcgttgggagattacaggttgttgaagatgcggttttcttgttcctgcacacgaataaacgtggtgcgcttggtcagctctttcaggcgtgaagccccaacgtatgtacaagctgaacgcaggccg        >  7:84911‑M1/1‑10 (MQ=255)
                                                                                                                           cGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGt         >  2:161752/1‑148 (MQ=255)
                                                                                                                              tGGCGTGTTGCTCCACGCCAGCGTAAACCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCACGTTATGATCAGCaa                                   >  7:150426/1‑119 (MQ=255)
                                                                                                                              tGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAAt                                                                                                   <  2:171047/55‑1 (MQ=255)
                                                                                                                              tGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAAt                                                                                                   >  1:171047/1‑55 (MQ=255)
                                                                                                                                 cGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATCCGc  <  1:161752/149‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                                 
CGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATCCGC  >  NC_000913/114373‑114650

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