New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3377925	32 (0.880)	3 (0.080)	3/222	3.5	8.4%	coding (283/393 nt)	rpsI	30S ribosomal subunit protein S9
	?	NC_000913	3377929 =	34 (0.950)					coding (279/393 nt)	rpsI	30S ribosomal subunit protein S9
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGATGTACAGGTCCAGTTTCTCAACCATGTCGAC  >  NC_000913/3377794‑3378045                                                                                                                                    |                                                                                                                         ttCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGAc                                                                                                                   <  1:402325/139‑1 (MQ=255) ttCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGAc                                                                                                                   <  1:54855/139‑1 (MQ=255)                 gCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTa                                                                                                           <  1:227391/131‑1 (MQ=255)                 gCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTa                                                                                                           >  2:227391/1‑131 (MQ=255)                    aaTTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCagag                                                                                                >  1:661004/1‑139 (MQ=255)                               ggAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATa                                                                                     >  1:462084/1‑139 (MQ=255)                                        cGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAgcgc                                                                                             <  2:482338/122‑1 (MQ=255)                                        cGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAgcgc                                                                                             >  1:482338/1‑122 (MQ=255)                                        cGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACgg                                                                            >  2:143723/1‑139 (MQ=255)                                        cGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACgg                                                                            >  1:366629/1‑139 (MQ=255)                                        cGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAAC‑‑‑GGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATc                                                                         >  2:501896/1‑139 (MQ=255)                                         ggACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGa                                                                           <  1:143723/139‑1 (MQ=255)                                                gacgTGCTTCACGCAGACCGACTTTCTTACGTTCTACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAAC‑‑‑GGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGc                                                                 <  1:501896/139‑1 (MQ=255)                                                     gCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCt                                                               <  1:290710/139‑1 (MQ=255)                                                         tACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGAcc                                                           <  2:661004/139‑1 (MQ=255)                                                                aCCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATa                                                    >  1:406807/1‑139 (MQ=255)                                                                        tCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATaccaccacc                                            >  2:334257/1‑139 (MQ=255)                                                                               ttCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACg                                     <  2:366629/139‑1 (MQ=255)                                                                                  aaCCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAggg                                                                                                                     <  2:137869/56‑1 (MQ=255)                                                                                  aaCCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAggg                                                                                                                     >  1:637543/1‑56 (MQ=255)                                                                                  aaCCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAggg                                                                                                                     >  1:137869/1‑56 (MQ=255)                                                                                  aaCCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAggg                                                                                                                     <  2:637543/56‑1 (MQ=255)                                                                                  aaCCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACtcgt                                                                                                               <  1:580763/62‑1 (MQ=255)                                                                                  aaCCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACtcgt                                                                                                               >  2:580763/1‑62 (MQ=255)                                                                                  aaCCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGA‑‑‑‑‑‑‑‑ACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGATGTACAg                          >  1:496560/1‑139 (MQ=255)                                                                                            gCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACgg                                                                            >  1:334980/1‑87 (MQ=255)                                                                                            gCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACgg                                                                            <  2:334980/87‑1 (MQ=255)                                                                                                        aCGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGATGTACAGGTCCAGTTTCTCaa            <  2:668107/139‑1 (MQ=255)                                                                                                         cGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCa                                                                      >  2:718974/1‑80 (MQ=255)                                                                                                         cGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCa                                                                      <  1:718974/80‑1 (MQ=255)                                                                                                            aGCCAGCTTTACGCAGTTGAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACgg                                                                            >  2:436104/1‑71 (MQ=255)                                                                                                            aGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACgg                                                                            <  1:436104/71‑1 (MQ=255)                                                                                                                  cTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGATGTACAGGTCCAGTTTCTCAACCATGTCGAc  <  1:334257/139‑1 (MQ=255)                                                                                                                            ttCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGATGTACAGGTCCa                     <  2:669090/110‑1 (MQ=255)                                                                                                                            ttCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGATGTACAGGTCCa                     >  1:669090/1‑110 (MQ=255)                                                                                                                                    |                                                                                                                         TTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGATGTACAGGTCCAGTTTCTCAACCATGTCGAC  >  NC_000913/3377794‑3378045

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