New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_pgi	= 2641725	16 (0.920)	4 (0.250)	4/256	1.8	20.0%	coding (80/1119 nt)	ispG	1‑hydroxy‑2‑methyl‑2‑(E)‑butenyl 4‑diphosphate synthase
	?	NC_000913_3_pae_pgi	= 2641730	17 (1.040)					coding (75/1119 nt)	ispG	1‑hydroxy‑2‑methyl‑2‑(E)‑butenyl 4‑diphosphate synthase
Rejected: Frequency below/above cutoff threshold.

TGAACGCTTCTGCCGCGTCCATCGTCGGTACGGATACACGGACGATATCAGCGCCAACGCGTTCCAGCGCCTTGATTTGATTGACCGTTGCTTCGACGTCTGTCGTACGCGTATTGGTCATGGACTGTACGGCGATGGGAGCACCATCGCCAATCGGCACATTCCCAACGTAAATACGTGTTGATTTTCTACGTTGAATTGGAGCCTGGTTATGCATGAAAAATCTCCCGCGTTACCCGTCTGTTACTGCGCCGGTGATTGTTCGGCATTGAGGGTCAGAC  >  NC_000913_3_pae_pgi/2641588‑2641868                                                                                                                                          |                                                                                                                                                tGAACGCTTCTGCCGCGTCCATCGACGGTACGGATACACGGACGATATCAGCGCCACCGCGTTCCAGCGCCTTGATTTGATTGACCGTTGCTTCGACGTCTGTCGTACGCGTATTGGTCATGGACTGTACGGCGATGGGAGCAccatcgc                                                                                                                                     >  3:88436‑M2/1‑143 (MQ=255)  gAACGCTTCTGCCGCGTCCATCGTCGGTACGGATACACGGACGATATCAGCGCCAACGCGTTCCAGCGCCTTGATTTGATTGACCGTTGCTTCGACGTCTGTCGTACGCGTATTGGTCATGGACTGTACGGCGATGGGAGCAccatcgcc                                                                                                                                    <  1:93191‑M2/150‑9 (MQ=255)     cGCTTCTGCCGCGTCCATCGTCGGTACGGATACACGGACGATATCAGCGCCAACGCGTTCCAGCGCCTTGATTTGATTGACCGTTGCTTCGACGTCTGTCGTACGCGTATTGGTCATGGACTGTACGGCGATGGtgctcccatcgccgta                                                                                                                                 <  2:183386‑M1/150‑17 (MQ=255)        ttCTGCCGCGTCCATCGTCGGTACGGATACACGGACGATATCAGCGCCAACGCGTTCCAGCGCCTTGATTTGATTGACCGTTGCTTCGACGTCTGTCGTACGCGTATTGGTCATGGACTGTACGGCGATGGGAGCACCATCGCCa                                                                                                                                   >  4:66668/1‑145 (MQ=255)        ttCTGCCGCGTCCATCGTCGGTACGGATACACGGACGATATCAGCGCCAACGCGTTCCAGCGCCTTGATTTGATTGACCGTTGCTTCGACGTCTGTCGTACGCGTATTGGTCATGGACTGTACGGCGATGGGAGCACCATCGCCa                                                                                                                                   <  3:66668/145‑1 (MQ=255)        ttCTGCCGCGTCCATCGTCGGTACGGATACACGGACGATATCAGCGCCAACGCGTTCCAGCGCCTTGATTTGATTGACCGTTGCTTCGACGTCTGTCGTACGCGTATTGGTCATGGACTGTACGGCGATGGGAGCACCATCGCCAATCgg                                                                                                                              <  3:10534/150‑1 (MQ=255)                      tcgtcgGTACGGATACACGGACGATATCAGCGCCAACGCGTTCCAGCGCCTTGATTTGATTGACCGTTGCTTCGACGTCTGTCGTACGCGTATTGGTCATGGACTGTACGGCGATGGGAGCACCATCGCCAATCGGCACATTCCCAACGt                                                                                                                >  2:15156/1‑150 (MQ=255)                            gTACGGATACACGGACGATATCAGCGCCAACGCGTTCCAGCGCCTTGATTTGATTGACCGTTGCTTCGACGTCTGTCGTACGCGTATTGGTCATGGACTGTACGGCGATGGtgctcccatcgccgtacagtccatgaccaatacgcgtac                                                                                                          <  1:190493‑M1/150‑40 (MQ=255)                            gTACGGATACACGGACGATATCAGCGCCAACGCGTTCCAGCGCCTTGATTTGATTGACCGTTGCTTCGACGTCTGTCGTACGCGTATTGGTCATGGACTGTACGGCGATGGGAGCACCATCGCCAATCGGCACATTCCCAACGTAAATAc                                                                                                          >  1:120969/1‑150 (MQ=255)                                       cGGACGATATCAGCGCCAACGCGTTCCAGCGCCTTGATTTGATTGACCGTTGCTTCGACGTCTGTCGTACGCGTATTGGTCATGGACTGTACGGCGATGGGAGCACCATCGCCAATCGGCACATTCCCAACGTAAATACGTGTTGAtttt                                                                                               <  4:76979/150‑1 (MQ=255)                                                                        ttgatttgatTGACCGTTGCTTCGACGTCTGTCGTACGCGTATTGGTCATGGACTGTACGGCGATGGGAGCACCATCGCCAATCGGCACATTCCCAACGTAAATACGTGTTGATTTTCTACGTTGAATTGGAGCCTGGTTATGCATGaaa                                                              <  2:120969/150‑1 (MQ=255)                                                                                 ttgaCCGTTGCTTCGACGTCTGTCGTACGCGTATTGGTCATGGACTGTACGGCGATGGGAGCACCATCGCCAATCGGCACATTCCCAACGTAAATACGTGTTGATTTTCTACGTTGAATTGGAGCCTGGTTATGCATGAAAAATCTCCcg                                                     >  3:91394/1‑150 (MQ=255)                                                                                       gTTGCTTCGACGTCTGTCGTACGCGTATTGGTCATGGACTGTACGGCGATGGGAGCACCATCGCCAATCGGCACATTCCCAACGTAAATACGTGTTGATTTTCTACGTTGAATTGGAGCCTGGTTATGCATGAAAAATCTCCCGCGTTAc                                               <  1:92981/150‑1 (MQ=255)                                                                                        ttGCTTCGACGTCTGTCGTACGCGTATTGGTCATGGACTATACGGCGATGGtgctcccatcgccgtacagtccatgacc                                                                                                                     <  1:16020‑M1/79‑29 (MQ=255)                                                                                        ttGCTTCGACGTCTGTCGTACGCGTATTGGTCATGGACTATACGGCGATGGtgctcccatcgccgtacagtccatgacc                                                                                                                     >  2:16020‑M1/1‑51 (MQ=255)                                                                                                         gTACGCGTATTGGTCATGGACTGTACGGCGATGGGAGCAccatcgccgtacagtccatgaccaatacgcgtacgacagacgtcgaagcaacggtcaatcaaatcaaggcgctggaacgcgttggcgctgatatcgtccgtgtatccgtac                             <  1:190493‑M2/150‑112 (MQ=255)                                                                                                          tACGCGTATTGGTCATGGACTGTACGGCGATGGGAGCACCATCGCCAATCGGCACATTCCCAACGTAAATACGTGTTGATTTTCTACGTTGAATTGGAGCCTGGTTATGCATGAAAAATCTCCCGCGTTACCCGTCTGTTACTGCGCCgg                            >  1:8600/1‑150 (MQ=255)                                                                                                             gcgTATTGGTCATGGACTGTACAGCGATGGGAGCACCATCGCCAATCGGCACATTCCCAACGTAAATACGTGTTGATTTTCTACGTTGAAt                                                                                    >  2:4821/1‑91 (MQ=255)                                                                                                             gcgTATTGGTCATGGACTGTACAGCGATGGGAGCACCATCGCCAATCGGCACATTCCCAACGTAAATACGTGTTGATTTTCTACGTTGAAt                                                                                    <  1:4821/91‑1 (MQ=255)                                                                                                               gTATTGGTCATGGACTGTACGGCGATGGGAGCACCATCGCCAATCGGCACATTCCCAACGTAAATACGTGTTGATTTTCTACGTTGAATTGGAGCCTGGTTATGCATGAAAAATCTCCCGCGTTACCCGTCTGTTACTGCGCCGGTGAtt                       >  1:169119/1‑150 (MQ=255)                                                                                                                   tGGTCATGGACTGTACGGCGATGGGAGCACCATCGCCAATCGGCACATTCCCAACGTAAATACGTGTTGATTTTCTACGTTGAATTGGAGCCTGGTTATGCATGAAAAATCTCCCGCGTTACCCGTCTGTTACTGCGCCGGTGATTGTTc                   <  1:15156/150‑1 (MQ=255)                                                                                                                    ggTCATGGACTGTACGGCGATGGGAGCAccatcgccgtatagtccatgaccaatacgcgtacgacagacgtcgaagcaa                                                                                         <  1:16020‑M2/79‑52 (MQ=255)                                                                                                                    ggTCATGGACTGTACGGCGATGGGAGCAccatcgccgtatagtccatgaccaatacgcgtacgacagacgtcgaagcaa                                                                                         >  2:16020‑M2/1‑28 (MQ=255)                                                                                                                            ccTGTACGGCGATGGGAGCAACATAGCAAATCGGAACATCCCCAACGTCAATACGTGTTGATTTTCTACGTTGAATTGGAGCCTGGTTATGTATGAAAAATCTCCCGCGTTACCCGTCTGTTACTGCGACGGCAATTGTTCTGCATTGAg          >  2:139217/2‑150 (MQ=255)                                                                                                                                tACGGCGATGGGAGCAccatcgccgtacagtccatgaccaatacgcgtacgacagacgtcgaagcaacggtcaatcaaatcaaggcgctggaacgcgttggcgctgatatcgtccgtgtatccgtaccgacgatggacgcggcagaagcg      <  2:183386‑M2/150‑135 (MQ=255)                                                                                                                                tACGGCGATGGGAGCACCATCGCCAATCGGCACATTCCCAACGTAAATACGTGTTGATTTTCTACGTTGAATTGGAGCCTGGTTATGCATGAAAAATCTCCCGCGTTACCCGTCTGTTACTGCGCCGGTGATTGTTCGGCATTGAGGGTc      >  3:51361/1‑150 (MQ=255)                                                                                                                                   ggCGATGGtgctcccatcgccgtacagtccatgaccaatacgcgtacgacagacgtcgaagcaacggtcaatcaaatcaaggcgctggaacgcgttggcgctgatatcgtccgtgtatccgtaccgacgatggacgcggcagaagcgttc   <  1:93191‑M1/150‑143 (MQ=255)                                                                                                                                    gCGATGGtgctcccatcgccgtacagtccatgaccaatacgcgtacgacagacgtcgaagcaacggtcaatcaaatcaaggcgctggaacgcggtggcgctgatatcgtccgtgtatccgtaccgtcgatggacgcggcagaagcgttca  >  3:88436‑M1/1‑7 (MQ=255)                                                                                                                                          |                                                                                                                                                TGAACGCTTCTGCCGCGTCCATCGTCGGTACGGATACACGGACGATATCAGCGCCAACGCGTTCCAGCGCCTTGATTTGATTGACCGTTGCTTCGACGTCTGTCGTACGCGTATTGGTCATGGACTGTACGGCGATGGGAGCACCATCGCCAATCGGCACATTCCCAACGTAAATACGTGTTGATTTTCTACGTTGAATTGGAGCCTGGTTATGCATGAAAAATCTCCCGCGTTACCCGTCTGTTACTGCGCCGGTGATTGTTCGGCATTGAGGGTCAGAC  >  NC_000913_3_pae_pgi/2641588‑2641868

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