New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3497094 =	29 (1.110)	3 (0.120)	3/274	2.8	9.7%	coding (92/807 nt)	nirC	nitrite transporter
	?	NC_000913	3497106 =	29 (1.200)					coding (104/807 nt)	nirC	nitrite transporter
Rejected: Frequency below cutoff threshold.

TTTAAACTTTTTGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTGGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTC  >  NC_000913/3496953‑3497242                                                                                                                                              |                                                                                                                                                     cggtgaacagttcagacccggcgataatcaccagcgttaaggcgataccaaaggtcgcgcccatcaccaaagggcgtacggatggatcgagcaaattaccgggcgtgaaaatcaggatgatcccaagccccacatacgcgcCCGCCATggc                                                                                                                                             >  1:169328‑M1/142‑151 (MQ=255)                       ttGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTggg                                                                                                                       <  1:349546/151‑1 (MQ=255)                        tGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGa                                                                                                                      <  1:129058/151‑1 (MQ=255)                          tttttATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGat                                                                                                                     <  2:15276/150‑1 (MQ=255)                              tATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATcc                                                                                                                <  1:397500/151‑1 (MQ=255)                              tATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATcc                                                                                                                <  1:255755/151‑1 (MQ=255)                                    gttaaggcgataccaaaggccgcgcccatcaccaaagggcgtacggatggatcgagcaaattaccgagcgtgaaaatcaggatgatcccaagacccacatactcgcCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGCttt                                                                                                          <  1:378014‑M1/45‑1 (MQ=255)                                           aTAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGc                                                                                                    >  1:83987/1‑150 (MQ=255)                                           aTAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCt                                                                                                   <  2:76268/151‑1 (MQ=255)                                                    tGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAAttt                                                                                          <  2:357104/151‑1 (MQ=255)                                                     gTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGCGGGTCTTGGGATCATCCCGATTTTCACGCCCGGGAATTTg                                                                                         >  1:160964/1‑151 (MQ=255)                                                                cTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCCTCCTGCTTTTCACGCTCGGTAATTTGCTCGatccat                                                                               >  2:10030/1‑150 (MQ=255)                                                                  accaaagggcgtacggatggatcgagcaaattcccgagcgtgaaaatcaggatgatcccaagacccacatacgcgcCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCg                                                                            >  1:335154‑M1/77‑151 (MQ=255)                                                                  aTTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCg                                                                            >  2:73897/1‑151 (MQ=255)                                                                  aTTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCg                                                                            >  1:190584/1‑151 (MQ=255)                                                                  aTTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGGAATTTGCTCGATCCATCCg                                                                            >  1:335888/1‑151 (MQ=255)                                                                  aTTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCCCCGCCATGGCGGGCGCGGATGTGGGTCTTGGGATCATCCTGATTTTCCCGCTCGGTAATTTGCTCGATCCATCCg                                                                            >  2:374114/1‑151 (MQ=255)                                                                  aTTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGGCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCg                                                                            <  2:69952/151‑1 (MQ=255)                                                                  aTTAATAAGGGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGGCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCg                                                                            >  2:254443/1‑151 (MQ=255)                                                                        aaGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGcc                                                                      <  2:312926/151‑1 (MQ=255)                                                                               cGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTGGt                                                               <  1:62730/151‑1 (MQ=255)                                                                                                 ttGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGGAATTTGCTCGATCCCTCCGTACGCCCTTTGGGGATGGGCGCGACCTTTgg                                             >  1:398256/1‑151 (MQ=255)                                                                                                 ttGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGGCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTGGTGATGGGCGCGACCTTTgg                                             >  1:233639/1‑151 (MQ=255)                                                                                                        ccTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGGATGGGGGTCTTGGGATCATCCTGATTTTCACGCCCGGTAATTTGCTCGATCCATCCGTACGCCCTTTTGTGATGGGCGCGGCCTTTGGTATCGc                                       >  2:13332/1‑150 (MQ=255)                                                                                                                cAAATAACCCGATCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTGGTGATTGGCGCGACCTTTGGTATCGCCTTAACGc                               <  2:142064/150‑1 (MQ=255)                                                                                                                     aaCCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTGGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGa                         <  1:197853/151‑1 (MQ=255)                                                                                                                            tCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTGGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGc                  <  1:242970/151‑1 (MQ=255)                                                                                                                               gacccacatacgcgcCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTGGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGGGATTATCGCCgg               <  1:200621‑M1/136‑1 (MQ=255)                                                                                                                               gCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGCTTTTCACGCTCGGTAATTTGCTCGATCCCTCCGTACGCCCTTTGGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGGGATTATCGCCgg               >  2:51054/1‑151 (MQ=255)                                                                                                                               gCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTGGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCgg               >  1:63051/1‑151 (MQ=255)                                                                                                                                        tCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTGGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAAc       <  1:191226/150‑1 (MQ=255)                                                                                                                                            cTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTGGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGt    <  2:124711/149‑1 (MQ=255)                                                                                                                                            cTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTGGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTc  <  2:334808/151‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                     TTTAAACTTTTTGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTGGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTC  >  NC_000913/3496953‑3497242

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