New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 1647433 =27 (0.860)5 (0.170) 5/278 3.4 15.2% coding (1289/1335 nt) oprB‑II carbohydrate‑selective porin
?NC_002947 1647473 = 30 (1.000)coding (1249/1335 nt) oprB‑II carbohydrate‑selective porin
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GGCTTCCTTGGTTCGCACCTTCAGTGCTCATCGTAGGTTCTCCGTCAGGTCAAAAAAGGGGTTAGAAACTGCTCTGGATCTTCAGGCCGCCGATCAGGGCGCCATCGACCTGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG  >  NC_002947/1647326‑1647612
                                                                                                                                                   |                                                                                                                                           
ggCTTCCTTGGTTCGCACCTTCAGTGCTCATCGTAGGTTCTCCGTCAGGTCAAAAAAGGGGTTAGAAACTGCTCTGGATCTTCAGGCCGCCGATCAGGGCGCCATCGACCTGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGtt                                                                                                                                             >  1:13039/1‑148 (MQ=255)
   ttCCTTGGTTCGCACCTTCAGTGCTCATCGTAGGTTCTCCGTCAGGTCAAAAAAGGGGTTAGAAACTGCTCTGGATCTTCAGGCCGCCGATCAGGGCGCCATCGACCTGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTggg                                                                                                                                          >  5:161439/1‑148 (MQ=255)
          gTTCGCACCTTCAGTGCTCATCGTAGGTTCTCCGTCAGGTCAAAAAAGGGGTTAGAAACTGCTCTGGATCTTCAGGCCGCCGATCAGGGCGCCATCGACCTGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTAc                                                                                                                                    <  2:112345/147‑1 (MQ=255)
               cACCTTCAGTGCTCATCGTAGGTTCTCCGTCAGGTCAAAAAAGGGGTTAGAAACTGCTCTGGATCTTCAGGCCGCCGATCAGGGCGCCATCGACCTGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGa                                                                                                                               <  8:40418/147‑1 (MQ=255)
                aCCTTCAGTGCTCATCGTAGGTTCTCCGTCAGGTCAAAAAAGGGGTTAGAAACTGCTCTGGATCTTCAGGCCGCCGATCAGGGCGCCATCGACCTGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGc                                                                                                                             <  1:18072/148‑1 (MQ=255)
                  cTTCAGTGCGCATCGTAGGTTCTCCGTCAGGTCAAAAAAGGGGTTAGAAACTGCACTGGATCTTCAGGCCGCCGCTCAGGGCGCCATCGACCTGCGCCACCACCCCCGGGTGGCGCATGTACTGCAGGATGGGGCGTCCCGTGAGCCa                                                                                                                           >  8:121712/1‑148 (MQ=255)
                           tCATCGTAGGTTCTCCGTCAGGTCAAAAAAGGGGTTAGAAACTGCTCTGGATCTTCAGGCCGCCGATCAGGGCGCCATCGACCTGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCaa                                                                                                                  <  4:145917/148‑1 (MQ=255)
                                  aGGTTCTCCGTCAGGTCAAAAAAGGGGTTAGAAACTGCTCTGGATCTTCAGGCCGCCGATCAGGGCGCCATCGACCTGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATg                                                                                                           >  3:158221/1‑148 (MQ=255)
                                           gTCAGGTCAAAAAAGGGGTTAGAAACTGCTCTGGATCTTCAGGCCGCCGATCAGGGCGCCATCGACCTGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATa                                                                                                   <  1:134731/147‑1 (MQ=255)
                                            tCAGGTCAAAAAAGGGGTTAGAAACTGCTCTGGATCTTCAGGCCGCCGATCAGGGCGCCATCGACCTGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACa                                                                                                 >  3:72174/1‑148 (MQ=255)
                                                 tCAAAAAAGGGGTTAGAAACTGCTCTGGATCTTCAGGCCGCCGATCAGGGCGCCATCGACCTGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCg                                                                                            >  1:11987/1‑148 (MQ=255)
                                                            gTTAGAAACTGCTCTGGATCTTCAGGCCGCCGATCAGGGCGCCATCGACCTGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTc                                                                                 >  1:32642/1‑148 (MQ=255)
                                                                 aaaCTGCTCTGGATCTTCAGGCCGCCGATCAGGGCGCCATCGACCTGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGgtg                                                                             <  4:158221/147‑1 (MQ=255)
                                                                    cTGCTCTGGATCTTCAGGCCGCCGATCAGGGCGCCATCGACCTGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCt                                                                         <  6:56910/148‑1 (MQ=255)
                                                                        tctGGATCTTCAGGCCGCCGATCAGGGCGCCATCGACCTGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCa                                                                      <  4:66499/147‑1 (MQ=255)
                                                                            gATCTTCAGGCCGCCGATCAGGGCGCCATCGACCTGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTg                                                                   >  6:126889/1‑146 (MQ=255)
                                                                                       cgccgATCAGGGCGCCATCGACCTGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCgg                                                      <  2:73817/148‑1 (MQ=255)
                                                                                       cgccgATCAGGGCGCCATCGACCTGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCgg                                                      >  3:160550/1‑148 (MQ=255)
                                                                                          cgATCAGGGCGCCATCGACCTGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGtt                                                   <  1:72523/148‑1 (MQ=255)
                                                                                           gATCAGGGCGCCATCGACCTGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTg                                                  >  8:57481/1‑148 (MQ=255)
                                                                                                  gcgcCATCGACCTGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGt                                           <  5:126889/148‑1 (MQ=255)
                                                                                                      cATCGACCTGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTc                                          >  3:35727/1‑145 (MQ=255)
                                                                                                                     ccccGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTCATAGAGGCCGGCGGCCTg                        <  4:151553/148‑1 (MQ=255)
                                                                                                                      cccgcccgGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTCATAGAGGCCGGCGGCCTgg                       >  5:166126/1‑148 (MQ=255)
                                                                                                                      cccgcccgGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTCATAGAGGCCGGCGGCCTgg                       <  1:180750/148‑1 (MQ=255)
                                                                                                                             gggTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTCATAGAGGCCGGCGGCCTGGTTGACCa                <  6:161439/148‑1 (MQ=255)
                                                                                                                                   cGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTCATAGAGGCCGGCGGCCTGGTTGACCAggcgg           <  3:79841/147‑1 (MQ=255)
                                                                                                                                   cGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTCATAGAGGCCGGCGGCCTGGTTGACCAggcgg           <  3:77191/147‑1 (MQ=255)
                                                                                                                                          aCTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCtt   <  7:57481/148‑1 (MQ=255)
                                                                                                                                            tGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTg  >  6:167350/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
GGCTTCCTTGGTTCGCACCTTCAGTGCTCATCGTAGGTTCTCCGTCAGGTCAAAAAAGGGGTTAGAAACTGCTCTGGATCTTCAGGCCGCCGATCAGGGCGCCATCGACCTGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG  >  NC_002947/1647326‑1647612

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.