New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_tpiA	3423672 =	NA (NA)	4 (0.130) +ACATCAA	4/270	3.9	4.5%	noncoding (116/120 nt)	rrfD	5S ribosomal RNA of rrnD operon
	?	NC_000913_3_pae_tpiA	= 4040683	89 (2.850)					intergenic (+47/+223)	rrfA/mobB	5S ribosomal RNA of rrnA operon/molybdopterin‑guanine dinucleotide biosynthesis protein B
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGACCG  >  NC_000913_3_pae_tpiA/3423524‑3423807                                                                                                                                                     |                                                                                                                                        gcgcTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGcc                                                                                                                                         <  1:304815/149‑1 (MQ=255)    ctccTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTgg                                                                                                                                      <  1:209852/146‑1 (MQ=255)    cTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACGAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTgg                                                                                                                                      <  1:363956/149‑1 (MQ=255)    cTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTgg                                                                                                                                      <  1:177348/149‑1 (MQ=255)     tACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGc                                                                                                                                     <  1:138931/149‑1 (MQ=255)     tACTGCCGCCAAACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTTCCAAGGATGCGCTCTACCAACTGAGCCATATCAGCACCCTAAATATAACGCCTCGc                                                                                                                                     >  2:520909/1‑149 (MQ=37)           cgccAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTcc                                                                                                                               >  1:222683/1‑149 (MQ=255)                      ttCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCa                                                                                                                    <  1:117659/149‑1 (MQ=255)                       ttgttgttAATCTGCCGGACTTGAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCATCCTTACCAAGGGTGCGCTCTACCAACTGCGCCATATCAGCACGCTAAATTTGATGCCTGGTAGTTCCCTACTCTCGCAt                                                                                                                   <  2:178438/142‑1 (MQ=37)                           ttCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGAGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAACTTTGATGCCTGGCAGTTCCCTACTCTCGCATgggg                                                                                                               >  2:82405/1‑149 (MQ=255)                                  cTGCCTAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGCGAACTCACCATTACCAAGGGTGCGCGCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGAcccc                                                                                                        <  1:255411/149‑1 (MQ=38)                                           tttAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATggg                                                                                                                <  1:146585/132‑1 (MQ=255)                                           tttAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATggg                                                                                                                >  2:146585/1‑132 (MQ=255)                                                                     ccAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCCCTACGGCTTTGCACTTCTgc                                                                     >  2:519023/1‑148 (MQ=255)                                                                             gAACTGGGGACCTCACCCATACCAAGGGTGCGCACTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTAGCATGAGTAGACCCCACACTACCATCAGCGCTACGGCGTTTCACTTCTGAGTTCGGCa                                                             >  1:456837/1‑149 (MQ=21)                                                                                   tggACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATgg                                                          <  1:319989/145‑1 (MQ=255)                                                                                   gggACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATgg                                                          >  2:319989/1‑146 (MQ=255)                                                                                         tCACCCTTACCAGGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACGCTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTg                                                 >  2:37402/1‑149 (MQ=255)                                                                                         tCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTg                                                 >  2:190334/1‑149 (MQ=255)                                                                                           aCCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTggg                                               >  2:46245/1‑149 (MQ=255)                                                                                               ttACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGacca                                           >  1:523959/1‑149 (MQ=255)                                                                                                            gcTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTAAGTTCGGCATGGGGTCAGATGGGACCACCGCACTACGgcc                              >  2:241831/1‑149 (MQ=255)                                                                                                             ctctACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTAGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTACGGTGAGAGCACCGCGCCACGgccc                             >  2:207523/1‑148 (MQ=17)                                                                                                                  ccAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGATGGGCCCACCGCGCTACGGCCGCCAgg                        >  2:34650/1‑149 (MQ=255)                                                                                                                   cAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATg                                                           <  1:261105/113‑1 (MQ=255)                                                                                                                   cAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATg                                                           >  2:261105/1‑113 (MQ=255)                                                                                                                   cAACTGAGCCATATCAGCACGCTAAATTTGATGCCCCGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGc                       <  1:207523/149‑1 (MQ=255)                                                                                                                           ccATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTg               >  1:208499/1‑149 (MQ=255)                                                                                                                                  cgtaaaatAAATTGGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAACTTCTGTTTTATc        >  2:506833/8‑149 (MQ=255)                                                                                                                                  aGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGttt            <  2:267254/145‑1 (MQ=37)                                                                                                                                  aGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGttt            >  1:267254/1‑145 (MQ=37)                                                                                                                                       gCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCACAACGGGCGCCAGGCGAATTCTGTTTTATTAGAcc   >  1:344143/1‑149 (MQ=11)                                                                                                                                        ttAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTGCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGACCg  <  2:222683/148‑1 (MQ=14)                                                                                                                                         tAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGt              <  1:479626/136‑1 (MQ=31)                                                                                                                                         tAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGt              >  2:479626/1‑136 (MQ=31)                                                                                                                                                     |                                                                                                                                        GCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGACCG  >  NC_000913_3_pae_tpiA/3423524‑3423807

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.