New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3587430 =	28 (0.970)	6 (0.230)	6/232	2.5	18.8%	coding (684/744 nt)	ugpQ	glycerophosphodiester phosphodiesterase, cytosolic
	?	NC_000913	3587464 =	26 (0.980)					coding (650/744 nt)	ugpQ	glycerophosphodiester phosphodiesterase, cytosolic
Rejected: Frequency below/above cutoff threshold.

ACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCA  >  NC_000913/3587316‑3587600                                                                                                                                                     |                                                                                                                                         aCGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGtt                                                                                                                                          >  1:19027/1‑149 (MQ=255) aCGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGtt                                                                                                                                          <  6:90057/149‑1 (MQ=255) aCGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGtt                                                                                                                                          >  5:86246/1‑149 (MQ=255)     cGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTAtt                                                                                                                                      <  6:67693/149‑1 (MQ=255)                   aGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAAcc                                                                                                                         <  2:149602/148‑1 (MQ=255)                   aGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCa                                                                                                                        <  2:19027/149‑1 (MQ=255)                                   aatcataagttactcaataaagcgcgagtgatgcagttgaaagacgccggactgcggattctggtttataccgtcaataaACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTAt                                                                                                                                       >  6:97034‑M1/81‑118 (MQ=255)                                   aatcataagttactcaataaagcgcgagtgatgcagttgaaagacgccggactgcggattctggtttataccgtcaataaACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTAt                                                                                                                                       <  5:97034‑M1/38‑1 (MQ=255)                                           ttCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTc                                                                                                 >  8:56443/1‑148 (MQ=255)                                              cGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACt                                                                                             <  6:86246/149‑1 (MQ=255)                                              cGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACt                                                                                             <  4:37372/149‑1 (MQ=255)                                                   aaaaCTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATc                                                                                        >  4:74079/1‑149 (MQ=255)                                                             ggCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCtt                                                                              >  2:116786/1‑149 (MQ=255)                                                                    aaGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGa                                                                        <  3:74079/148‑1 (MQ=255)                                                                       ttCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAGTaa                                                                    >  2:53903/1‑149 (MQ=255)                                                                                aTCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGttt                                                                                                                                         >  8:40726/1‑71 (MQ=255)                                                                                aTCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGttt                                                                                                                                         <  7:40726/71‑1 (MQ=255)                                                                                aTCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAGTAACTTATGAtt                                                           <  1:53903/149‑1 (MQ=255)                                                                                   aCGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAGTAACTTATGATTgag                                                        <  1:116786/149‑1 (MQ=255)                                                                                       cAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAAc                                                                                              >  5:91883/1‑107 (MQ=255)                                                                                       cAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAAc                                                                                              <  6:91883/107‑1 (MQ=255)                                                                                                            aaTCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCACCTGCATCACTCGCGCTTTATTGCGTAACTTATGATTGAGAt                                                      >  3:80519/1‑126 (MQ=255)                                                                                                            aaTCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAGTAACTTATGATTGAGAt                                                      <  4:80519/126‑1 (MQ=255)                                                                                                               ataaaccccagcgcgcagcagagttgctgcgctggggtTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAGTAACTTATGAtt                                                           <  5:97034‑M2/80‑1 (MQ=255)                                                                                                               ataaaccccagcgcgcagcagagttgctgcgctggggtTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAGTAACTTATGAtt                                                           >  6:97034‑M2/39‑118 (MQ=255)                                                                                                                 acacCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTAtt                                                                          >  6:49819/1‑101 (MQ=255)                                                                                                                 acacCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTAtt                                                                          <  5:49819/101‑1 (MQ=255)                                                                                                                 acacCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGt                          >  4:68306/1‑149 (MQ=255)                                                                                                                   accccagcgcgcagcagagttgctgcgctggggtTTATTGACGGTATAAACCAGAATCCGCATTCCGGCGTCTTTCTACTGCATCTCTCGCGCTTTATTGTGTAACTTATGATTGAGAt                                                      <  1:122500‑M2/85‑1 (MQ=255)                                                                                                                   accccagcgcgcaccagagttgctgcgctggggtTTATAGACGGTATAAACCAGAATCCGAAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTATATAGAGTAACTTATGATTGAGAt                                                      >  2:122500‑M2/35‑119 (MQ=255)                                                                                                                     cccagcgcgcagcagagttgctgcgctggggtTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAGTAACTTATGATTGAGATGAATAGAGAc                                            >  1:60077‑M2/33‑127 (MQ=255)                                                                                                                     cccagcgcgcagcagagttgctgcgctggggtTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAGTAACTTATGATTGAGATGAATAGAGAc                                            <  2:60077‑M2/95‑1 (MQ=255)                                                                                                                     CCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGcccagccgcgcgaactga                      <  3:27410‑M1/149‑7 (MQ=255)                                                                                                                                       ctgcGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCAGCACTCGCGCTTTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCgcgc    <  7:56443/149‑1 (MQ=255)                                                                                                                                          cgcgcTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAGTAACTTATGATGGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCa  <  3:68306/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                         ACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCA  >  NC_000913/3587316‑3587600

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