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.

AGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCTGGATAGAGAGTCAGACGACTTTTGCGGCCGCTGACGATATCGTAACTGTCGAGTGAGTTTTTGCTGTTCTGTAACTCCACCTCA  >  NC_000913/306926‑307215
                                                                                                                                                    |                                                                                                                                             
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)
                                                                                                              cGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGTCGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCTGGATAGAGAGTCAGACGACTTTTGCGGCCGGTGACGATATCGTAACTGTCga                                  <  8:135128/148‑1 (MQ=255)
                                                                                                                 gCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCTGGATAGAGAGTCAGACGACTTTTGCGGCCGCTGACGATATCGTAACTGTCgagtga                              <  4:74354/149‑1 (MQ=255)
                                                                                                                        ccAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCTGGATAGAGAGTCAGACGACTTTTGCGGCCGCTGACGATATCGTAACTGTCGAGTGAGTTTTTg                       <  5:94129/149‑1 (MQ=255)
                                                                                                                         cagcagTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCTGGATAGAGAGTCa                                                                         >  3:4657/1‑98 (MQ=255)
                                                                                                                         cagcagTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCTGGATAGAGAGTCa                                                                         <  4:4657/98‑1 (MQ=255)
                                                                                                                           gcagTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCTGGATagagag                                                                            <  6:29664/93‑1 (MQ=255)
                                                                                                                           gcagTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCTGGATagagag                                                                            >  5:29664/1‑93 (MQ=255)
                                                                                                                                             cACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCTGGATAGAGAGTCAGACGACTTTTGCGGCCGCTGACGATATCGTAACTGTCGAGTGAGTTTTTGCTGTTCTGTAACTCCACCTCa  >  8:58940/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
AGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCTGGATAGAGAGTCAGACGACTTTTGCGGCCGCTGACGATATCGTAACTGTCGAGTGAGTTTTTGCTGTTCTGTAACTCCACCTCA  >  NC_000913/306926‑307215

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.