New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2631969	38 (1.160)	6 (0.190)	6/230	3.0	14.3%	coding (567/1578 nt)	guaA	GMP synthetase (glutamine aminotransferase)
	?	NC_000913	= 2631973	36 (1.160)					coding (563/1578 nt)	guaA	GMP synthetase (glutamine aminotransferase)
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CACACCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGCTGGCGTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCT  >  NC_000913/2631822‑2632112                                                                                                                                                    |                                                                                                                                                cacaCCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGCTGGCGTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATAcccg                                                                                                                                                >  5:127559/1‑148 (MQ=255)                     ggATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGCTGGCGTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTcc                                                                                                                            >  6:122526/1‑149 (MQ=255)                      gATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGCTGGCGTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCg                                                                                                                           <  2:23255/149‑1 (MQ=255)                                               tCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGCTGGCGTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCgg                                                                                                                                           >  8:49436/1‑108 (MQ=255)                                               tCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGCTGGCGTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCgg                                                                                                                                           <  7:49436/108‑1 (MQ=255)                                                 gcgGATGCGAGCTACAGCATCGTCGATAATTTTCGCTGGCGTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGtt                                                                                                <  5:122526/149‑1 (MQ=255)                                                         gAGCTACAGCATCGTCGATAATTTTCGCTGGCGTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCg                                                                                        >  1:180458/1‑149 (MQ=255)                                                              aCAGCATCGTCGATAATTTTCGCTGGCGTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGaa                                                                                                     <  7:33008/131‑1 (MQ=255)                                                              aCAGCATCGTCGATAATTTTCGCTGGCGTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGaa                                                                                                     >  8:33008/1‑131 (MQ=255)                                                                       tcgATAATTTTCGCTGGCGTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTa                                                                                     >  5:100826/1‑138 (MQ=255)                                                                       tcgATAATTTTCGCTGGCGTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTa                                                                                     <  6:100826/138‑1 (MQ=255)                                                                          aTAATTTTCGCTGGCGTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCaaa                                                                       <  2:180458/149‑1 (MQ=255)                                                                              ttttCGCTGGCGTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTttctt                                                                                          <  5:164277/126‑1 (MQ=255)                                                                              ttttCGCTGGCGTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTttctt                                                                                          >  6:164277/1‑126 (MQ=255)                                                                               tttCGCTGGCGTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGaa                                                                                                                    <  3:100400/99‑1 (MQ=255)                                                                               tttCGCTGGCGTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGaa                                                                                                                    >  4:100400/1‑99 (MQ=255)                                                                               tttCGCTGGCGTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGc                                                                  >  4:80168/1‑149 (MQ=255)                                                                                            cacaGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGaa                                                                                                     <  2:61761/101‑1 (MQ=255)                                                                                            cacaGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGaa                                                                                                     >  1:61761/1‑101 (MQ=255)                                                                                             acaGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTg                                                    >  3:40081/1‑149 (MQ=255)                                                                                              caGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTgg                                                   >  3:138761/1‑149 (MQ=255)                                                                                                           tGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGtt                                                                                                <  3:89970/91‑1 (MQ=255)                                                                                                           tGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGtt                                                                                                >  4:89970/1‑91 (MQ=255)                                                                                                           tGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGtgatg                                        >  2:89411/1‑147 (MQ=255)                                                                                                           tGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGcgatgtt                                      >  7:72060/1‑147 (MQ=255)                                                                                                             gCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATa                                                                                                        <  4:34651/81‑1 (MQ=255)                                                                                                             gCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATa                                                                                                        >  3:34651/1‑81 (MQ=255)                                                                                                             gCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGcgatgtt                                      <  8:72060/147‑3 (MQ=255)                                                                                                                  tatCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTcggac                               <  4:40081/149‑1 (MQ=255)                                                                                                                   atCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCg                                                          <  6:150306/121‑1 (MQ=255)                                                                                                                   atCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCg                                                          >  5:150306/1‑121 (MQ=255)                                                                                                                          acaAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAAt                                                                             <  2:28489/95‑1 (MQ=255)                                                                                                                          acaAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAAt                                                                             >  1:28489/1‑95 (MQ=255)                                                                                                                               aCGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTaa                  <  3:80168/149‑1 (MQ=255)                                                                                                                                         aTGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACgg                                             <  7:41628/112‑1 (MQ=255)                                                                                                                                         aTGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACgg                                             >  8:41628/1‑112 (MQ=255)                                                                                                                                             gcATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGtt                                                                                                <  4:33034/57‑1 (MQ=255)                                                                                                                                             gcATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGtt                                                                                                >  3:33034/1‑57 (MQ=255)                                                                                                                                               aTACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCt  >  4:102904/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                CACACCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGCTGGCGTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCT  >  NC_000913/2631822‑2632112

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