New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2925373	32 (0.970)	6 (0.190)	6/232	3.0	17.3%	coding (26/849 nt)	queF	7‑cyano‑7‑deazaguanine reductase (NADPH‑dependent)
	?	NC_000913	= 2925381	27 (0.860)					coding (34/849 nt)	queF	7‑cyano‑7‑deazaguanine reductase (NADPH‑dependent)
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GTAGCGTGCAGTAAAATCTTTCAGGGCTTGTGCGGTCAAATCGTCCACTAATTTCTCTCTTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGATAACCTGCCTTTTCACGGTACGGATATCTGG  >  NC_000913/2925233‑2925518                                                                                                                                             |                                                                                                                                                  gTAGCGTGCAGTAAAATCTTTCAGGGCTTGTGCGGTCAAATCGTCCACTAATTTCTCTCTTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCgg                                                                                                                                           <  5:108112/149‑1 (MQ=255)     cGTGCAGTAAAATCTTTCAGGGCTTGTGCGGTCAAATCGTCCACTAATTTCTCTCTTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCt                                                                                                                                       <  3:75269/149‑1 (MQ=255)               aaTCTTTCAGGGCTTGTGCGGTCAAATCGTCCACTAATTTCTCTCTTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCt                                                                                                                                       <  8:5380/139‑1 (MQ=255)               aaTCTTTCAGGGCTTGTGCGGTCAAATCGTCCACTAATTTCTCTCTTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTATTCTTATGCAAACCATCAGGCACTTGCGGTCCt                                                                                                                                       >  7:5380/1‑139 (MQ=255)                     tCAGGGCTTGTGCGGTCAAATCGTCCACTAATTTCTCTCTTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATc                                                                                                                        <  8:95062/148‑1 (MQ=255)                       aGGGCTTGTGCGGTCAAATCGTCCACTAATTTCTCTCTTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAAc                                                                                                                     <  1:11958/149‑1 (MQ=255)                                    tCAAATCGTCCACTAATTTCTCTCTTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGAt                                                                                                        <  4:28677/149‑1 (MQ=255)                                      aaaTCGTCCACTAATTTCTCTCTTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATAc                                                                                                      >  4:90736/1‑149 (MQ=255)                                            tCCACTAATTTCTCTCTTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTg                                                                                                                                              <  8:51380/103‑1 (MQ=255)                                            tCCACTAATTTCTCTCTTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTg                                                                                                                                              >  7:51380/1‑103 (MQ=255)                                                       ctctctTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGa                                                                                                         >  7:92703/1‑129 (MQ=255)                                                       ctctctTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGa                                                                                                         <  8:92703/129‑1 (MQ=255)                                                       ctctctTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTa                                                                                      <  3:90736/148‑1 (MQ=255)                                                       ctctctTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTAc                                                                                     >  2:51830/1‑149 (MQ=255)                                                                                                    tttATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCaa                                                                                <  3:87259/109‑1 (MQ=255)                                                                                                    tttATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCaa                                                                                >  4:87259/1‑109 (MQ=255)                                                                                                     ttATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGa                                        >  8:121609/1‑148 (MQ=255)                                                                                                                      gTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTcc                                                                        >  4:98279/1‑99 (MQ=255)                                                                                                                         ttcttATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGATAACCTGCCtttt                    >  5:20378/1‑148 (MQ=255)                                                                                                                           cttATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGATAACCTGCCTTTTCa                  <  7:121609/148‑1 (MQ=255)                                                                                                                             tATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGATAACCTGCCTTTTCACgg               >  8:79300/1‑149 (MQ=255)                                                                                                                              aTGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCg                                                                            <  1:94422/87‑1 (MQ=255)                                                                                                                              aTGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCg                                                                            >  2:94422/1‑87 (MQ=255)                                                                                                                               tGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAg                                                                               <  7:109239/83‑1 (MQ=255)                                                                                                                               tGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAg                                                                               >  8:109239/1‑83 (MQ=255)                                                                                                                                gCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGATAACCTGCCTTTTCACGGTAc            >  2:162014/1‑149 (MQ=255)                                                                                                                                    aCCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCAc                                                                      >  1:178512/1‑87 (MQ=255)                                                                                                                                    aCCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCAc                                                                      <  2:178512/87‑1 (MQ=255)                                                                                                                                       aTCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGATAACCTGCCTTTTCACGGTACGGatat      <  6:128254/148‑1 (MQ=255)                                                                                                                                       aTCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGATAACCTGCCTTTTCACGGTACGGatat      >  5:128254/1‑148 (MQ=255)                                                                                                                                        tCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGATAACCTGCCTTTTCACGGTACGGATATCt    >  8:68824/1‑149 (MQ=255)                                                                                                                                          aGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGATAACCTGCCTTTTCACGGTACGGATATCTgg  <  1:51830/149‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                  GTAGCGTGCAGTAAAATCTTTCAGGGCTTGTGCGGTCAAATCGTCCACTAATTTCTCTCTTCACTTAAACCAGATACACTTGTCGTTTAGTTTATCTGGTTTATGACGGTGAAACATGTCTTCTTATGCAAACCATCAGGCACTTGCGGGCCTGACTCTTGGAAAATCAACCGATTACCGGGATACCTATGACGCCAGCCTACTGCAAGGCGTTCCACGCAGCCTGAATCGCGACCCGCTGGGTCTGAAAGCGGATAACCTGCCTTTTCACGGTACGGATATCTGG  >  NC_000913/2925233‑2925518

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