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.

AGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGT  >  NC_000913/3587289‑3587576                                                                                                                                              |                                                                                                                                                   agcagcATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAgc                                                                                                                                             >  3:37372/1‑149 (MQ=255)  gcagcaTATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAgcg                                                                                                                                            <  5:143989/149‑1 (MQ=255)    agcaTATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCgca                                                                                                                                          <  7:106585/149‑1 (MQ=255)      caTATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCgcag                                                                                                                                         >  5:67693/1‑148 (MQ=255)                            aCGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGtt                                                                                                                  >  1:19027/1‑149 (MQ=255)                            aCGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGtt                                                                                                                  >  5:86246/1‑149 (MQ=255)                            aCGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGtt                                                                                                                  <  6:90057/149‑1 (MQ=255)                                cGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTAtt                                                                                                              <  6:67693/149‑1 (MQ=255)                                              aGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAAcc                                                                                                 <  2:149602/148‑1 (MQ=255)                                              aGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCa                                                                                                <  2:19027/149‑1 (MQ=255)                                               gtctctattcatctcaatcataagttactcaataaagcgcgagtgatgcagttgaaagacgccggactgcggattctggtttataccgtcaataaACCCCAGCGCAGCAACTCTGCTGCGCGCTggg                                                                                                                     <  2:60077‑M1/32‑1 (MQ=255)                                               gtctctattcatctcaatcataagttactcaataaagcgcgagtgatgcagttgaaagacgccggactgcggattctggtttataccgtcaataaACCCCAGCGCAGCAACTCTGCTGCGCGCTggg                                                                                                                     >  1:60077‑M1/96‑127 (MQ=255)                                                         atctcaatcataagttactctatatagcgcgagtgatgcagttgaaagacgccggacttcggattctggtttataccgtctataaACCCCAGCGCAGCAACTCTGGTGCGCGCTGGGGt                                                                                                                   >  2:122500‑M1/86‑119 (MQ=255)                                                         atctcaatcataagttacacaataaagcgcgagagatgcagtagaaagacgccggaatgcggattctggtttataccgtcaataaACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGt                                                                                                                   <  1:122500‑M1/34‑1 (MQ=255)                                                              aatcataagttactcaataaagcgcgagtgatgcagttgaaagacgccggactgcggattctggtttataccgtcaataaACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTAt                                                                                                               <  5:97034‑M1/38‑1 (MQ=255)                                                              aatcataagttactcaataaagcgcgagtgatgcagttgaaagacgccggactgcggattctggtttataccgtcaataaACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTAt                                                                                                               >  6:97034‑M1/81‑118 (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                                                                                                                 <  7:40726/71‑1 (MQ=255)                                                                                                           aTCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGttt                                                                                                                 >  8:40726/1‑71 (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)                                                                                                                                            acacCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTAtt                                                  >  6:49819/1‑101 (MQ=255)                                                                                                                                            acacCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTAtt                                                  <  5:49819/101‑1 (MQ=255)                                                                                                                                            acacCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGt  >  4:68306/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                   AGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGT  >  NC_000913/3587289‑3587576

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.