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.

AACAGCGTCAGTAAATCCGAGCGTGAAAGCATTATCGCCGCGCTGCACGGGCAGTCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCA  >  NC_000913/250066‑250358
                                                                                                                                                    |                                                                                                                                                
aaCAGCGTCAGTAAATCCGAGCGTGAAAGCATTATCGCCGCGCTGCACGGGCAGTCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGga                                                                                                                                                  <  3:91206/149‑1 (MQ=255)
  cAGCGTCAGTAAATCCGAGCGTGAAAGCATTATCGCCGCGCTGCACGGGCAGTCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGgag                                                                                                                                                 >  7:38024/1‑148 (MQ=255)
          gTAAATCCGAGCGTGAAAGCATTATCGCCGCGCTGCACGGGCAGTCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAgcg                                                                                                                                        >  1:2493/1‑149 (MQ=255)
               tCCGAGCGTGAAAGCATTATCGCCGCGCTGCACGGGCAGTCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGa                                                                                                                                   >  2:190494/1‑149 (MQ=255)
                ccGAGCGTGAAAGCATTATCGCCGCGCTGCACGGGCAGTCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAgcgc                                                                                                                                       <  1:166415/144‑1 (MQ=255)
                ccGAGCGTGAAAGCATTATCGCCGCGCTGCACGGGCAGTCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAgcgc                                                                                                                                       >  2:166415/1‑144 (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                                                                                            <  2:15588/149‑1 (MQ=255)
                                                      tCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGc                                                                                            >  7:31114/1‑149 (MQ=255)
                                                                aCGGTGGTGGGCTGTCGCCGTTAGATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCtt                                                                                  <  8:8491/149‑1 (MQ=255)
                                                                 cGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTg                                                                                 <  7:46540/149‑1 (MQ=255)
                                                                       tggGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAgcg                                                                                                                                        <  6:156004/88‑1 (MQ=255)
                                                                       tggGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAgcg                                                                                                                                        >  5:156004/1‑88 (MQ=255)
                                                                              tCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCa                                                                      <  1:190494/147‑1 (MQ=255)
                                                                                  cGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAg                                                                 >  3:13077/1‑148 (MQ=255)
                                                                                      aaataaaATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATc                                                             <  8:31114/148‑1 (MQ=255)
                                                                                                          aTCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCa                                                                                          >  2:77966/1‑99 (MQ=255)
                                                                                                          aTCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCa                                                                                          <  1:77966/99‑1 (MQ=255)
                                                                                                          aTCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATgcg                                        >  2:153740/1‑149 (MQ=255)
                                                                                                             cGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCa                                      >  2:186989/1‑148 (MQ=255)
                                                                                                             cGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCa                                      >  7:151092/1‑148 (MQ=255)
                                                                                                                         ccACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAAt                                                                                                                            <  1:48863/50‑1 (MQ=255)
                                                                                                                         ccACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAAt                                                                                                                            <  1:51557/50‑1 (MQ=255)
                                                                                                                         ccACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAAt                                                                                                                            >  2:48863/1‑50 (MQ=255)
                                                                                                                         ccACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAAt                                                                                                                            >  2:51557/1‑50 (MQ=255)
                                                                                                                            ccGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCtg                      <  6:58345/149‑1 (MQ=255)
                                                                                                                                       gcGGATGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCtggagatggaga            <  6:57883/147‑1 (MQ=255)
                                                                                                                                         ggAAGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTg         >  1:85252/1‑149 (MQ=255)
                                                                                                                                              agaCGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGCAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTccc    >  5:172928/1‑149 (MQ=255)
                                                                                                                                              agaCGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTcc     <  4:5661/148‑1 (MQ=255)
                                                                                                                                                aCGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCa  <  4:13077/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
AACAGCGTCAGTAAATCCGAGCGTGAAAGCATTATCGCCGCGCTGCACGGGCAGTCGATTTTTAACGGTGGTGGGCTGTCGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCA  >  NC_000913/250066‑250358

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