New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2004132 =	34 (1.010)	6 (0.180)	6/262	3.1	14.3%	coding (261/1407 nt)	fliD	flagellar filament capping protein
	?	NC_000913	2004173 =	39 (1.190)					coding (302/1407 nt)	fliD	flagellar filament capping protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCTGACGACTTTCCAGACCGCCAATACTGCATTGTCTAAAGCCGATCTTTTTTCCGCCACCAGCACCACCAGCAGCACCACCGCGTTCAGTGCCACCACTGCGGGTAACGCCATCGCCGGGAAATACACCATCAGCGTCACCCATCTGGCGCAGGCGCAAACCCTGACCACGCGCACCACCAGAGACGATACGAAAACGGCGATCGCCACCAGCGACAGTAAACTCACCATTCAACAAGGCGG  >  NC_000913/2003993‑2004272                                                                                                                                            |                                                                                                                                             ccGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCTGACGACTTTCCAGACCGCCAATACTGCATTGTCTAAAGCCGATCTTTTTTCCGCCACCAGCACCACCAGCAGCACCACCGCGTTCAGTGCCACCACTGCGGGTAACGc                                                                                                                                      >  5:27739/1‑148 (MQ=255)      aaaCTTAGCGCCTACGGAACGCTGAAAAGCGCGCTGACGACTTTCCAGACCGCCAATACTGCATTGTCTAAAGCCGATCTTTTTTCCGCCACCAGCACCACCAGCAGCACCACCGCGTTCAGTGCCACCACTGCGGGTAACGCCATCGc                                                                                                                                <  1:161232/149‑1 (MQ=255)             gcgcCTACGGTACGCTGAAAAGCGCGCTGACGACTTTCCAGACCGCCAATACTGCATTGTCTAAAGCCGATCTTTTTTCCGCCACCAGCACCACCAGCAGCACCACCGCGTTCAGTGCCACCACTGCGGGTAACGCCATCGCCGGGaaa                                                                                                                         <  5:107539/149‑1 (MQ=255)                       tACGCTGAAAAGCGCGCTGACGACTTTCCAGACCGCCAATACTGCATTGTCTAAAGCCGATCTTTTTTCCGCCACCAGCACCACCAGCAGCACCACCGCGTTCAGTGCCACCACTGCGGGTAACGCCATCGCCGGGAAATACACCATCa                                                                                                               <  6:54051/149‑1 (MQ=255)                         cGCTGAAAAGCGCGCTGACGACTTTCCAGACCGCCAATACTGCATTGTCTAAAGCCGATCTTTTTTCCGCCACCAGCACCACCAGCAGCACCACCGCGTTCAGTGCCACCACTGCGGGTAACGCCATCGCCGGGAAATACACCATCAGc                                                                                                             <  6:27739/149‑1 (MQ=255)                                               tttCCAGACCGCCAATACTGCATTGTCTAAAGCCGATCTTTTTTCCGCCACCAGCACCACCAGCAGCACCACCGCGTTCAGTGCCACCACTGCGGGTAACGCCATCGCCGGGAAATACACCATCAGCGTCACCCATCTggcgcaggcg                                                                                        >  8:128101/1‑148 (MQ=255)                                               tttCCAGACCGCCAATACTGCATTGTCTAAAGCCGATCTTTTTTCCGCCACCAGCACCACCAGCAGCACCACCGCGTTCAGTGCCACCACTGCGGGTAACGCCATCGCCGGGAAATACACCATCAGCGTCACCCATCTggcgcaggcg                                                                                        >  4:40328/1‑148 (MQ=255)                                                    aGACCGCCAATACTGCATTGTCTAAAGCCGATCTTTTTTCCGCCACCAGCACCACCAGCAGCACCACCGCGTTCAGTGCCACCACTGCGGGTAACGCCATCGCCGGGAAATACACCATCAGCGTCACCCATCTGGCGCAGGCGCAAAcc                                                                                  >  1:84470/1‑149 (MQ=255)                                                    aGACCGCCAATACTGCATTGTCTAAAGCCGATCTTTTTTCCGCCACCAGCACCACCAGCAGCACCACCGCGTTCAGTGCCACCACTGCGGGTAACGCCATCGCCGGGAAATACACCATCAGCGTCACCCATCTGGCGCAGGCGCAAAcc                                                                                  >  2:27607/1‑149 (MQ=255)                                                      aCCGCCAATACTGCATTGTCTAAAGCCGATCTTTTTTCCGCCACCAGCACCACCAGCAGCACCACCGCGTTCAGTGCCACCACTGCGGGTAACGCCATCGCCGGGAAATACACCATCAGCGTCACCCATCTGGCGCAGGCGCAAACCCt                                                                                <  7:128101/149‑1 (MQ=255)                                                      aCCGCCAATACTGCATTGTCTAAAGCCGATCTTTTTTCCGCCACCAGCACCACCAGCAGCACCACCGCGTTCAGTGCCACCACTGCGGGTAACGCCATCGCCGGGAAATACACCATCAGCGTCACCCATCTGGCGCAGGCGCAAACCCt                                                                                <  1:27607/149‑1 (MQ=255)                                                        cGCCAATACTGCATTGTCTAAAGCCGATCTTTTTTCCGCCACCAGCACCACCAGCAGCACCACCGCGTTCAGTGCCACCACTGCGGGTAACGCCATCGCCGGGAAATACACCATCAGCGTCACCCATCTGGCGCAGGCGCAAACCCTGa                                                                              >  8:26467/1‑149 (MQ=255)                                                         gCCAATACTGCATTGTCTAAAGCCGATCTTTTTTCCGCCACCAGCACCACCAGCAGCACCACCGCGTTCAGTGCCACCACTGCGGGTAACGCCATCGCCGGGAAATACACCATCAGCGTCACCCATCTGGCGCAGGCGCAAACCCTGAc                                                                             >  1:84220/1‑149 (MQ=255)                                                           cAATACTGCATTGTCTAAAGCCGATCTTTTTTCCGCCACCAGCACCACCAGCAGCACCACCGCGATCAGTGCCACCACTGCCGGTAAAGCCATCGCCGGGATATACACCATCAGCGTCCCACATCTGGCGCAGGCGCAACCCCTGACCa                                                                           >  3:96896/1‑149 (MQ=255)                                                                               ccGATCTTTTTTCCGCCACCAGCACCACCAGCAGCACCACCGCGTTCAGTGCCACCACTGCGGGTAACGCCATCGCCGGGAAATACACCATCAGCGTCACCCATCTGGCGCAGGCGCAAACCCTGACCACGCGCACCACCAGAGACGa                                                        >  7:21946/1‑148 (MQ=255)                                                                                  aTCTTTTTTCCGCCACCAGCACCACCAGCAGCACCACCGCGTTCAGTGCCACCACTGCGGGTAACGCCATCGCCGGGAAATACACCATCAGCGTCACCCATCTGGCGCAGGCGCAAACCCTGACCACGCGCACCACCAGAGACGATa                                                      <  2:117782/147‑1 (MQ=255)                                                                                  aTCTTTTTTCCGCCACCAGCACCACCAGCAGCACCACCGCGTTCAGTGCCACCACTGCGGGTAACGCCATCGCCGGGAAATACACCATCAGCGTCACCCATCTGGCGCAGGCGCAAACCCTGACCACGCGCACCACCAGAGACGATa                                                      >  1:117782/1‑147 (MQ=255)                                                                                    cTTTTTTCCGCCACCAGCACCACCAGCAGCACCACCGCGTTCAGTGCCACCACTGCGGGTAACGCCATCGCCGGGAAATACACCATCAGCGTCACCCATCTGGCGCAGGCGCAAACCCTGAc                                                                             >  2:85099/1‑122 (MQ=255)                                                                                     ttttttCCGCCACCAGCACCACCAGCAGCACCACCGCGTTCAGTGCCCCCACTGCGGGTAACGCCATCGCCGGGAAATACACCATCAGCGTCACCCATCTGGCGCAGGCGCAAACCCTGACCACGCGCACCACCAGAGACGATACGaa                                                  >  5:35833/1‑148 (MQ=255)                                                                                       ttttCCGCCACCAGCACCACCAGCAGCACCACCGCGTTCAGTGCCACCACTGCGGGTAACGCCATCGCCGGGAAATACACCATCAGCGTCACCCATCTGGCGCAGGCGCAAACCCTGACCACGCGCACCACCAGAGACGATAc                                                     <  2:84470/143‑1 (MQ=255)                                                                                              ccaccagcaccaccagcaGCACCACCGCGTTCAGTGCCACCACTGCGGGTAACGCCATCGCCGGGAAATACACCATCAGCGTCACCCATCTGGCGCAGGCGCAAACCCTGACCACGCGCACCACCAGAGACGATACGAAAACGGCGATc                                        <  2:84220/149‑1 (MQ=255)                                                                                                 ccagcaccaccagcaGCACCACCGCGTTCAGTGCCACCACTGCGGGTAACGCCATCGCCGGGAAATACACCATCAGCGTCACCCATCTGGCGCAGGCGCAAACCCTGa                                                                              <  7:35341/108‑1 (MQ=255)                                                                                                 ccagcaccaccagcaGCACCACCGCGTTCAGTGCCACCACTGCGGGTAACGCCATCGCCGGGAAATACACCATCAGCGTCACCCATCTGGCGCAGGCGCAAACCCTGa                                                                              >  8:35341/1‑108 (MQ=255)                                                                                                     caccaccagcaGCACCACCGCGTTCAGTGCCACCACTGCGGGTAACTCCATCGCCGGGAAATACACCATCAGCGTCACCCATCTGGCGCAGGCGCAAACCCTGACCACGCGCACCACCAGAGACGATACGAAAACGGCGATCGccacca                                 <  8:21946/149‑1 (MQ=255)                                                                                                        caccagcaGCACCACCGCGTTCAGTGCCACCACTGCGGGTAACGCCATCGCCGGGAAATACACCATCAGCGTCACCCATCTggc                                                                                               >  2:118258/1‑84 (MQ=255)                                                                                                        caccagcaGCACCACCGCGTTCAGTGCCACCACTGCGGGTAACGCCATCGCCGGGAAATACACCATCAGCGTCACCCATCTGGCGCAGGCGCAAACCCTGACCACGCGca                                                                     <  2:123012/110‑1 (MQ=255)                                                                                                        caccagcaGCACCACCGCGTTCAGTGCCACCACTGCGGGTAACGCCATCGCCGGGAAATACACCATCAGCGTCACCCATCTGGCGCAGGCGCAAACCCTGACCACGCGca                                                                     >  1:123012/1‑110 (MQ=255)                                                                                                        caccagcaGCACCACCGCGTTCAGTGCCACCACTGCGGGTAACGCCATCGCCGGGAAATACACCATCAGCGTCACCCATCTGGCGCAGGCGCAAACCCTGACCACGCGCACCACCAGAGACGATACGAAAACGGCGATCGCCACCAGc                               >  3:46355/1‑148 (MQ=255)                                                                                                                    caccGCGTTCAGTGCCACCACTGCGGGTAACGCCATCGCCGGGAAATACACCATCAGCGTCACCCATCTGGCGCAGGCGCAAACCCTGACCACGCGCACCACCAGAGACGATACGAAAACGGCGATCGCCACCAGCGACAGTAAACTCa                  <  3:40328/149‑1 (MQ=255)                                                                                                                       cgcgTTCATTGCCACCACTGCGGTTACCGCCATCGCCGTGAACTACACCATCAGCTTCACCCAACTGGCTCATGAGCATATCCTGACCACGCGCACCACCAGTGACGATCCGAATACGTCTTGCGCCACCAGCGACAGTAAACTCACCa               <  4:96896/149‑1 (MQ=255)                                                                                                                               gTGCCACCACTGCGGGTAACGCCATCGCCGGGAAATACACCATCAGCGTCACCCATCTGGCGCAGGCGCAAACCCTGACCACGCGCACCACCAGAGACGATACGAAAACGGCGATCGCCACCAGCGACAGTAAACTCACCATTcaac         >  3:27481/1‑147 (MQ=255)                                                                                                                               gTGCCACCACTGCGGGTAACGCCATCGCCGGGAAATACACCATCAGCGTCACCCATCTGGCGCAGGCGCAAACCCTGACCACGCGCACCACCAGAGACGATACGAAAACGGCGATCGCCACCAGCGACAGTAAACTCACCATTc            <  2:122245/144‑1 (MQ=255)                                                                                                                               gTGCCACCACTGCGGGTAACGCCATCGCCGGGAAATACACCATCAGCGTCACCCATCTGGCGCAGGCGCAAACCCTGACCACGCGCACCACCAGAGACGATACGAAAACGGCGATCGCCACCAGCGACAGTAAACTCACCATTc            >  1:122245/1‑144 (MQ=255)                                                                                                                                    accacTGCGGGTAACGCCATCGCCGGGAAATACACCATCAGCGTCACCCATCTGGCGCAGGCGCAAACCCTGACCACGCGCACCACCAGAGACGATACGAAAACGGCGATCGCCACCAGCGACAGTAAACTCACCATTCAACAAggcgg  >  5:8751/1‑149 (MQ=255)                                                                                                                                            |                                                                                                                                             CCGCTAAACTTAGCGCCTACGGTACGCTGAAAAGCGCGCTGACGACTTTCCAGACCGCCAATACTGCATTGTCTAAAGCCGATCTTTTTTCCGCCACCAGCACCACCAGCAGCACCACCGCGTTCAGTGCCACCACTGCGGGTAACGCCATCGCCGGGAAATACACCATCAGCGTCACCCATCTGGCGCAGGCGCAAACCCTGACCACGCGCACCACCAGAGACGATACGAAAACGGCGATCGCCACCAGCGACAGTAAACTCACCATTCAACAAGGCGG  >  NC_000913/2003993‑2004272

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