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.

CGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCATTGTGAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACAACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAACGGTTTGCGGCACAGCGGCACAGAAACG  >  NC_000913/3372503‑3372783                                                                                                                                            |                                                                                                                                              cGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCATTGTGAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCa                                                                                                                                      >  1:432609/1‑149 (MQ=255)        tGCGCGTCATTGTGCACGGTTGAGAGTGCGGTACCACGGGATATTCTGGGTTGTAGTACTCATTGTGAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGTTGCGTTTATTTCCCCTCACCCGGTa                                                                                                                               <  1:231347/149‑1 (MQ=255)             gcCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCATTGTGAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCGCCCGGTAGGGgc                                                                                                                          >  2:323131/1‑149 (MQ=255)             gcCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCATTGTGAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGAGGGGgc                                                                                                                          >  2:21390/1‑149 (MQ=255)              ccATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCATTGTGAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGgcg                                                                                                                         >  2:365477/1‑149 (MQ=255)                    ggCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCATTGTGAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGcgagcgag                                                                                                                   <  2:415161/148‑1 (MQ=255)                      cACGGTTGATATGGCGATACCACGGGAGATTCGGGGTTGTAGTACTCATTGTGAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAggg                                                                                                                 <  1:170617/149‑1 (MQ=255)                        cGGTTGAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCATTGTGAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGa                                                                                                               <  1:435207/149‑1 (MQ=255)                                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)                                                                                                                                            |                                                                                                                                              CGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCATTGTGAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACAACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAACGGTTTGCGGCACAGCGGCACAGAAACG  >  NC_000913/3372503‑3372783

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