New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2871564 =	23 (0.740)	5 (0.170)	5/228	3.1	16.5%	coding (217/480 nt)	ispF	2C‑methyl‑D‑erythritol 2,4‑cyclodiphosphate synthase
	?	NC_000913	2871614 =	29 (0.990)					coding (167/480 nt)	ispF	2C‑methyl‑D‑erythritol 2,4‑cyclodiphosphate synthase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TAGTGACATCGACGTTGCCAAGGGTATAACCCTTCGCCTGAATACGACGCCAGGCTTCGCGTAGCAGCTCGCGGCTATCGGCACCTTTAAATGCCGGATCGGTATCCGGGAACAGCTTGCCGATATCCCCCAGCGCCGCCGCGCCAAGCAATGCATCGGTCAACGCATGGAGCGCCACGTCGCCATCAGAATGCGCCAGCAATCCTTTTTCGTAAGGAATGCGTACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGT  >  NC_000913/2871489‑2871758                                                                                                                              |                                                                                                                                                 tAGTGACATCGACGTTGCCAAGGGTATAACCCTTCGCCTGAATACGACGCCAGGCTTCGCGTAGCAGCTCGCGGCTATCGGCACCTTTAAATGCCGGATCGGTATCCGGGAACAGCTTGCCGATATCCCCCAGCGCCGCCGCGCCAAGc                                                                                                                           <  6:92502/149‑1 (MQ=255)        aTCGACGTTGCCAAGGGTATAACCCTTCGCCTGAATACGACGCCAGGCTTCGCGTAGCAGCTCGCGGCTATCGGAACCTTTAAATGCCGGATCGGTATCCGGGAACAGCTTGCCGATATCCCCCAGCGCCGCCGCGCCAAGCAATGCAt                                                                                                                    >  8:40803/1‑149 (MQ=255)         tCGACGTTGCCAAGGGTATAACCCTTCGCCTGAATACGACGCCAGGCTTCGCGTAGCAGCTCGCGGCTATCGGCACCTTTAAATGCCGGATCGGTATCCGGGAACAGCTTGCCGATATCCCCCAGCGCCGCCGCGCCAAGCAATGCATc                                                                                                                   >  8:135732/1‑149 (MQ=255)                    aaGGGTATAACCCTTCGCCTGAATACGACGCCAGGCTTCGCGTAGCAGCTCGCGGCTATCGGCACCTTTAAATGCCGGATCGGTATCCGGGAACAGCTTGCCGATATCCCCCAGCGCCGCCGCGCCAAGCAATGCATCGGTCAACGCAt                                                                                                        >  3:98976/1‑149 (MQ=255)                            aaCCCTTCGCCTGAATACGACGCCAGGCTTCGCGTAGCAGCTCGCGGCTATCGGCACCTTTAAATGCCGGATCGGTATCCGGGAACAGCTTGCCGATATCCCCCAGCGCCGCCGCGCCAAGCAATGCATCGGTCAACGCATGGAgcgc                                                                                                 >  1:133850/1‑148 (MQ=255)                                      cTGAATACGACGCCAGGCTTCGCGTAGCAGCTCGCGGCTATCGGCACCTTTAAATGCCGGATCGGTATCCGGGAACAGCTTGCCGATATCCCCCAGCGCCGCCGCGCCAAGCaa                                                                                                                         >  7:105566/1‑114 (MQ=255)                                      cTGAATACGACGCCAGGCTTCGCGTAGCAGCTCGCGGCTATCGGCACCTTTAAATGCCGGATCGGTATCCGGGAACAGCTTGCCGATATCCCCCAGCGCCGCCGCGCCAAGCaa                                                                                                                         <  8:105566/114‑1 (MQ=255)                                        gAATACGACGCCAGGCTTCGCGTAGCAGCTCGCGGCTATCGGCACCTTTAAATGCCGGATCGGTATCCGGGAACAGCTTGCCGATATCCCCCAGCGCCGCCGCGCCAAGCAATGCATCGGTCAACGCATGGAGCGCCACGTCGCCATCa                                                                                    <  4:98976/149‑1 (MQ=255)                                                  ccAGGCTTCGCGTAGCAGCTCGCGGCTATCGGCACCTTTAAATGCCGGATCGGTATCCGGGAACAGCTTGCCGATATCCCCCAGCGCCGCCGCGCCAAGCAATGCATCGGTCAACGCATGGAGCGCCACGTCGCCATCAGAATGCGCCa                                                                          <  7:135732/149‑1 (MQ=255)                                                      gCTTCGCGTAGCAGCTCGCGGCTATCGGCACCTTTAAATGCCGGATCGGTATCCGGGAACAGCTTGCCGATATCCCCCAGCGCCGCCGCGCCAAGCAATGCATCGGTCAACGCATGGAGCGCCACGTCGCCATCAGAATGCGCCAGCaa                                                                      <  5:34611/149‑1 (MQ=255)                                                        ttCGCGTAGCAGCTCGCGGCTATCGGCACCTTTAAATGCCGGATCGGTATCCGGGAACAGCTTGCCGATATCCCCCAGCGCCGCCGCGCCAAGCAATGCATCGGTCAACGCATGGAGCGCCa                                                                                               >  2:61553/1‑122 (MQ=255)                                                                gcagcTCGCGGCTATCGGCACCTTTAAATGCCGGATCGGTATCCGGGAACAGCTTGCCGATATCCCCCAGCGCCGCCGCGCCAAGCAATGCATCGGTCAACGCATGGAGCGCCACGTCGCCATCAGAATGCGCCAGCAATCCTTTTTc                                                             <  2:133850/148‑1 (MQ=255)                                                                        cgGCTATCGGCACCTTTAAATGCCGGATCGGTATCCGGGAACAGCTTGCCGATATCCCCCAGCGCCTCCGCGCCAAGCAATGCATCGGTCAACGCATGGAGCGCCACGTCGCCATCAGAATGCGCCAGCAATCCTTTTTCGTAAGGAAt                                                    <  7:40803/149‑1 (MQ=255)                                                                         ggCTATCGGCACCTTTAAATGCCGGATCGGTATCCGGGAACAGCTTGCCGATATCCCCCAGCGCCGCCGCGCCAAGCAATGCATCGGTCAACGCATGGAGCGCCACGTCGCCATCAGAATGCGCCAGCAATCCTTTTTCGTAAGGAATg                                                   >  1:1902/1‑149 (MQ=255)                                                                            tATCGGCACCTTTAAATGCCGGATCGGTATCCGGGAACAGCTTGCCGATATCCCCCAGCGCCGCCGCGCCAAGCAATGCATCGGTCAACGCATGGAGCGCCACGTCGCCATCAGAATGCGCCAGCAATCCTTTTTCGTAAGGAAt                                                    <  5:90339/145‑1 (MQ=255)                                                                            tATCGGCACCTTTAAATGCCGGATCGGTATCCGGGAACAGCTTGCCGATATCCCCCAGCGCCGCCGCGCCAAGCAATGCATCGGTCAACGCATGGAGCGCCACGTCGCCATCAGAATGCGCCAGCAATCCTTTTTCGTAAGGAAt                                                    >  6:90339/1‑145 (MQ=255)                                                                                      tttAAATGCCGGATCGGTATCCGGGAACAGCTTGCCGATATCCCCCAGCGCCGCCGCGCCAAGCAATGCATCGGTCAACGCATGGAGCGCCACGTCGCCATCAGAATGCGCCAGCAATCCTTTTTCGTAAg                                                        >  2:112884/1‑131 (MQ=255)                                                                                      tttAAATGCCGGATCGGTATCCGGGAACAGCTTGCCGATATCCCCCAGCGCCGCCGCGCCAAGCAATGCATCGGTCAACGCATGGAGCGCCACGTCGCCATCAGAATGCGCCAGCAATCCTTTTTCGTAAg                                                        <  1:112884/131‑1 (MQ=255)                                                                                         aaaTGACGGCTCTTTATCCGGGAACAGCTTGCCGAAAACCCCCAGCACCTCCCCGCCAACCAATGCAACGTTCAACGCATGGAGCGCCCCGTCGCCATAAGAATGCGCCAGCAATCCATAATCGTAAGCAATGCGTACGCCACCAATg                                    >  8:75426/1‑148 (MQ=255)                                                                                            tGCCGGATCGGTATCCGGGAACAGCTTGCCGATATCCCCCAGCGCCGCCGCGCCAAGCAATGCATCGGTCAACGCATGGAGCGCCACGTCGCCATCAGAATGCGCCAg                                                                         <  6:74512/108‑1 (MQ=255)                                                                                            tGCCGGATCGGTATCCGGGAACAGCTTGCCGATATCCCCCAGCGCCGCCGCGCCAAGCAATGCATCGGTCAACGCATGGAGCGCCACGTCGCCATCAGAATGCGCCAg                                                                         >  5:74512/1‑108 (MQ=255)                                                                                              ccGGATCGGTATCCGGGAACAGCTTGCCGATATCCCCCAGCGCCGCCGCGCCAAGCAATGCATCGGTCAACGCATGGAGCGCCACGTCGcc                                                                                        >  2:28469/1‑91 (MQ=255)                                                                                              ccGGATCGGTATCCGGGAACAGCTTGCCGATATCCCCCAGCGCCGCCGCGCCAAGCAATGCATCGGTCAACGCATGGAGCGCCACGTCGcc                                                                                        <  1:28469/91‑1 (MQ=255)                                                                                                     ggTATCCGGGAACAGCTTGCCGATATCCCCCAGCGCCGCCGCGCCAAGCAATGCATCGGTCAACGCATGGAGCGCCACGTCGCCATCAGAATGCGCCAGCAATCCTTTTTCGTAAGGaa                                                     <  7:60823/119‑1 (MQ=255)                                                                                                     ggTATCCGGGAACAGCTTGCCGATATCCCCCAGCGCCGCCGCGCCAAGCAATGCATCGGTCAACGCATGGAGCGCCACGTCGCCATCAGAATGCGCCAGCAATCCTTTTTCGTAAGGaa                                                     >  8:60823/1‑119 (MQ=255)                                                                                                           cGGGAACAGCTTGCCGATATCCCCCAGCGCCGCCGCGCCAAGCAATGCATCGGTCAACGCATGGAGCGCCACGTCGCCATCAGAATGCGCCAGCAATCCTTTTTCGTAAGGAATGCGTACGCCACCAATGATAAtt                              >  6:134291/1‑136 (MQ=255)                                                                                                           cGGGAACAGCTTGCCGATATCCCCCAGCGCCGCCGCGCCAAGCAATGCATCGGTCAACGCATGGAGCGCCACGTCGCCATCAGAATGCGCCAGCAATCCTTTTTCGTAAGGAATGCGTACGCCACCAATGATAAtt                              <  5:134291/136‑1 (MQ=255)                                                                                                                     ttGCCGATATCCCCCAGCGCCGCCGCGCCAAGCAATGCATCGGTCAACGCATGGAGCGCCACGTCGCCATCAGAATGCGCCAGCAATCCtttt                                                               <  6:13860/93‑1 (MQ=255)                                                                                                                     ttGCCGATATCCCCCAGCGCCGCCGCGCCAAGCAATGCATCGGTCAACGCATGGAGCGCCACGTCGCCATCAGAATGCGCCAGCAATCCtttt                                                               >  5:13860/1‑93 (MQ=255)                                                                                                                          gATATCCCCCAGCGCCGCCGCGCCAAGCAATGCATCGGTCAACGCATGGAGCGCCACGTCGCCATCAGAAt                                                                                <  7:5448/71‑1 (MQ=255)                                                                                                                          gATATCCCCCAGCGCCGCCGCGCCAAGCAATGCATCGGTCAACGCATGGAGCGCCACGTCGCCATCAGAAt                                                                                >  8:5448/1‑71 (MQ=255)                                                                                                                          gATATCCCCCAGCGCCGACGCGCCAAGCAATGCATCGGTCAACTCATGGAGCGCCACGTCGCCATCAGAATTCGCCAGCAATCCTTTTTCGTAAGGAATGCGTACGCCACGAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGt  <  7:75426/149‑1 (MQ=255)                                                                                                                              |                                                                                                                                                 TAGTGACATCGACGTTGCCAAGGGTATAACCCTTCGCCTGAATACGACGCCAGGCTTCGCGTAGCAGCTCGCGGCTATCGGCACCTTTAAATGCCGGATCGGTATCCGGGAACAGCTTGCCGATATCCCCCAGCGCCGCCGCGCCAAGCAATGCATCGGTCAACGCATGGAGCGCCACGTCGCCATCAGAATGCGCCAGCAATCCTTTTTCGTAAGGAATGCGTACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGT  >  NC_000913/2871489‑2871758

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