BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,137,772	G→A	intergenic (‑115/+160)	yijF ← / ← gldA	DUF1287 domain‑containing protein YijF/L‑1,2‑propanediol dehydrogenase/glycerol dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,137,772	0	G	A	100.0%	67.6 / NA	20	intergenic (‑115/+160)	yijF/gldA	DUF1287 domain‑containing protein YijF/L‑1,2‑propanediol dehydrogenase/glycerol dehydrogenase
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (9/11); total (9/11)

GCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGTGAAGCTTTCATCTGGTTCGTCCGTGAAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATAAAAAATTCCTGAACAACGTGGGAGTTTCAGATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCC  >  NC_000913/4137539‑4138012

gCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGTGAAGCTTTCATCTGGTTCGTCCGTGAAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATCGTTCTGTAGAGAAAATGTCCAGAAGCGGATAAAAAATTCCTGAACAACGTGGGAGTTTCAAATTACTGGCACagag                                                                                                                                                                                                                                   <  2:178690/249‑1 (MQ=255)

             tGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGTGAAGCTTTCATCTGGTTCGTCCGTGAAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATAAAAAATTCCTGAACAACGTGGGAGTTTCAAATTACTGGCACAGAGTTCCGGATGCGGAg                                                                                                                                                                                                                     <  2:356815/250‑1 (MQ=255)

                           aTTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGTGAAGCTTTCATCTGGTTCGTCCGTGAAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATAAAAAATTCCTGAACAACGTGGGAGTTTCAAATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCg                                                                                                                                                                                                      <  2:185460/251‑1 (MQ=255)

                                    aCGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGTGAAGCGTTCATCTGGTTCGTCCGTGAAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATAAAAAATTCCTGAACAACGTGGGAGTTTCAAATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTCCaaa                                                                                                                                                                                             <  2:13952/251‑1 (MQ=255)

                                                    ggggAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGTGAAGCTTTCATCTGGTTCGTCCGTGAAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATAAAAAATTCCTGAACAACGTGGGAGTTTCAAATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAat                                                                                                                                                                             <  2:286712/250‑1 (MQ=255)

                                                              tAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGTGAAGCTTTCATCTGGTTCGTCCGTGAAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATAAAAAATTCCTGAACAACGTGGGAGTTTCAAATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCag                                                                                                                                                                    <  2:31735/250‑1 (MQ=255)

                                                                                                          gTGAAGCTTTCATCTGGTTCGTCCGTGAAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATAAAAAATTCCTGAACAACGTGGGAGTTTCAAATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCg                                                                                                                       >  2:114621/1‑251 (MQ=255)

                                                                                                                     aTCTGGTTCGTCCGTGAAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATAAAAAATTCCTGAACAACGTGGGAGTTTCAAATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTccc                                                                                                            >  2:70782/1‑251 (MQ=255)

                                                                                                                              gTCCGTGAAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATAAAAAATTCCTGAACAACGTGGGAGTTTCAAATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGcc                                                                                                   >  2:286167/1‑251 (MQ=255)

                                                                                                                                             agTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATAAAAAATTCCTGAACAACGTGGGAGTTTCAAATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGt                                                                                      >  1:263226/1‑249 (MQ=255)

                                                                                                                                             agTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATAAAAAATTCCTGAACAACGTGGGAGTTTCAAATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGt                                                                                      >  1:73235/1‑249 (MQ=255)

                                                                                                                                                  caGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATAAAAAATTCCTGAACAACGTGGGAGTTTCAAATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTAt                                                                               <  2:78895/251‑1 (MQ=255)

                                                                                                                                                               ccGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATAAAAAATTCCTGAACAACGTGGGAGTTTCAAATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGc                                                                   >  1:106979/1‑250 (MQ=255)

                                                                                                                                                                                                   aaGCGGATAAAAAATTCCTGAACAACGTGGGAGTTTCAAATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCagag                              <  1:163051/251‑1 (MQ=255)

                                                                                                                                                                                                                aTTCCTGAACAACGTGGGAGTTTCAAATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTg                 >  1:115565/1‑251 (MQ=255)

                                                                                                                                                                                                                           acGTGGGAGTTTCAAATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTcg      <  2:246646/251‑1 (MQ=255)

                                                                                                                                                                                                                           acGTGGGAGTTTCAAATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTcg      <  2:102100/251‑1 (MQ=255)

                                                                                                                                                                                                                            cGTGGGAGTTTCAAATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTcgc     >  1:9911/1‑251 (MQ=255)

                                                                                                                                                                                                                              tGGGAGTTTCAAATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCgcgc   >  2:210994/1‑251 (MQ=255)

                                                                                                                                                                                                                               gggAGTTTCAAATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGGAAACCTGATCTGGCGTCGcgcc  <  2:342318/251‑1 (MQ=255)

GCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGTGAAGCTTTCATCTGGTTCGTCCGTGAAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATAAAAAATTCCTGAACAACGTGGGAGTTTCAGATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCC  >  NC_000913/4137539‑4138012

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