New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3372642 =	30 (1.110)	3 (0.120)	3/278	3.7	9.8%	intergenic (‑70/+39)	nanT/nanA	sialic acid transporter/N‑acetylneuraminate lyase
	?	NC_000913	3372677 =	27 (1.060)					intergenic (‑105/+4)	nanT/nanA	sialic acid transporter/N‑acetylneuraminate lyase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

ATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCATTGTGAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACAACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAACGGTTTGCGGCACAGCGGCACAGAAACGACATCCATATAATGGAGGACAGTTTTCAGGCCGCGGAATACG  >  NC_000913/3372534‑3372825                                                                                                                                                |                                                                                                                                                     aTGGCGATACCACGGGATATTCTGGGTTGTAGTACTCATTGTGAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACAACTc                                                                                                                                                 >  2:311459/1‑149 (MQ=255)                 atatTCTGGGTTGTAGTACTCAGTGTGAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACAACTCACCCGCGCTCTTGCAt                                                                                                                                 <  1:121594/149‑1 (MQ=255)                           ttGTAGTACTCATTGTGAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACAACTCACCCGCGCTCTTGCATCAACTGCTgg                                                                                                                       >  2:76065/1‑149 (MQ=255)                                gtaCTCATTGTGAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACAACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAg                                                                                                                  >  1:182393/1‑149 (MQ=255)                                   cTCATTGTGAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACAACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAgcgc                                                                                                               >  1:111055/1‑149 (MQ=255)                                      aTTGTGAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACAACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCtt                                                                                                            >  1:34413/1‑149 (MQ=255)                                             gTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTAAAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGAGAGGGGAAACAACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCCGTTCt                                                                                                     >  1:355384/1‑149 (MQ=255)                                              tctcCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACAACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTg                                                                                                    >  1:160835/1‑149 (MQ=255)                                                        tGGGCGAGGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACAACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATAttt                                                                                          <  1:105071/149‑1 (MQ=255)                                                            cGATGCCCACACGCTTTGGTAAGATTGTTATAGGGTACAGATGCGGTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGCAACAACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGGCTTCCGTTCCGGCAGAGAtttttga                                                                                      >  2:250050/1‑147 (MQ=255)                                                             caTGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACAACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCAt                                                                                     <  2:351371/148‑1 (MQ=255)                                                                 cccACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACAACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTAc                                                                                 <  2:68619/149‑1 (MQ=255)                                                                  ccACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAAGGGAAACAACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTAcc                                                                                >  1:44517/1‑149 (MQ=255)                                                                       cGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACAACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTcc                                                                           <  1:365477/149‑1 (MQ=255)                                                                        gCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACAACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCa                                                                          >  2:78081/1‑149 (MQ=255)                                                                            cgttaggAAAATTGTAGGGTACAGCTGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACAACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTTGCAGATATTTTTCATCTACCGGTCCAAACg                                                                      <  2:2728/143‑1 (MQ=255)                                                                                aTGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACAACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAACGGttt                                                                  >  2:353684/1‑149 (MQ=255)                                                                                  gAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACAACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAACGGTTTgc                                                                >  1:294860/1‑149 (MQ=255)                                                                                   aaaaTTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACAACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAACGGTTTgcg                                                               <  1:184749/149‑1 (MQ=255)                                                                                    aaaTTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACAACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAACGGTTTgcgg                                                              <  1:76065/149‑1 (MQ=255)                                                                                             ggTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACAACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAACCGGTTTgcggcacagcggc                                                     >  1:21492/1‑149 (MQ=255)                                                                                                      gCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACAACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGGCCAAACGGTTTGCGGCACAGCGGCACAGAAACg                                            <  2:107620/149‑1 (MQ=255)                                                                                                                ccccTCACCCGGTAGGGGCGAGCGAGGGGAAACAACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAACGGTTTGCGGCACAGCGGCACAGAAACGACATCCatat                                  >  2:261599/1‑149 (MQ=255)                                                                                                                             agcagcgggcgaGAGGGGAGAAGTCACCGGCGCTCGTGCAGCAACGGCTGGGCCAGCGCATTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAACGGTTTGCGGCACAGCGGCACCGAAACGACATCCATATAATGGAGGACAGt                     >  2:329251/5‑149 (MQ=255)                                                                                                                                ggCGAGCGAGGGGAAACAACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAACGGTTTGCGGCACAGCGGCACAGAAACGACATCCATATAATGGAGGACAGtttt                  <  2:294860/149‑1 (MQ=255)                                                                                                                                ggCGAGCGAGGGGAAACAACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAACGGTTTGCGGCACAGCGGCACAGAAACGACATCCATATAATGGAGGACAGtttt                  <  2:443293/149‑1 (MQ=255)                                                                                                                                      cgagGGGAAACAACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAACGGTTTGCGGCACAGCGGCACAGAAACGACATCCATATAATGGAGGACAGCTTTCAGGcc            >  2:287021/1‑149 (MQ=255)                                                                                                                                                caacTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAACGGTTTGCGGCACAGCGGCACAGAAACGACATCCATATAATGGAGGACAGTTTTCAGGCCGCGGAATACg  <  1:416054/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                     ATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCATTGTGAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACAACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAACGGTTTGCGGCACAGCGGCACAGAAACGACATCCATATAATGGAGGACAGTTTTCAGGCCGCGGAATACG  >  NC_000913/3372534‑3372825

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.