Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	NC_002947	4,712,487	Δ33 bp	38.7%	coding (673‑705/2115 nt)	PP_4171 →	hypothetical protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 4712486	14 (0.720)	9 (0.480)	9/250	NT	38.7%	coding (672/2115 nt)	PP_4171	hypothetical protein
	?	NC_002947	4712520 =	15 (0.800)					coding (706/2115 nt)	PP_4171	hypothetical protein

CTTGGGCTGGGACGGTGGCGACATGCGCCTGTCGCTGTCGGGCCAACAGGACGACGGCTTCGATGAGAACCAGTTCGGCCAGGACTACCGCGACAGCCGCCGCCTGACCCGCTTCAACCTCAGCGCCAACCACAACCTGGCGGTCGACCAGACCCTGGAGTGGCAACTGGCGGCCAAGGAAGGCAGCAATCAGCGGCCCTACACCTACCAACCGGTGTTCCAGCAGTACGAGTCCGGGAACAATGCCGACGTCAACGCCAAGGACTATGCGGGCTCGTTGCGCTGGA  >  NC_002947/4712381‑4712667                                                                                                                                            |                                                                                                                                                    ggacggtatcaacgacttctatgccagccacggcttgggctgggacggtggcgacatgcgcctgtcgctgtcgggccaacaggacgacggcttcgatgagaaccagttcggccaggactaccgcgacagccgccgcctgGCGGTCGACCa                                                                                                                                           <  1:144342‑M2/11‑1 (MQ=255)       tatcaacgacttctatgccagccacggcttgggctgggacggtggcgacatgcgcctgtcgctgtcgggccaacaggacgacggcttcgatgagaaccagttcggccaggactaccgcgacagccgccgcctgGCGGTCGACCAGACCCt                                                                                                                                     >  1:234121‑M2/134‑150 (MQ=255)         tcaacgacttctatgccagccacggcttgggctgggacggtggcgacatgcgcctgtcgctgtcgggccaacaggacgacggcttcgatgagaaccagttcggccaggactaccgcgacagccgccgcctgGCGGTCGACCAGACCCTgg                                                                                                                                   >  1:435485‑M2/132‑150 (MQ=255)                    gACATGCGCCTGTCGCTGTCGGGCCAACAGGACGACGGCTTCGATGAGAACCAGTTCGGCCAGGACTACCGCGACAGCCGCCGCCTGACCCGCTTCAACCTCAGCGCCAACCACAACCTGGCGGTCGACCAGACCCTGGAGTGGCAACTg                                                                                                                        >  2:111313/1‑150 (MQ=255)                        tGCGCCTGTCGCTGTCGGGCCAACAGGACGACGGCTTCGATGAGAACCAGTTCGGCCAGGACTACCGCGACAGCCGCCGCCTGACCCGCTTCAACCTCAGCGCCAACCACAACCTGGCGGTCGACCAGACCCTGGAGTGGCAACTggcgg                                                                                                                    <  2:246189/150‑1 (MQ=255)                           gcatgtcgctgtcgGGCCAACAGGACGACGGCTTCGATGAGAACCAGTTCGGCCAGGACTACCGCGACAGCCGCCGCCTGACCCGCTTCAACCTCAGCGCAAACCACAACCTGGCGGTCGACCAGACCCTGGAGTGGCAACTGGCGGCCa                                                                                                                 >  2:269479/4‑150 (MQ=255)                               gtcgctgtcgGGCCAACAGGACGACGGCTTCGATGAGAACCAGTTCGGCCAGGACTACCGCGACAGCCGCCGCCTGACCCGCTTCAACCTCAGCGCCAACCACAACCTGGCGGTCGACCAGACCCTGGAGTGGCAACTGGCGGCCaagga                                                                                                             >  1:198206/1‑150 (MQ=255)                                        gggCCAACAGGACGACGGCTTCGATGAGAACCAGTTCGGCCAGGACTACCGCGACAGCCGCCGCCTGACCCGCTTCAACCTCAGCGCCAACCACAACCTGGCGGTCGACCAGACCCTGGAGTGGCAACTGGCGGCCAAGGAAGGCAGCaa                                                                                                    >  1:19111/1‑150 (MQ=255)                                                    cgacatgcgcctgtcgctgtcgggccaacaggacgacggcttcgatgagaaccagttcggccaggactaccgcgacagccgccgcctgGCGGTCGACCAGACCCTGGAGTGGCAACTGGCGGCCAAGGAAGGCAGCAATCAGCGGCCCTa                                                                                        <  2:435485‑M2/62‑1 (MQ=255)                                                      acatgcgcctgtcgctgtcgggccaacaggacgacggcttcgatgagaaccagttcggccaggactaccgcgacaggaccacccgtGCGGTCGACCAGACCCTGGCGTGGAAACCGGCGGCCGAGGAAGGGAGCAATCAGCGGCCCTaca                                                                                      >  1:86129‑M2/87‑150 (MQ=255)                                                           ttCAATGAGAAAACGTTCGGCCAGGACTACCGCGACAGCCGCCGCCTGACCCGCTTCAACCTCAGCGCCAACCACAACCTGGCGGTCGACCAGACCCTGGAGTGGCAACTGGCGGCCAAGGAAGGCAGCAATCAGCGGCCCTACACCTAc                                                                                 <  2:397358/150‑1 (MQ=255)                                                                               acaggacgacggcttcgatgagaaccagttcggccaggactaccgcgacagccgccgcctgGCGGTCGACCAGACCCCGGAGCGGCAACTGGCGGCGAAGGAAGGCAGCAATCAGCGGCCTTACACCTACAAACCGGTGCTCCAGCAGTa                                                             >  1:471037‑M2/62‑150 (MQ=255)                                                                                          gcgACAGCCGCCGCCTGACCCGCTTCAACCTCAGCGCCAACCACAACCTGGCGGTCGACCAGACCCTGGAGTGGCAACTGGCGGCCAAGGAAGGCAGCAATCAGCGGCCCTACACCTACCAACCGGTGTTCCAGCAGTACGAGTCCGGGa                                                  >  2:420242/1‑150 (MQ=255)                                                                                           cgACAGCCGCCGCCTGACCCGCTTCAACCTCAGCGACAACCCCAAACTGGCGGTCTACCAGACCCTGGAGTGATCACGGTGGGCCAAGGAAGGCAGCAATCAGCGGCCCTACACCTACCAACCGGTGTTCCAGCAGTACGAGTCCGGGaa                                                 <  1:111313/150‑1 (MQ=255)                                                                                              cAGCCGCCGCCTGACCCGCTTCAACCTCAGCGCCAACCACAACCTGGCGGTCGACCAGACCCTGGAGTGGCAACTGGCGGCCAAGGAAGGCAGCAATCAGCGGCCCTACACCTACCAACCGGTGTTCCAGCAGTACGAGTCCGGGAACaa                                              >  2:403652/1‑150 (MQ=255)                                                                                               gatgagaaccagttcggccaggactaccgcgacagccgccgcctgGCGGTCGACCAGACCCTGGAGTGGCAACTGGCGGCCAAGGAAGGCAGCAATCAGCGGCCCTACACCTACCAACCGGTGTTCCAGCAGTACGAGTCCGGGAACAAt                                             >  1:72930‑M2/46‑150 (MQ=255)                                                                                                    ccgccTGACCCGCTTCAACCTCAGCGCCAACCACAACCTGGCGGTCGACCAGACCCTGGAGTGGCAACTGGCGGCCACGGAAGGCAGCAATCAGCGGCCCTACACCTACCAACCGGTGTTCCAGCAGTACGAGTCCGGGAACAATGCCGa                                        <  1:403652/150‑1 (MQ=255)                                                                                                     cgccTGACCCGCTTCAACCTCAGCGCCAACCACAACCTGGCGGTCGACCAGACCCTGGAGTGGCAACTGGCGGCCAAGGAAGGCAGCAATCAGCGGCCCTACACCTACCAACCGGTGTTCCAGCAGTACGAGTCCGGGAACAATGCCGAc                                       <  2:198206/150‑1 (MQ=255)                                                                                                         tGACCCGCTTCAACCTCAGCGCCAACCACAACCTGGCGGTCGACCAGACCCTGGAGTGGCAACTGGCGGCCAAGGAAGGCAGCAATCAGCGGCCCTACACCTACCAACCGGTGTTCCAGCAGTACGAGTCCGGGAACaa                                              <  1:486552/139‑1 (MQ=255)                                                                                                         tGACCCGCTTCAACCTCAGCGCCAACCACAACCTGGCGGTCGACCAGACCCTGGAGTGGCAACTGGCGGCCAAGGAAGGCAGCAATCAGCGGCCCTACACCTACCAACCGGTGTTCCAGCAGTACGAGTCCGGGAACaa                                              >  2:486552/1‑139 (MQ=255)                                                                                                                ccaggactaccgcgacagccgccgcctgGCGGTCGACCAGACCCTGGAGTGGCAACTGGCGGCCAAGGAAGGCAGCAATCAGCGGCCCTACACCTACCAACCGGTGTTCCAGCAGTACGAGTCCGGGAACAATGCCGACGTCAACGCCaa                            <  2:279425‑M2/122‑1 (MQ=255)                                                                                                                    gactaccgcgacagccgccgcctgGCGGTCGACCAGACCCTGGAGTGGCAACTGGCGGCCAAGGAAGGCAGCAATCAGCGGCCCTACACCTACCAACCGGTGTTCCAGCAGTACGAGTCCGGGAACAATGCCGACGTCAACGCCAAGGAc                        <  2:72930‑M2/126‑1 (MQ=255)                                                                                                                             gcCAACCACAACCTGGCGGTCGACCAGACCCTGGAGTGGCAACTGGCGGCCAAGGAAGGCAGCAATCAGCGGCCCTACACCTACCAACCGGTGTTCCAGCAGTACGAGTCCGGGAACAATGCCGACGTCAACGCCAAGGACTATGCGGGc               >  2:395148/1‑150 (MQ=255)                                                                                                                                 aCCACAACCTGGCGGTCGACCAGACCCTGGAGTGGCAACTGGCGGCCAAGGAAGGCAGCAATGAGTGGCCCTACACCTACCAACCGGTGTTCCAGCAGTACGAGTCCGGGCACAATGCCGACGTCAACTCCAAGGACTATGCGGGCTCGt           >  2:47821/1‑150 (MQ=255)                                                                                                                                    acaACCTGGCGGTCGACCAGACCCTGGAGTGGCAACTGGCGGCCAAGGAAGGCAGCAATCAGCGGCCCTACACCTACCAACCGGTGTTCCAGCAGTACGAGTCCGGGAACAATGCCGACGTCAACGCCAAGGACTATGCGGGCTCGTTgc        <  2:212094/150‑1 (MQ=255)                                                                                                                                          tgGCGGTCGACCAGACCCTGGAGTGGCAACTGGCGGCCAAGGAAGGCAGCAATCAGCGGCCCTACACCTACCAACCGGTGTTCCAGCAGTACGAGTCCGGGAACAATGCCGACGTCAACGCCAAGGACTATGCGGGCTCGTTGCGCTCGa  <  2:19111‑M2/148‑1 (MQ=255)                                                                                                                                            |                                                                                                                                                    CTTGGGCTGGGACGGTGGCGACATGCGCCTGTCGCTGTCGGGCCAACAGGACGACGGCTTCGATGAGAACCAGTTCGGCCAGGACTACCGCGACAGCCGCCGCCTGACCCGCTTCAACCTCAGCGCCAACCACAACCTGGCGGTCGACCAGACCCTGGAGTGGCAACTGGCGGCCAAGGAAGGCAGCAATCAGCGGCCCTACACCTACCAACCGGTGTTCCAGCAGTACGAGTCCGGGAACAATGCCGACGTCAACGCCAAGGACTATGCGGGCTCGTTGCGCTGGA  >  NC_002947/4712381‑4712667

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