New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 250165 =32 (0.940)7 (0.220) 6/230 3.1 18.4% pseudogene (94/756 nt) lafU pseudogene, lateral flagellar motor protein fragment
?NC_000913 250214 = 32 (1.000)pseudogene (143/756 nt) lafU pseudogene, lateral flagellar motor protein fragment
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TCGGAAATACGCTCGGCCCCTTTGGTGACGACAATAAAGTTGATTGTCAACAGCGTCAGTAAATCCGAGCGTGAAAGCATTATCGCCGCGCTGCACGGGCAGTCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATC  >  NC_000913/250018‑250299
                                                                                                                                                   |                                                                                                                                      
tCGGAAATACGCTCGGCCCCTTTGGTGACGACAATAAAGTTGATTGTCAACAGCGTCAGTAAATCCGAGCGTGAAAGCATTATCGCCGCGCTGCACGGGCAGTCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGccc                                                                                                                                       >  2:130550/1‑149 (MQ=255)
   gAAATACGCTCGGCCCCTTTGGTGACGACAATAAAGTTGATTGTCAACAGCGTCAGTAAATCCGAGCGTGAAAGCATTATCGCCGCGCTGCACGGGCAGTCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCCCGtc                                                                                                                                    <  1:130550/149‑1 (MQ=255)
                 cccTTTGGTGACGACAATAAAGTTGATTGTCAACAGCGTCAGTAAATCCGAGCGTGAAAGCATTATCGCCGCGCTGCACGGGCAGTCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAAc                                                                                                                      <  7:42678/149‑1 (MQ=255)
                             gacAATAAAGTTGATTGTCAACAGCGTCAGTAAATCCGAGCGTGAAAGCATTATCGCCGCGCTGCACGGGCAGTCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTgg                                                                                                          >  1:15588/1‑149 (MQ=255)
                             gacAATAAAGTTGATTGTCAACAGCGTCAGTAAATCCGAGCGTGAAAGCATTATCGCCGCGCTGCACGGGCAGTCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTgg                                                                                                          >  3:2739/1‑149 (MQ=255)
                                  tAAAGTTGATTGTCAACAGCGTCAGTAAATCCGAGCGTGAAAGCATTATCGCCGCGCTGCACGGGCAGTCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCTCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTg                                                                                                     >  6:141049/1‑149 (MQ=255)
                                   aaaGTTGATTGTCAACAGCGTCAGTAAATCCGAGCGTGAAAGCATTATCGCCGCGCTGCACGGGCAGTCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGc                                                                                                    >  3:128422/1‑149 (MQ=255)
                                              tCAACAGCGTCAGTAAATCCGAGCGTGAAAGCATTATCGCCGCGCTGCACGGGCAGTCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACg                                                                                         <  4:128422/149‑1 (MQ=255)
                                              tCAACAGCGTCAGTAAATCCGAGCGTGAAAGCATTATCGCCGCGCTGCACGGGCAGTCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACg                                                                                         >  5:58345/1‑149 (MQ=255)
                                                aaCAGCGTCAGTAAATCCGAGCGTGAAAGCATTATCGCCGCGCTGCACGGGCAGTCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGga                                                                                       <  3:91206/149‑1 (MQ=255)
                                                  cAGCGTCAGTAAATCCGAGCGTGAAAGCATTATCGCCGCGCTGCACGGGCAGTCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGgag                                                                                      >  7:38024/1‑148 (MQ=255)
                                                   aGCGTCAGTAAATCCGAGCGTGAAAGCATTATCGCCGCGCTGCACGGGCAGTCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAgaga                                                                                           <  3:4252/142‑1 (MQ=255)
                                                   aGCGTCAGTAAATCCGAGCGTGAAAGCATTATCGCCGCGCTGCACGGGCAGTCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAgaga                                                                                           >  4:4252/1‑142 (MQ=255)
                                                          gTAAATCCGAGCGTGAAAGCATTATCGCCGCGCTGCACGGGCAGTCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAgcg                                                                             >  1:2493/1‑149 (MQ=255)
                                                               tCCGAGCGTGAAAGCATTATCGCCGCGCTGCACGGGCAGTCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGa                                                                        >  2:190494/1‑149 (MQ=255)
                                                                ccGAGCGTGAAAGCATTATCGCCGCGCTGCACGGGCAGTCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAgcgc                                                                            >  2:166415/1‑144 (MQ=255)
                                                                ccGAGCGTGAAAGCATTATCGCCGCGCTGCACGGGCAGTCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAgcgc                                                                            <  1:166415/144‑1 (MQ=255)
                                                                  gAGCGTGAAAGCATTATCGCCGCGCTGCACGGGCAGTCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCTCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGAtgt                                                                     <  5:141049/149‑1 (MQ=255)
                                                                  gAGCGTGAAAGCATTATCGCCGCGCTGCACGGGCAGTCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGAtgt                                                                     <  8:38024/149‑1 (MQ=255)
                                                                      gTGAAAGCATTATCGCCGCGCTGCACGGGCAGTCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAAt                                                                 <  4:2739/149‑1 (MQ=255)
                                                                        gAAAGCATTATCGCCGCGCTGCACGGGCAGTCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAGATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAAt                                                                 >  7:8491/1‑147 (MQ=255)
                                                                         aaaGCATTATCGCCGCGCTGCACGGGCAGTCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAATga                                                               <  2:2493/148‑1 (MQ=255)
                                                                                                gggCAGTCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTgaagaaga                                        >  8:46540/1‑148 (MQ=255)
                                                                                                      tCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGc                                 >  7:31114/1‑149 (MQ=255)
                                                                                                      tCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGc                                 <  2:15588/149‑1 (MQ=255)
                                                                                                                aCGGTGGTGGGCTGTCGCCGTTAGATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCtt                       <  8:8491/149‑1 (MQ=255)
                                                                                                                 cGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTg                      <  7:46540/149‑1 (MQ=255)
                                                                                                                       tggGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAgcg                                                                             >  5:156004/1‑88 (MQ=255)
                                                                                                                       tggGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAgcg                                                                             <  6:156004/88‑1 (MQ=255)
                                                                                                                              tCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCa           <  1:190494/147‑1 (MQ=255)
                                                                                                                                  cGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAg      >  3:13077/1‑148 (MQ=255)
                                                                                                                                      aaataaaATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATc  <  8:31114/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                      
TCGGAAATACGCTCGGCCCCTTTGGTGACGACAATAAAGTTGATTGTCAACAGCGTCAGTAAATCCGAGCGTGAAAGCATTATCGCCGCGCTGCACGGGCAGTCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATC  >  NC_000913/250018‑250299

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