New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AM260480	2577612 =	53 (1.800)	6 (0.230)	5/192	1.8	13.4%	coding (827/1665 nt)	nirS	cytochrome cd1 nitrite reductase (NirS)
	?	AM260480	2577663 =	30 (1.130)					coding (878/1665 nt)	nirS	cytochrome cd1 nitrite reductase (NirS)
Rejected: Frequency below/above cutoff threshold.

TACGCCGTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATGGTCAACTACTCGGACTTGTCCAACCTGAAGACCACCACCATCGATTCGGCCAAGTTCCTG  >  AM260480/2577533‑2577793                                                                                                                                   |                                                                                                                                   tACGCCGTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTggcggc                                                                                                                                 >  1:391273/1‑134 (MQ=255) tACGCCGTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTggcggc                                                                                                                                 >  2:472000/1‑134 (MQ=255)       gTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCAtcgt                                                                                                                           >  1:280853/1‑134 (MQ=255)                        gccgccGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATc                                                                                                          >  2:556160/1‑134 (MQ=255)                             cgcAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCAtcg                                                                                                                            >  1:821706/1‑111 (MQ=255)                             cgcAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCAtcg                                                                                                                            <  2:821706/111‑1 (MQ=255)                              cttggcgttgatcacgaactccggatggaaatggctggcgacgatagccgCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGTCGCCCATCGTCACCAGCCATTTCCATCCggaat                                                                                                    <  1:784315‑M1/84‑3 (MQ=255)                                 gTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCg                                                                                                 <  1:806138/134‑1 (MQ=255)                                 gTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCg                                                                                                 <  2:391273/134‑1 (MQ=255)                                 gTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCg                                                                                                 <  2:387662/134‑1 (MQ=255)                                 gTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCg                                                                                                 <  1:526983/134‑1 (MQ=255)                                   cgttgatcacgaactccggatggaaatggctggcgacgatggccgccacgcgcggctcggggtggtattcgttgtccacggtcatgccgcgcgtggCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTg                                                                                               <  1:551996‑M2/38‑1 (MQ=255)                                          cacgaactccggatggaaatggctggcgacgatggccgCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACg                                                                                        <  1:551996‑M1/96‑1 (MQ=255)                                           aTGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACg                                                                                        <  1:674348/133‑1 (MQ=255)                                               attccggatggaaatggctggtgacgatgggcgacacgcgcggctcggggtggtattcgttgtccacggtcatgccgcgcgtggCGGCTATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAg                                                                                   <  1:784315‑M2/50‑1 (MQ=255)                                                     acacGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGAcc                                                                             <  2:280853/134‑1 (MQ=255)                                                       acGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGAcc                                                                             <  1:930982/132‑1 (MQ=255)                                                       acGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGAcc                                                                             >  2:930982/1‑132 (MQ=255)                                                       acGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCgg                                                                           <  2:135402/134‑1 (MQ=255)                                                                          tggtCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATGGTCAac                                                        <  1:472000/134‑1 (MQ=255)                                                                           ggtCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATGGTCAac                                                        >  1:909325/1‑133 (MQ=255)                                                                                 CACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCg                                                                                                 <  2:493141‑M1/86‑1 (MQ=255)                                                                                 CACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCg                                                                                                 >  1:493141‑M1/1‑86 (MQ=255)                                                                                 CACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATGGTCAACTACTCg                                                  <  2:919800‑M1/133‑1 (MQ=255)                                                                                   CGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATGGTCAACTACTCGGAc                                               >  1:948884‑M1/1‑134 (MQ=255)                                                                                      GCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATGGTCAACTACTCGGACTTg                                            <  2:282037‑M1/134‑1 (MQ=255)                                                                                        GGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATGGTCAACTACTCGGACTTGTc                                          <  1:287594‑M1/134‑1 (MQ=255)                                                                                           tcggggtggtattcgttgtccacggtcatgccgcgcgtggCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATgg                                                             >  2:909964‑M2/41‑112 (MQ=255)                                                                                           tcggggtggtattcgttgtccacggtcatgccgcgcgtggCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATgg                                                             <  1:909964‑M2/72‑1 (MQ=255)                                                                                           tcggggtggtattcgttgtccacggtcatgccgcgcgtggCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATGGTCAACTACTCGGACTTGTCCaa                                       >  2:699486‑M2/41‑134 (MQ=255)                                                                                           aTGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATGGTCAACTACTCGGACTTGTCCaa                                       <  2:909325/134‑1 (MQ=255)                                                                                           aTGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATGGTCAACTACTCGGACTTGTCCaa                                       >  2:66986/1‑134 (MQ=255)                                                                                              aCCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATGGTCAACTACTCGGACTTGTCCAACCt                                    <  2:784315/134‑1 (MQ=255)                                                                                               ccGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATGGTCAACTACTCGGACt                                              >  2:973140/1‑123 (MQ=255)                                                                                               ccGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATGGTCAACTACTCGGACt                                              <  1:973140/123‑1 (MQ=255)                                                                                                      cAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATgg                                                             >  1:604517/1‑101 (MQ=255)                                                                                                      cAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATgg                                                             <  2:604517/101‑1 (MQ=255)                                                                                                      cAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATGGTCAACTACTCGGACTTGTCCAACCTGAAGacca                            >  2:1051138/1‑134 (MQ=255)                                                                                                        aCGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATGGTCAACTACTCGGACTTGTCCAACCTGAAGaccac                           >  2:763852/1‑133 (MQ=255)                                                                                                         gttgtccacggtcatgccgcgcgtggCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATGGTCAACTACTCGGACTTGTCCAACCTGAAGaccacca                         <  1:699486‑M2/108‑1 (MQ=255)                                                                                                             tACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATGGTCAactact                                                    <  2:389265/103‑1 (MQ=255)                                                                                                             tACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATGGTCAactact                                                    >  1:389265/1‑103 (MQ=255)                                                                                                             tACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATGGTCAACTACTCGGACTTGTCCAACCTGAAGACCACCACCAt                     >  1:677492/1‑134 (MQ=255)                                                                                                                    ccGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGAACGGCAAGATCCTGATGGTCAACTACTCGGACTTGTCCAACCTGAAGACCACCACCATCGATTCg              <  1:556160/134‑1 (MQ=255)                                                                                                                                tggCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATGGTCAACTACTCGGACTTGTCCAACCTGAAGACCACCACCATCGATTCGGCCAAGTTCCTg  >  1:714410‑M2/4‑134 (MQ=255)                                                                                                                                 ggCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATGGTCAACTACTCGGACTTGTCCAAc                                      >  2:588545‑M2/3‑97 (MQ=255)                                                                                                                                 ggCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATGGTCAACTACTCGGACTTGTCCAAc                                      <  1:588545‑M2/95‑1 (MQ=255)                                                                                                                                   |                                                                                                                                   TACGCCGTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATGGTCAACTACTCGGACTTGTCCAACCTGAAGACCACCACCATCGATTCGGCCAAGTTCCTG  >  AM260480/2577533‑2577793

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