New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1547494 =	11 (0.430)	3 (0.120)	3/320	2.6	16.2%	coding (94/3048 nt)	fdnG	formate dehydrogenase N subunit alpha
	?	NC_000913	= 3828512	20 (0.780)					coding (154/1206 nt)	gltS	glutamate:sodium symporter
Rejected: Frequency below/above cutoff threshold.

GCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAGGGACATATCAAAGTTGACTTCCCAGCCCATGCTTTTTTTCAGTACTAGTAGCGCCAGCGCCACCAACAAACCACCCGCAACAGGTTCCGGTATGGTGTATTTCTTCAAAAAGGAGACGGAATGGACCAACTTACGCCCGAGCAGCAACGTCAGCGTTGCGGCAACAAGCGTTGCTAAAGTATCGAGATGAAACATAGTTACTCCTTTGTATCCGCATCACTCTTTCATACTCGTCATACTTCATCGCGCAGTGGCAACTTCATGAACTATGACGGGTA  >  NC_000913/3828238‑3828748                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                             gCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCACTGCCAATGGCATTTTGCATCACCCACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATCAACGGATCGCGCAGGGACATATCAAAGTTGACTTCCCAGCCCATGCTTTTTTTCAGc                                                                                                                                                                                                                                             <  1:33936‑M2/276‑2 (MQ=255)                                 ccaccccccccccAACAGCACCATCGACGGATCAAACCATAACAAACTAGCACTACCAATTCCAATGGCCGTTTGAATCACCAACATACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGGTCGCGCAGGGACATATCAAAGTTGACTTCCCAGCCCATGCTTTTTTTCAGTACTAGTAGCGCCAGCGCCACCAACAAtccacc                                                                                                                                                                                                             <  2:33936/265‑1 (MQ=255)                                                                                                                    aCAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGAGCAGGGACATATCAAAGTTGACTTCCCAGCCCATGCTTTTTTTCAGTACTAGTAGCGCGAGCGCTACTAACAAACCACCCGCAACAGGTTCCg                                                                                                                                                                                               >  1:445161/1‑207 (MQ=255)                                                                                                                                                     ccACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAGGGACATATCAAAGTTGACTTCCCAGCCCATGCTTTTTTTCAGTACTAGTAGCGCCAGCGCCACCAACAAACCACCCGCAACAGGTTCCGGTATGGTGTATTTCTTCAAAAAgg                                                                                                                                                                       >  2:210902/1‑198 (MQ=255)                                                                                                                                                        cacGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAGGGACATATCAAAGTTGACTTCCCAGCCCATGCTTTTTTTCAGTACTAGTAGCGCCAGCGCCACCAACAAACCACCCGCAACAGGTTCCGGTATGGTGTATTTCTTCAAAAAGgaga                                                                                                                                                                    <  1:210903/198‑1 (MQ=255)                                                                                                                                                                         aaaCTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAGGGACATATCAAAGTTGACTTCCCAGCCCATGCTTTTTTTCAGTACTAGTAGCGCCAGCGCCACCAACAAACCACCCGCAACGGGTTCCGGTATGGTGTATTTCTTCAAAAAGGAGACGGAATGGACCAACTTACGCCCGAGCAGCAACGTCAGCGTTGCGGcaac                                                                                                                   >  1:37545/1‑230 (MQ=255)                                                                                                                                                                           aCTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGACGAAAGCCAGCATTCACGGATCGCGCAGGGACATATCAAAGTTGACTTCCCAGCCCATCCTTTTTTTCAGTACTAGTAGCGCCAGCGCCACCAACAAACCACCCGCAACAgg                                                                                                                                                                                                    >  1:332156/1‑147 (MQ=255)                                                                                                                                                                           aCTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAGGGACATATCAAAGTTGACTTCCCAGCCCATGCTTTTTTTCAGTACTAGTAGCGCCAGCGCCACCAACAAACCACCCGCAAc                                                                                                                                                                                                       >  1:33840/1‑144 (MQ=255)                                                                                                                                                                           aCTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAGGGACATATCAAAGTTGACTTCCCAGCCCATGCTTTTTTTCAGTACTAGTAGCGCCAGCGCCACCAACAAACCACCCGCAACAGGTTCCg                                                                                                                                                                                               >  2:16065/1‑152 (MQ=255)                                                                                                                                                                            cTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAGGGACATATCAAAGTTGACTTCCCAGCCCATGCTTTTTTTCAGTACTAGTAGCGCCAGCGCCACCAACAAACCACCCGCAACGGGTTCCGGTATGGTGTATTTCTTCAAAAAGGAGACGGAATGGACCAACTTACGCCCGAGCAGCAACGTCAGCGTTGCGGCAACAAg                                                                                                                <  2:37545/230‑1 (MQ=255)                                                                                                                                                                              ggCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAGGGACATATCAAAGTTGACTTCCCAGCCCATGCTTTTTTTCAGTACTAGTAGCGCCAGCGCCACCAACAAACCACCCGCAACAgg                                                                                                                                                                                                    <  2:33840/144‑1 (MQ=255)                                                                                                                                                                              ggCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAGGGACATATCAAAGTTGACTTCCCAGCCCATGCTTTTTTTCAGTACTAGTAGCGCCAGCGCCACCAACAAACCACCCGCAACAGGTTCCGGTa                                                                                                                                                                                            <  1:16065/152‑1 (MQ=255)                                                                                                                                                                                                                        cATTAACGGATCGCGCAGGGACATATCAAAGTTGACTTCCCAGCCCATGCTTTTTTTCAGTACTAGTAGCGCCAGCGCCACCAACAAacc                                                                                                                                                                                                                >  2:64769/1‑90 (MQ=255)                                                                                                                                                                                                                           tAACGGATCGCGCAGGGACATATCAAAGTTGACTTCCCAGCCCATGCTTTTTTTCAGTACTAGTAGCGCCAGCGCCACCAACAAaccacc                                                                                                                                                                                                             <  1:64769/90‑1 (MQ=255)                                                                                                                                                                                                                             aCGGATCGCGCAGGGACATATCAAAGTTGACTTCCCAGCCCATGCTTTTTTTCAGTACTAGTAGCGCCAGCGCCACCAACAAACCACCCGCAACAGGTTCCGGTATGGTGTATTTCTTCAAAAAGGAGACGGAATg                                                                                                                                                             >  1:44650/1‑136 (MQ=255)                                                                                                                                                                                                                                gctCGCGCAGGGACATATCAAAGTTGACTTCCCAGCCCATGCTTTTTTTCAGTACTAGTGGCGCCAGCGCCACCAACAAACCACCCGCAACAGGTTCCGGTATGGTGTATTTCTTCAAAAAGGAGACGGAATGGAc                                                                                                                                                          <  2:44650/134‑1 (MQ=255)                                                                                                                                                                                                                                  tCGCGCAGGGACATATCAAAGTTGACTTCCCAGCCCATGCTTTTTTTCAGTACTAGTAGCGCCAGCGCCACCAACAAACCACCCGCAACAGGTTCCGGTATGGTGTATTTCTTCAAAAAGGAGACGGAATGGACCAACTTACGCCCGAGCAGCAACGTCAGCGTTGCGGCAACAAGCGTTGCTAAAGTa                                                                                                   >  2:150453/1‑189 (MQ=255)                                                                                                                                                                                                                                     cgcAGGGACATATCAAAGTTGACTTCCCAGCCCATGCTTTTTTTCAGTACTAGTAGCGCCAGCGCCACCAACAAACCACCCGCAACAGGTTCCGGTATGGTGTATTTCTTCAAAAAGGAGACGGAATGGACCAACTTACGCCCGag                                                                                                                                           >  1:365386/1‑146 (MQ=255)                                                                                                                                                                                                                                     cgcAGGGACATATCAAAGTTGACTTCCCAGCCCATGCTTTTTTTCAGTACTAGTAGCGCCAGCACCACCAACAAACCACCCGCAACAGGTTCCGGTATGGTGTATTTCTTCAAAAAGGAGACGGAATGGACCAACTGACGCCCGAGCAGCAACGTCAGCGTTGCGGCAACAAGCGTTGCTAAAGTATCg                                                                                                <  1:150453/189‑1 (MQ=255)                                                                                                                                                                                                                                        aGGGACATATCAAAGTTGACTTCCCAGCCCATGCTTTTTTTCAGTACTAGTAGCGCCAGCGCCACCAACAAACCACCCGCAACAGGTTCCGGTATGGTGTATTTCTTCAAAAAGGAGACGGAATGGACCAACTTACGCCCGagcag                                                                                                                                        <  2:365382/146‑1 (MQ=255)                                                                                                                                                                                                                                            aCATATCAAAGTTGACTTCCCAGCACATGCTTTTTTTCAGTACTAGTAGCGCCAGCGCCACCAACAAACCACCCGCAACAGGTTCCGGTATGGTGTATTTCTTCAAAAAGGAGACGGAATGGACCAACTTACGCCCGAGCAGCAACGTCAGCGTTGCGGCAACAAGCGTTGCTAAAGGATCGAGATGAAACATAGTTACTCCTTTGTATCCGCATCACTCTTTCATACTCGTCATACTTCATCGCGCAGTGGCAACTTCATGAACTATGACGGGTa  <  1:438523/276‑1 (MQ=255)                                                                                                                                                                                                                                                        tGACTTCCCAGCCCATGCTTTTTTTCAGctggctcaggcgcgaaactacaaattatt                                                                                                                                                                                                                 <  2:427606‑M2/57‑30 (MQ=255)                                                                                                                                                                                                                                                           cTTCCCAGCCCATGCTTTTTTTCAGctggctcaggcgcgaaactacaaattattacg                                                                                                                                                                                                              >  1:427610‑M2/1‑25 (MQ=255)                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                             GCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAGGGACATATCAAAGTTGACTTCCCAGCCCATGCTTTTTTTCAGTACTAGTAGCGCCAGCGCCACCAACAAACCACCCGCAACAGGTTCCGGTATGGTGTATTTCTTCAAAAAGGAGACGGAATGGACCAACTTACGCCCGAGCAGCAACGTCAGCGTTGCGGCAACAAGCGTTGCTAAAGTATCGAGATGAAACATAGTTACTCCTTTGTATCCGCATCACTCTTTCATACTCGTCATACTTCATCGCGCAGTGGCAACTTCATGAACTATGACGGGTA  >  NC_000913/3828238‑3828748

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