BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	1,978,503	Δ776 bp		insB5–insA5	insB5, insA5

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1978503–1979267	1979278	12–776	11 [0]	[0] 11	insB5–insA5	insB5,insA5

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1978502	0 (0.000)	11 (1.050)	10/294	0.1	100%	intergenic (‑305/+16)	flhD/insB5	DNA‑binding transcriptional dual regulator FlhD/IS1 family protein InsB
*	?	NC_000913	1979279 =	0 (0.000)	11 (1.050)	10/294	0.1	100%	intergenic (‑64/‑474)	insA5/uspC	IS1 family protein InsA/universal stress protein C

ATGGCCATATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGG  >  NC_000913/1979034‑1979359

gccgggaagccccggtaaaaaataattagcattagaatagatgcgataaaatgcaataagcagaactacatttttggtttaatatgtccttacaaatagaactgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtata                                                                               >  1:117890‑M2/246‑249 (MQ=255)

                  aaaataattagcattagaatagttgcgataagctgcaataagcagaaccacctttttggtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTggg                                                           <  1:209116‑M2/24‑1 (MQ=255)

                                                  ctgcaataagcagaaccacctttttggtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCa                                                  >  1:181930‑M2/196‑228 (MQ=255)

                                                    gcaataagcagaaccacctttttggtttaatatgtccttacaaatagaaatgggtctttactcttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTAGAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacg                         <  1:175182‑M2/58‑1 (MQ=255)

                                                     caataagcagaaccacctttttggtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaccggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattccttaacataaaagtatgacttatccatttatgttaag                                                           >  1:241577‑M2/193‑195 (MQ=255)

                                                                          tggtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgccactgtactcgtcactacccgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattccttaacgaaaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTc                               >  1:148242‑M2/172‑223 (MQ=255)

                                                                            gtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtacgcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGTGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAgg  <  1:193312‑M2/81‑1 (MQ=255)

                                                                            gtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatgtttcctaaatcgacgcaactgtactcggcactacacgcacatacagcggagtggggctgcgattttcaataatgcgtgatgcagatcacacaaaagacgcaagtacttaacataaatgTATAAGTCAGACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGt                                    <  1:230021‑M2/47‑1 (MQ=255)

                                                                                                                                   ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaa                                         >  1:13862‑M2/115‑156 (MQ=255)

                                                                                                                                                                                                tcaatnntgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaa                                         <  2:13862‑M2/42‑1 (MQ=255)

                                                                                                                                                                                                                     cacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTc                               >  1:121850‑M2/33‑84 (MQ=255)

ATGGCCATATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGG  >  NC_000913/1979034‑1979359

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 16 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`