New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_tpiA	= 4511425	41 (1.120)	5 (0.140)	4/270	4.6	11.3%	coding (906/957 nt)	fecD	ferric citrate ABC transporter permease
	?	NC_000913_3_pae_tpiA	= 4511432	39 (1.110)					coding (899/957 nt)	fecD	ferric citrate ABC transporter permease
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCG  >  NC_000913_3_pae_tpiA/4511277‑4511570                                                                                                                                                     |                                                                                                                                                  gggCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCg                                                                                                                                                   <  2:112934/149‑1 (MQ=255)       tGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTCGCGGGTGTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGc                                                                                                                                             >  1:474119/1‑149 (MQ=255)        gATCGTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTTTCCCTGAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCa                                                                                                                                            <  2:95910/149‑1 (MQ=255)                   ttGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCggg                                                                                                                                 >  2:471396/1‑149 (MQ=255)                       cAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCt                                                                                                                             >  1:279078/1‑149 (MQ=255)                          tGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAGGCAACCCGACAAACCCCGGCGCCCCGATAATGGCGGTCAGCACGCCAACCGGGAGCTccg                                                                                                                          >  1:250326/1‑148 (MQ=255)                                tGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTgggg                                                                                                                    <  1:232633/149‑1 (MQ=255)                                   aaCGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGAt                                                                                                                 <  2:636662/149‑1 (MQ=255)                                           tAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAAtt                                                                                                         <  1:370992/149‑1 (MQ=255)                                           tAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAAtt                                                                                                         <  2:279078/149‑1 (MQ=255)                                              gTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCAGTCTCACAAGCAACCAGACACACCACGGCGCACCGATAATGGCGGTCGGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATtctc                                                                                                      >  1:275720/1‑149 (MQ=255)                                                            ccGTAACTGACCGTCAGATTTACAGTTCGTAAAGGAATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCCCGCCCACCGGGAGCCACAGTGGGGGATGACTAATTCTCGCCGGCAGATCGGc                                                                                        >  2:99798/1‑149 (MQ=255)                                                               tAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAac                                                                                     <  1:211123/149‑1 (MQ=255)                                                                   tGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATtctc                                                                                                      <  1:475334/128‑1 (MQ=255)                                                                   tGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATtctc                                                                                                      >  2:475334/1‑128 (MQ=255)                                                                   tGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCCGACCGCCAgccccc                                                                                 >  2:57963/1‑149 (MQ=255)                                                                         tCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCaaaaac                                                                           >  1:124513/1‑145 (MQ=255)                                                                                  cAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCt                                                                  <  1:334529/149‑1 (MQ=255)                                                                                    gTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGGGGGGGATGAATAACACTCGCCAGCAGATCGGGCACCACCAACAGCAACGCAcccac                                                                >  2:90839/1‑146 (MQ=255)                                                                                     ttCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCCCCAACAGCAACGCACCTGTc                                                               >  2:593929/1‑149 (MQ=255)                                                                                     ttCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTc                                                               >  1:159290/1‑149 (MQ=255)                                                                                     ttCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACAGGGAGCTCCAGCGGGGGATGAACAAGTCTCGCCAGCAGCTCCGCCACCACCACCACGAACGAACCtctc                                                               >  2:97336/1‑149 (MQ=255)                                                                                         tAAAGCCATTTATCGCATTCTCACAAGCAACCAGACAACACACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAggg                                                           <  2:474119/149‑1 (MQ=255)                                                                                              tCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGCCAGGGCCGac                                                      >  2:596125/1‑148 (MQ=255)                                                                                                        cATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGgcagca                                            >  2:184950/1‑149 (MQ=255)                                                                                                          cgacctcaAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCCCCAGTGGGGCATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGt                                          <  1:296486/143‑1 (MQ=255)                                                                                                             tcACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTg                                       >  1:586113/1‑149 (MQ=255)                                                                                                               acaAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCg                                     <  2:159290/149‑1 (MQ=255)                                                                                                               acaAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCACAAGCTCCAGGGGGGGATGAATAATTCTCGCCAGCAGATCGGCCACCACCAACAGCAACGCACCTGCCAGCGCCCAAACAGGccgcagcccccg                                     >  2:42878/1‑143 (MQ=255)                                                                                                                   gCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCTGAGGTGTAGTAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGa                                 <  1:471396/149‑1 (MQ=255)                                                                                                                      aCCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGt                              >  2:53086/1‑149 (MQ=255)                                                                                                                        cAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAacca                                                                                   >  1:278845/1‑94 (MQ=255)                                                                                                                        cAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAacca                                                                                   <  2:278845/94‑1 (MQ=255)                                                                                                                             aaaCCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCgg                       <  1:53086/149‑1 (MQ=255)                                                                                                                                      cgcACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTAc              <  1:497081/149‑1 (MQ=255)                                                                                                                                      cgcACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTAc              <  1:323427/149‑1 (MQ=255)                                                                                                                                      cgcACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTAc              <  1:611635/149‑1 (MQ=255)                                                                                                                                          ccGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCCGGCTGCGGTGACGTCCACCGGTGATGGTACGcat          >  1:94215/1‑149 (MQ=255)                                                                                                                                           cGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGcctc         >  1:298354/1‑149 (MQ=255)                                                                                                                                           cGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCCGTGACGACCACCGGCGATGCAACCcatc         >  2:36203/1‑149 (MQ=255)                                                                                                                                                aTGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATg    <  1:518718/149‑1 (MQ=255)                                                                                                                                                  ggCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCg  <  1:231630/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCG  >  NC_000913_3_pae_tpiA/4511277‑4511570

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