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.

CATGGAGGCGCGCTCGGTCGAGACCTCCAAGTACAAGGGCTATGAAGACAAGTACGCCGTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATGGTCAAC  >  AM260480/2577481‑2577739                                                                                                                                    |                                                                                                                                caggaacttggccgaatcgatggtggtggtcttcaggttggacaagtccgagtagttgaccatcaggatcttgccggtctccttggcgttgatcacgaactccggatggaaatggctggcgacgatggccgCCa                                                                                                                               >  1:714410‑M1/132‑134 (MQ=255) cATGGAGGCGCGCTCGGTCGAGACCTCCAAGTACAAGGGCTATGAAGACAAGTACGCCGTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCa                                                                                                                               <  1:720885/134‑1 (MQ=255)    ggAGGCGCGCTCGGTCGAGACCTCCAAGTACAAGGGCTATGAAGACAAGTACGCCGTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCAcgc                                                                                                                            >  2:948884/1‑134 (MQ=255)     gAGGCGCGCTCGGTCGAGACCTCCAAGTACAAGGGCTATGAAGACAAGTACGCCGTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCAcgcg                                                                                                                           >  1:1018397/1‑134 (MQ=255)      aGGCGCGCTCGGTCGAGACCTCCAAGTACAAGGGCTATGAAGACAAGTACGCCGTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCAcgcgc                                                                                                                          <  2:823405/134‑1 (MQ=255)      aGGCGCGCTCGGTCGAGACCTCCAAGTACAAGGGCTATGAAGACAAGTACGCCGTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCAcgcgc                                                                                                                          <  1:656109/134‑1 (MQ=255)      aGGCGCGCTCGGTCGAGACCTCCAAGTACAAGGGCTATGAAGACAAGTACGCCGTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCAcgcgc                                                                                                                          <  1:643764/134‑1 (MQ=255)      aGGCGCGCTCGGTCGAGACCTCCAAGTACAAGGGCTATGAAGACAAGTACGCCGTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCAcgcgc                                                                                                                          <  1:384657/134‑1 (MQ=255)      aGGCGCGCTCGGTCGAGACCTCCAAGTACAAGGGCTATGAAGACAAGTACGCCGTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCAcgcgc                                                                                                                          <  1:29215/134‑1 (MQ=255)                 gTCGAGACCTCCAAGTACAAGGGCTATGAAGACAAGTACGCCGTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGt                                                                                                               >  2:641943/1‑134 (MQ=255)                        tggtggtcttcaggttggacaagtccgagtagttgaccatcaggatcttgccggtctccttggcgttgatcacgaactccggatggaaatggctggcgacgatggccgCCACGCGCGGCATGACCGTGGACAAc                                                                                                        <  1:699486‑M1/26‑1 (MQ=255)                         cTCCAAGTACAAGGGCTATGAAGACAAGTACGCCGTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACg                                                                                                       <  2:1018397/134‑1 (MQ=255)                         cTCCAAGTACAAGGGCTATGAAGACAAGTACGCCGTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACg                                                                                                       <  2:934599/134‑1 (MQ=255)                                     gttggacaagtccgagtagttgaccatcaggatcttgccggtctccttggcgttgatcacgaactccggatggaaatggctggcgacgatggccgcc                                                                                                                                >  2:588545‑M1/96‑97 (MQ=255)                                     gttggacaagtccgagtagttgaccatcaggatcttgccggtctccttggcgttgatcacgaactccggatggaaatggctggcgacgatggccgcc                                                                                                                                <  1:588545‑M1/2‑1 (MQ=255)                                      ttggacaagtccgagtagttgaccatcaggatcttgccggtctccttggcgttgatcacgaactccggatggaaatggctggcgacgatggccgCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGa                                                                                          >  2:699486‑M1/95‑134 (MQ=255)                                       gCTATGAAGACAAGTACGCCGTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCg                                                                                                                >  2:757930/1‑111 (MQ=255)                                       gCTATGAAGACAAGTACGCCGTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCg                                                                                                                <  1:757930/111‑1 (MQ=255)                                        cTATGAAGACAAGTACGCCGTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCAcgcgc                                                                                                                          >  1:944131/1‑100 (MQ=255)                                        cTATGAAGACAAGTACGCCGTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCAcgcgc                                                                                                                          <  2:944131/100‑1 (MQ=255)                                             aaGACAAGTACGCCGTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCcgcgc                                                                                   >  1:387662/1‑134 (MQ=255)                                              aGACAAGTACGCCGTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCg                                                                                                                >  2:565589/1‑104 (MQ=255)                                              aGACAAGTACGCCGTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCg                                                                                                                <  1:565589/104‑1 (MQ=255)                                              aGACAAGTACGCCGTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACa                                                                                                          >  2:450250/1‑110 (MQ=255)                                              aGACAAGTACGCCGTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACa                                                                                                          <  1:450250/110‑1 (MQ=255)                                              aGACAAGTACGCCGTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCcgcgcg                                                                                  >  2:526983/1‑134 (MQ=255)                                              aGACAAGTACGCCGTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCcgcgcg                                                                                  >  2:649030/1‑134 (MQ=255)                                                     tACGCCGTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGtctc                                                                                                                                 >  1:604514/1‑80 (MQ=255)                                                     tACGCCGTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGtctc                                                                                                                                 <  2:604514/80‑1 (MQ=255)                                                     tACGCCGTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCAcgcgc                                                                                                                          <  2:4634/87‑1 (MQ=255)                                                     tACGCCGTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCAcgcgc                                                                                                                          >  1:4634/1‑87 (MQ=255)                                                     tACGCCGTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTggcggc                                                                           >  2:472000/1‑134 (MQ=255)                                                     tACGCCGTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTggcggc                                                                           >  1:391273/1‑134 (MQ=255)                                                           gTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTg                                                                                                              >  1:238151/1‑93 (MQ=255)                                                           gTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTg                                                                                                              <  2:238151/93‑1 (MQ=255)                                                           gTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCAtcgt                                                                     >  1:280853/1‑134 (MQ=255)                                                            ccatcaggatcttgccggtctccttggcgttgatcacgaactccggatggaaatggctggcgacgatggccgCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGa                                                                                          >  2:909964‑M1/73‑112 (MQ=255)                                                            ccatcaggatcttgccggtctccttggcgttgatcacgaactccggatggaaatggctggcgacgatggccgCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGa                                                                                          <  1:909964‑M1/40‑1 (MQ=255)                                                                            gccgccGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACg                                                                                                       >  1:554561/1‑83 (MQ=255)                                                                            gccgccGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACg                                                                                                       <  2:554561/83‑1 (MQ=255)                                                                            gccgccGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATc                                                    >  2:556160/1‑134 (MQ=255)                                                                                 cgcAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCcg                                                                                      >  2:661025/1‑95 (MQ=255)                                                                                 cgcAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCcg                                                                                      <  1:661025/95‑1 (MQ=255)                                                                                 cgcAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCAtcg                                                                      <  2:821706/111‑1 (MQ=255)                                                                                 cgcAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCAtcg                                                                      >  1:821706/1‑111 (MQ=255)                                                                                  cttggcgttgatcacgaactccggatggaaatggctggcgacgatagccgCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGTCGCCCATCGTCACCAGCCATTTCCATCCggaat                                              <  1:784315‑M1/84‑3 (MQ=255)                                                                                     gTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGtctc                                                                                                                                 >  1:27163/1‑48 (MQ=255)                                                                                     gTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGtctc                                                                                                                                 <  2:27163/48‑1 (MQ=255)                                                                                     gTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCg                                           <  1:806138/134‑1 (MQ=255)                                                                                     gTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCg                                           <  2:391273/134‑1 (MQ=255)                                                                                     gTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCg                                           <  1:526983/134‑1 (MQ=255)                                                                                     gTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCg                                           <  2:387662/134‑1 (MQ=255)                                                                                              cacgaactccggatggaaatggctggcgacgatggccgCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACg                                  <  1:551996‑M1/96‑1 (MQ=255)                                                                                               aTGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACg                                  <  1:674348/133‑1 (MQ=255)                                                                                                    aGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGAcc                                                                                                                 >  2:501198/1‑49 (MQ=255)                                                                                                    aGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGAcc                                                                                                                 <  1:501198/49‑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)                                                                                                                                    |                                                                                                                                CATGGAGGCGCGCTCGGTCGAGACCTCCAAGTACAAGGGCTATGAAGACAAGTACGCCGTGGCCGGCTCCTACTGGCCGCCGCAGTACGTGATCATGGAAGGCGACACGCTCAAGCCGCTCAAGGTGGTCTCCACGCGCGGCATGACCGTGGACAACGAATACCACCCCGAGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATGGTCAAC  >  AM260480/2577481‑2577739

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.