New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4422633 =	25 (0.740)	3 (0.090)	3/586	4.2	9.9%	coding (448/651 nt)	ulaD	3‑keto‑L‑gulonate 6‑phosphate decarboxylase
	?	NC_000913	4422782 =	30 (0.910)					coding (597/651 nt)	ulaD	3‑keto‑L‑gulonate 6‑phosphate decarboxylase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below cutoff threshold.

CCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAACTGACCGGTTACTGGACCTGGGAACAGGCGCAACAGTGGCGCGATGCAGGCATTGGGCAGGTGGTTTATCACCGCAGCCGTGACGCGCAGGCCGCAGGCGTGGCGTGGGGCGAAGCGGACATCACCGCGATCAAACGTCTTTCCGATATGGGCTTCAAAGTCACCGTCACCGGAGGCCTGGCGCTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGC  >  NC_000913/4422373‑4422928                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                        ccAAAATTTGCAATGCAGGCAAAAGCCCTTGGCGTATGTGGTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTGCTGTGCGGAGATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAACTGACCGGGTACTGGACCTGGGAACAGGCGCAACAGTGGCGCGATGCAGGCATTGGGCAGGTGGTTTATCACCGCAGCCGTGACGCGCAGGCCGCAGGCGTGGCGTGGGGCGAAGCGGACATCACCGCGATCAAACGTCTTTCCGATATggg                                                                                                                                                                                                                                                                 <  1:154051/301‑1 (MQ=255)       ttGACGATTCAGGCACGATCCTTTCGGGGATGAGCTGGGAAGCCAACGCTGACTGGGTGACGGTAATTTGCTGGGCGGATATCAACACCGCCAAAAACGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAACTGACCAGTTACTGGACCTGGGAACGGGCGCAACAGTGGCGCGATGCAGGCATTGGGCAGGTGGTTTATCACCGCAGCCGTGACGCGCAGGCCGCAGGCGTGGCGTGGGGCGAAGCGGACATCACCGCGATCAAACGTCTTTCCGATATGGGCTTCaa                                                                                                                                                                                                                                                           <  1:147223/301‑1 (MQ=255)             aTGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAACTGACCGGTTACTGGACCTGGGAACAGGCGCAACAGTGGCGCGATGCAGGCATTGGGCAGGTGGTTTATCACCGCAGCCGTGACGCGCAGGCCGCAGGCGTGGCGTGGGGCGAAGCGGACATCACCGCGATCAAACGTCCTTCCGATATGGGCTTCAAAgtca                                                                                                                                                                                                                                                      >  2:91393/1‑300 (MQ=255)                            tttCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTGCGTTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAACTGACCGGTTACTGGACCTGGGAACAGGCGCAACAGTGGCGCGATGCAGGCATTGGGCAGGTGGTTTATCACCGCAGCCGTGACGCGCAGGCCGCAGGCGTGGCGTGGGGCGAAGCGGACATCACCGCGATCAAACGTCTTTCCGATATGGGCTTCAAAGTCACCGTCACCGGAGGCCt                                                                                                                                                                                                                                      <  1:91390/301‑1 (MQ=255)                               cacggtcgtgCCTAGAAGGCAACGAGGGCTTGGGGACGGTAATGTGCTGGGGGGATATCAACACCGGCAAAGGCGCGCCGGACGGGGCAAAAGAGTGTAACGGCGACGTGCAGATCGAACTGACCGGTTACTGGACCTGGGAACAGGCGCAACAGTGGCGCGATGCAGGCATTGGGCAGGTGGTTTATCACCGCAGCCGTGACGCGCAGGCCGCAGGCGTGGCGTGGGGCGAAGCGGACATCACCGCGATCAAACGTCTTTCCGATATGGGCTTCAAAGTCACCGTCACCGGAGGCCTGgc                                                                                                                                                                                                                                   <  2:86317/294‑1 (MQ=255)                                                                          ttGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAACTGACCGGTTACTGGACCTGGGAACAGGCGCAACAGTGGCGCGATGCAGGCATTGGGCAGGTGGTTTATCACCGCAGCCGTGACGCGCAGGCCGCAGGCGTGGCGTGGGGCGAAGCGGACATCACCGCGATCAAACGTCTTTCCGATATGGGCTTCAAAGTCACCGTCACCGGAGGCCTGGCGCTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTc                                                                                                                                                                                        <  2:305818/301‑1 (MQ=255)                                                                                         cAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAACTGACCGGTTACTGGACCTGGGAACAGGCGCAACAGTGGCGCGATGCAGGCATTGGGCAGGTGGTTTATCACCGCAGCCGTGACGCGCAGGCCGCAGGCGTGGCGTGGGGCGACGCGGACATCACCGCGATCAAACGTCTTTCCGATATGGGCTTCAAAGTCACCGCCCCCGGAGGCCTGGCGCCGGAAAGTCTGCCGCTGTCCAAGGGTTTTCCGTTTCGCGCTTTTAGCGCTGGCCgt                                                                                                                                                                         >  1:41867/1‑301 (MQ=255)                                                                                              ccGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAACTGACCGTTTACTGGACCTGGGAACAGGCGCAACAGTGGAGCGATGCAAGCATTGGGCAGGTGGTTCATCACCGCAACCGTGACGCGCAGGCCGCAGGCGTGGCGCGGGGCGAAGCGGACCACACCGCGATCAAACGTTTTTCCGATATGGGTGTTAAAGTCACCGCCACCGGCGGTCCTGCCCTGGAAGATCTTCCCCTGTTTACGGGTGTACCGATTTACCTCTCTATCCCGTTCCgtagt                                                                                                                                                                      >  1:282612/1‑299 (MQ=255)                                                                                              ccGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAACTGACCGGTTACTGGACCTGGGAACAGGCGCAACAGTGGCGCGATGCAGGCATTGGGCAGGTGGTTTATCACCGCAGCCGTGACGCGCAGGCCGCAGGCGTGGCGTGGGGCGAAGCGGACATCACCGCGATCAAACGTCTTTCCGATATGGGCTTCAAAGTCACCGTCACCGGAGGCCTGGCGCTGGAAGATCCGCCGCAGTTCAAGGGTATTCCGATTCAGTGCTTTATCGCGGCCCgtctctt                                                                                                                                                                    >  1:96727/1‑296 (MQ=255)                                                                                                                         aatagTTTAACGGCGACGTGCAGATCGAAGTGACCGGTTACTGGACCTGGGGACAGGCGCAACAGTGGCGCGATGCAGGCATTGTGCAGGTGGTTTATCACCGCAGCCGTGACGCGCAGGCCGCAGGCGTGGCGTGGGGAGAAGCGGACATCACCGCGATCAAACGTCTTTCCGATATGGGCTTCAAAGTCACCGTCACCGGAGGCCTGGCGCTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGTCGCTTCTCCGGTGGAAGc                                                                                                                                         <  2:235859/298‑1 (MQ=255)                                                                                                                            agTTTAACGGCGACGTGCAGATCGAACTGACCGGTTACTGGACCTGGGAACAGGCGCAACAGTGGCGCGATGCAGGCATTGGGCAGGTGGTTTATCACCGCAGCCGTGACGCGCAGGCCGCAGGCGTGGCGTGGGGCGAAGCGGACATCACCGCGATCAAACGTCTTTCCGCTATGGGCTTCAAAGTCACCGTCACCGGAGGCCTGGCGCTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTACCCGTGATGCCCCTTCTCCGGTGGAAGCCg                                                                                                                                       >  2:152947/1‑300 (MQ=255)                                                                                                                                            gCAGATCGAACTGACCGGTTACTGGACCTGGGAACAGGCGCAACAGTGGCGCGATGCAGGCATTGGGCAGGTGGTTTATCACCGCAGCCGTGACGCGCAGGCCGCAGGCGTGGCGTGGGGCGAAGCGGACATCACCGCGATCAAACGTCTTTCCGATATGGGCTTCAAAGTCACCGTCACCGGAGGCCTGGCGCTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGt                                                                                                                      >  2:159449/1‑301 (MQ=255)                                                                                                                                                aTCGAACTGACCGGTTACTGGACCGGGGAACGTGCGCAACAGTGGCACGATGCAGGCATTGGGCAGGTTGTTTATCACCGCAGCCGTAACGCGCAGGCCGCAGGCGTGGCGTGGGGCGAAGCGGACATCACCGCGATCAAACGTCTTTCCGATATGGGCTGCAAAGTCACCGTCACCGGAGGCCTGGCGCTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTc                                                                                                                    <  2:301581/299‑1 (MQ=255)                                                                                                                                                         aCCGGTTACTGGACCTGGGAACAGGCGCAACAGTGGCGCGATGCAGGCATTGGGCAGGTGGTTTATCACCGCAGCCGTGACGCGCAGGCCGCAGGCGTGGCGTGGGGCGAAGCGGACATCACCGCGATCAAACGTCTTTCCGATATGGGCTTCAAAGTCACCGTCACCGGAGGCCTGGCGCTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAAc                                                                                                         <  2:95837/301‑1 (MQ=255)                                                                                                                                                              ttACTGGACCTGGGAACAGGCGCAACAGTGGCGCGATGCAGGCAGTGGGCAGGTGGTTTATCACCGCAGCCGTGACGCGCAGGCGGCAGGCGTGGCGTGGGGCTAAGCGGACATCACCGCGATCAAACGTCTTTCCGATATGGGCTTCAAAGTCACCGTCACCGGAGGCCTGGCGCTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTgg                                                                                                    <  2:145284/301‑1 (MQ=255)                                                                                                                                                               tACTGGACCTGGGAACAGGCGCAACAGTGGCGCGATGCAGGCATTGGGCAGGTGGTTTATCACCGCAGCCGTGACGCGCAGGCCGCAGGCGTGGCGTGGGGCGAAGCGGACATCACCGCGATCAAACGTCTTTCCGATATGGGCTTCAAAGTCACCGTCACCGGAGGCCTGGCGCTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTggg                                                                                                   <  2:303305/301‑1 (MQ=255)                                                                                                                                                                      ccTGGGAACAGGCGCAACAGTGGCGCGATGCAGGCATTGGGCAGGTGGTTTATCACCGCAGCCGTGACGCGCAGGCCGCAGGCGTGGCGTGGGGTGAAGCGGACATCACCGCGATCAAACGTCTTTCCGATATGGGCTTCAAAGTCACCGTCACCGGAGGCCTGGCGCTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAgg                                                                                            >  2:100866/1‑301 (MQ=255)                                                                                                                                                                      ccTGGGAACAGGCGCAACAGTGGCGCGATGCAGGCATTGGGCAGGTGGTTTATCACCGCAGCCGTGACGCGCAGGCCGCAGGCGTGGCGTGGGGCGAAGCGGACATCACCGCGATCAAACGTCTTTCCGATATGGGCTTCAAAGTCACCGTCACCGGAGGCCTGGCGCTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAgg                                                                                            <  2:41868/301‑1 (MQ=255)                                                                                                                                                                                 gcgcAACAGTGGCGCGATGCAGGCATTGGGCAGGTGGTTTATCACGGCAGCCGTGACGCGCAGGCCGCAGGCGTGGCGTGGGGCGAAGCGGACATCACCACGATCAAACGTCTTTACGATATGGGCTTCAAAGTCACCGTCACCGGAGGCCTGGCGCTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGata                                                                                    <  2:151494/298‑1 (MQ=255)                                                                                                                                                                                                     aGGCATTGGGCAGGGAGTTTATCAACGAAGCCGGGACGCGCAGGCCGCAGGAGTGGGGTGGGGCGAAGCGGACATCACCGCGATCAAACGTCTGTCCGATATGGGCTTCAAAGTCACCGTCACCGGAGGCCTGGCGCTGGAAGAGCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCt                                                              <  2:286252/300‑1 (MQ=255)                                                                                                                                                                                                      tgcattggtgagtgatgttATCCCAGCAGCCGTGGCGAGAAGGCAGCGGGCGTGGCGTGGGGCTAAACGGACAGCCGCGAGATCAAACGTCTTTCCGAAATGGGGTTCATAGGCGCCGGTGCCGGTGGCCAGGCGCTAGAAGATCTGCCGCTGTTCAAGTTTTTTCCGATTCACGTCTTGATCGCGGGCCGTAGTTTGCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCAACGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCGAACAAATCCCGCtt                                                             <  1:152944/283‑1 (MQ=255)                                                                                                                                                                                                                                    ccgatctcgcaggccgcaggcGTGGCGGGGGGCGAAGCGGACATCACCGCGATCAAACGTCTTTCCGATATGGGCTTCAAAGTCACCGTCACCGGAGGCCTGGCGCTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCgggg                              <  2:326347/294‑1 (MQ=255)                                                                                                                                                                                                                                           cgcaggccgcaggcGTGGCGTGGGGCGAAGCGGACATCACCGCGATCAAACGTCTTTCCGATATGGGCTTCAAAGTCACCGTCACCGGAGGCCTGGCGCTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATCTGTTGTCCAAACAAATCCCGCTTGGCCTCTTATCAACAGCCCTCCCCCCCGTGGAgctccg                       >  1:326338/1‑296 (MQ=255)                                                                                                                                                                                                                                                          cggggtggGGCGAAGCGGACAGCCCCGCGAGCAAACGTCTGTCCGATATGGGCTTCCAAGGCACCGTCACCGGAGGCCTGGCGCTGGAAGATCTGCCGCTGGTCAAGGGTATTCCGATTCACGTCTTTATCGGGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGc         <  1:217867/296‑1 (MQ=255)                                                                                                                                                                                                                                                                 gggCGAAGCGGACATCCCCGCGATCAAACGTCTTTCCGATATGGGCTTCAAAGTCACCGTCACCGGAGGCCTGGCGCTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGc  <  1:100863/300‑1 (MQ=255)                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                        CCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAACTGACCGGTTACTGGACCTGGGAACAGGCGCAACAGTGGCGCGATGCAGGCATTGGGCAGGTGGTTTATCACCGCAGCCGTGACGCGCAGGCCGCAGGCGTGGCGTGGGGCGAAGCGGACATCACCGCGATCAAACGTCTTTCCGATATGGGCTTCAAAGTCACCGTCACCGGAGGCCTGGCGCTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGC  >  NC_000913/4422373‑4422928

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