New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 327138	26 (0.960)	3 (0.110)	3/288	3.8	10.5%	coding (110/1671 nt)	betA	choline dehydrogenase, a flavoprotein
	?	NC_000913	= 327162	26 (0.990)					coding (86/1671 nt)	betA	choline dehydrogenase, a flavoprotein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAG  >  NC_000913/326993‑327283                                                                                                                                                  |                                                                                                                                                  cAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGGCAAAGCGCTAGTCCGGgc                                                                                                                                                >  2:192082/1‑149 (MQ=255)       tttACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGgccgccgg                                                                                                                                          >  1:290581/1‑147 (MQ=255)           ccGCGTCCGCACTCCATGCGGCGGTTCTTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTCCCCTGTAGCGGGAATGCCAGGGCAGCGGGCACCTCGTTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTc                                                                                                                                      <  2:409623/149‑1 (MQ=255)            cgcgTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCa                                                                                                                                     >  2:348061/1‑149 (MQ=255)             gcgTCCGCACTCCATGCGGCGGTTATTCATACACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCaa                                                                                                                                    <  1:367103/149‑1 (MQ=255)             gcgTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGTCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCaa                                                                                                                                    >  1:361085/1‑149 (MQ=255)                                      ttCATAAACGGTTCAGGCTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATc                                                                                                           <  1:56803/149‑1 (MQ=255)                                        cATACACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTACTCggcgctt                                                                                                         >  2:27767/1‑144 (MQ=255)                                         aTAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTtc                                                                                                        >  1:190498/1‑149 (MQ=255)                                          tAATCGGTCCAGGTTGCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTGGCGGGAATGCCAGGGCAGCGGGCATCTGGGGGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTtca                                                                                                       <  1:60453/149‑1 (MQ=255)                                          tAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTtca                                                                                                       <  1:209085/149‑1 (MQ=255)                                                   cAGGTTCCGTTTCCTAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACgg                                                                                              <  1:201600/149‑1 (MQ=255)                                                        tCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGCGGGGttgg                                                                                         <  2:56012/149‑4 (MQ=255)                                                                    gCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGcc                                                                             <  1:348061/149‑1 (MQ=255)                                                                     cccGGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCACGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCg                                                                            <  2:287403/149‑1 (MQ=255)                                                                     cccAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAACGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGGGGTATGCGGAGCTTCAGTCAGACGGGTAGCGAGAACGTTGCCg                                                                            >  2:31545/1‑149 (MQ=255)                                                                      ccAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCgg                                                                           >  1:251986/1‑149 (MQ=255)                                                                               gCGGTTACCCGGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGGGCGGACGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGCCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCgg                                                                  <  2:278636/149‑1 (MQ=255)                                                                                          gTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCCGAGCCGGCACCAATAatg                                                       >  2:20452/1‑149 (MQ=255)                                                                                                                      cggggTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAAcctcct                           <  1:148119/147‑1 (MQ=255)                                                                                                                        gggTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGAAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTgg                         >  2:421990/1‑149 (MQ=255)                                                                                                                              ggAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAACACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGaa                   >  1:44893/1‑149 (MQ=255)                                                                                                                                         cGATGGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGa        <  2:255040/149‑1 (MQ=255)                                                                                                                                         cGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGa        <  2:130333/149‑1 (MQ=255)                                                                                                                                               tCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACAGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGACTATGGGCTGGAATTTAg  >  2:285549/1‑149 (MQ=255)                                                                                                                                               tCCGGGCCGCCCGCTTCAAGCAGCAGCAAGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAg  <  1:83401/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                  CAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAG  >  NC_000913/326993‑327283

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