New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 307033 =27 (0.790)4 (0.130) 4/242 3.3 13.5% coding (2300/2526 nt) ecpC ECP production outer membrane protein
?NC_000913 307074 = 26 (0.820)coding (2259/2526 nt) ecpC ECP production outer membrane protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCTGGATAGAGAGTCAGACGACTTTTGCGGCCGCTGACGATATCGTAACTGT  >  NC_000913/306896‑307180
                                                                                                                                         |                                                                                                                                                   
gcgcCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTa                                                                                                                                          <  5:58477/149‑1 (MQ=255)
       gTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCcag                                                                                                                                     >  5:131530/1‑147 (MQ=255)
                cAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGt                                                                                                                          >  3:74354/1‑149 (MQ=255)
                     ccACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTcc                                                                                                                     >  7:13571/1‑149 (MQ=255)
                       aCTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTcc                                                                                                                     <  8:85234/147‑1 (MQ=255)
                       aCTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTcc                                                                                                                     >  7:85234/1‑147 (MQ=255)
                           cGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACg                                                                                                                >  3:43790/1‑148 (MQ=255)
                             cAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGAt                                                                                                             <  8:13571/149‑1 (MQ=255)
                              aGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATa                                                                                                            <  8:37625/149‑1 (MQ=255)
                              aGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCAGTCTTCCGCACGGATa                                                                                                            >  7:37625/1‑149 (MQ=255)
                               ggTTTTATTGACACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGGTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATAc                                                                                                           <  8:99180/149‑1 (MQ=255)
                                     aTTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCg                                                                                                      <  6:131530/148‑1 (MQ=255)
                                     aTTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCg                                                                                                      <  4:43790/148‑1 (MQ=255)
                                                 gCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGga                                                                                                    <  2:12471/138‑1 (MQ=255)
                                                 gCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGga                                                                                                    >  1:12471/1‑138 (MQ=255)
                                                               tGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACctct                                                                           >  4:19009/1‑149 (MQ=255)
                                                                             tCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGc                                                             <  3:19009/149‑1 (MQ=255)
                                                                               cACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATa                                                                                                            <  5:551/100‑1 (MQ=255)
                                                                               cACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATa                                                                                                            >  6:551/1‑100 (MQ=255)
                                                                                aCGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGa                                                           >  6:124299/1‑148 (MQ=255)
                                                                                       tGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATc                                                                                       <  5:106347/113‑1 (MQ=255)
                                                                                           aaaCTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCTGGa                                                >  7:135128/1‑148 (MQ=255)
                                                                                                cGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCAt                                                                                        <  4:116037/103‑1 (MQ=255)
                                                                                                cGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCAt                                                                                        >  3:116037/1‑103 (MQ=255)
                                                                                                  ccGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCTGGATAGAGAGt                                        <  5:124299/149‑1 (MQ=255)
                                                                                                                                  gttgttAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCTGGATAGAGAGTCAGACGACTTTTGCGGCCGCTGACGATATCg         >  6:142872/1‑148 (MQ=255)
                                                                                                                                        aaTCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCTGGATAGAGAGTCAGACGACTTTTGCGGCCGCTGACGATATCGTAACTGt  >  6:94129/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                                   
GCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCTGGATAGAGAGTCAGACGACTTTTGCGGCCGCTGACGATATCGTAACTGT  >  NC_000913/306896‑307180

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