BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_007779	3,261,134	Δ1,199 bp		insH–[tdcD]	insH, [tdcD]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_007779	3261134–3262328	3262332	5–1199	29 [0]	[0] 29	insH–[tdcD]	insH,[tdcD]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_007779	= 3261133	0 (0.000)	29 (0.650)	26/192	0.7	100%	intergenic (‑59/‑69)	tdcE/insH	pyruvate formate‑lyase 4/2‑ketobutyrate formate‑lyase/IS5 element protein
*	?	NC_007779	3262333 =	0 (0.000)	29 (0.650)	26/192	0.7	100%	coding (1183/1296 nt)	tdcD	propionate kinase/acetate kinase C, anaerobic

CAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAACTTTGCCTAAATGAATGGCATCCAAAGCAATCATTTTTTCTTCGTTAGTCGGAATAACGGCACAAATGACACGCGCATTTTCACTGGAAACAATTC  >  NC_007779/3262235‑3262428

ggcgtacagcttatcgctggtatcaatatctaccttcatgaaaaataatctctctacaatacttcaactaaattatgcaaattctgcgggcgcgttaaCtt                                                                                               >  1:482823‑M2/99‑101 (MQ=255)

 gcgtacagcttatcgctggtatcaatatctaccttcatgaaaaataatctctctacaatacttcaactaaattatgcaaattctgcgggcgcgttaaCttt                                                                                              <  3:1359145‑M2/4‑1 (MQ=255)

   gtacagcttatcgctggtatcaatatctaccttcatgaaaaataatctctctacaatacttcaactaaattatgcaaattctgcgggcgcgttaaCTTTGc                                                                                            <  3:907958‑M2/6‑1 (MQ=255)

        gcttctcgctggtatcactctctaccttcatgaaaaataatctctctacaatacttcaactaaattatgcaaattctgcgggcgcgttaaCTTTGCCTaaa                                                                                       >  4:127414‑M2/91‑101 (MQ=255)

          ttatcgctggtatcaatatctaccttcatgaaaaataatctctctacaatacttcaactaaattatgcaaattctgcgggcgcgttaaCTTTGCCTAaatg                                                                                     <  3:1355454‑M2/13‑1 (MQ=255)

              cgctggtatcaatatctaccttcatgaaaaataatctctctacaatacttcaactaaattatgcaaattctgcgggcgcgttaaCTTTGCCTAaatgaatg                                                                                 <  3:558290‑M2/17‑1 (MQ=255)

                   gtatcaatatctaccttcatgaaaaataatctctctacaatacttcaactaaattatgcaaattctgcgggcgcgttaaCTTTGCCTAAATGAATGGCATc                                                                            >  3:191394‑M2/80‑101 (MQ=255)

                       caatatctaccttcatgaaaaataatctctctacaatacttcaactaaattatgcaaattctgcgggcgcgttaaCTTTGCCTAAATGAATGGCATCCaaa                                                                        <  3:826502‑M2/26‑1 (MQ=255)

                       caatatctaccttcatgaaaaataatctctctacaatacttcaactaaattatgcaaattctgcgggcgcgttaaCTTTGCCTAAATGAATGGCATCCaaa                                                                        <  3:787587‑M2/26‑1 (MQ=255)

                            tctaccttcatgaaaaataatctctctacaatacttcaactaaattatgcaaattctgcgggcgcgttaaCTTTGCCTAAATGAATGGCATCCAAAGCAAt                                                                   <  3:1263050‑M2/31‑1 (MQ=255)

                                ccttcatgaaaaataatctctctacaatacttcaactaaattatgcaaattctgcgggcgcgttaaCTTTGCCTAAATGAATGGCATCCAAAGCAATCAtt                                                               <  3:1246230‑M2/35‑1 (MQ=255)

                                 cttcatgaaaaataatctctctacaatacttcaactaaattatgcaaattctgcgggcgcgttaaCTTTGCCTAAATGAATGGCATCCAAAGCAATCAttt                                                              <  1:1478007‑M2/36‑1 (MQ=255)

                                 cttcatgaaaaataatctctctacaatacttcaactaaattatgcaaattctgcgggcgcggtaaCTTTGCCTAAATGAATGGCATCCAAAGCAATCAttt                                                              >  3:284854‑M2/66‑101 (MQ=255)

                                  ttcatgaaaaataatctctctacaatacttcaactaaattatgcaaattctgcgggcgcgttaaCTTTGCCTAAATGAATGGCATCCAAAGCAATCAtttt                                                             >  2:91564‑M2/65‑101 (MQ=255)

                                  ttcatgaaaaataatctctctacaatacttcaactaaattatgcaaattctgcgggcgcgttaaCTTTGCCTAAATGAATGGCATCCAAAGCAATCAtttt                                                             <  1:947515‑M2/37‑1 (MQ=255)

                                  ttcatgaaaaataatctctctacaatacttcaactaaattatgcaaattctgcgggcgcgttaaCTTTGCCTAAATGAATGGCATCCAAAGCAATCAtttt                                                             <  3:484371‑M2/37‑1 (MQ=255)

                                     atgaaaaataatctctctacaatacttcaactaaattatgcaaattctgcgggcgcgttaaCTTTGCCTAAATGAATGGCATCCAAAGCAATCATTTTttc                                                          <  1:1890‑M2/40‑1 (MQ=255)

                                             taatctctctacaatacttcaactaaattatgcaaattctgcgggcgcgttaaCTTTGCCTAAATGAATGGCATCCAAAGCAATCATTTTTTCTTCGTTAg                                                  >  3:1513017‑M2/54‑101 (MQ=255)

                                                        caatacttcaactaaattatgcaaattctgcgggcgcgttaaCTTTGCCTAAATGAATGGCATCCAAAGCAATCATTTTTTCTTCGTTAGTCGGAATAACg                                       >  3:435358‑M2/43‑101 (MQ=255)

                                                                          atgcaaattctgcgggcgcgttaaCTTTGCCTAAATGAATGGCATCCAAAGCAATCATTTTTTCTTCGTTAGTCGGAATAACGGCACAAATGACACGCGCa                     <  3:1489769‑M2/77‑1 (MQ=255)

                                                                             caaattctgcgggcgcgttaaCTTTGCCTAAATGAATGGCATCCAAAGCAATCATTTTTTCTTCGTTAGTCGGAATAACGGCACAAATGACACGCGCAttt                  >  3:752321‑M2/22‑101 (MQ=255)

                                                                             caaattctgcgggcgcgttaaCTTTGCCTAAATGAATGGCATCCAAAGCAATCATTTTTTCTTCGTTAGTCGGAATAACGGCACAAATGACACGCGCAttt                  >  3:742151‑M2/22‑101 (MQ=255)

                                                                                attctgcgggcgcgttaaCTTTGCCTAAATGAATGGCATCCAAAGCAATCATTTTTTCTTCGTTAGTCGGAATAACGGCACAAATGACACGCGCATTTTCa               >  3:1150035‑M2/19‑101 (MQ=255)

                                                                                    tgcgggcgcgttaaCTTTGCCTAAATGAATGGCATCCAAAGCAATCATTTTTTCTTCGTTAGTCGGAATAACGGCACAAATGACACGCGCATTTTCACTgg           >  1:935178‑M2/15‑101 (MQ=255)

                                                                                     gcgggcgcgttaaCTTTGCCTAAATGAATGGCATCCAAAGCAATCATTTTTTCTTCGTTAGTCGGAATAACGGCACAAATGACACGCGCATTTTCACTGGa          >  1:637672‑M2/14‑101 (MQ=255)

                                                                                        ggcgcgttaaCTTTGCCTAAATGAATGGCATCCAAAGCAATCACTTTTTCTTCGTTAGTCGGAATAACGGCACAAATGACACGCGCATTTTCACTGGAAAc       <  3:306678‑M2/91‑1 (MQ=255)

                                                                                           gcgttaaCTTTGCCTAAATGAATGGCATCCAAAGCAATCATTTTTTCTTCGTTAGTCGGAATAACGGCACAAATGACACGCGCATTTTCACTGGAAACAAt    >  1:1199629‑M2/8‑101 (MQ=255)

                                                                                            cgttaaCTTTGCCTAAATGAATGGCATCCAAAGCAATCATTTTTTCTTCGTTAGTCGGAATAACGGCACAAATGACACGCGCATTTTCACTGGAAACAAtt   <  3:565228‑M2/95‑1 (MQ=255)

                                                                                             gttaaCTTTGCCTAAATGAATGGCATCCAAAGCAATCATTTTTTCTTCGTTAGTCGGAATAACGGCACAAATGACACGCGCATTTTCACTGGAAACAATTc  >  3:1441435‑M2/6‑101 (MQ=255)

CAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAACTTTGCCTAAATGAATGGCATCCAAAGCAATCATTTTTTCTTCGTTAGTCGGAATAACGGCACAAATGACACGCGCATTTTCACTGGAAACAATTC  >  NC_007779/3262235‑3262428

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 22 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`