BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	71,584	T→C	G78G (GGT→GGC)	yabI →	DedA family protein YabI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	71,584	0	T	C	90.9%	58.9 / ‑4.7	22	G78G (GGT→GGC)	yabI	DedA family protein YabI
Reads supporting (aligned to +/- strand): ref base T (2/0); new base C (6/14); total (8/14)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GGAACACTTTATTACCCAATCCACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTATTGGCTGCTTGATGGGCGACTGGAT‑TTCTTTCTGGCTGGGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGT  >  NC_000913/71365‑71810

ggAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGAT‑CTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTg                                                                                                                                                                                                       <  1:166219/249‑1 (MQ=255)

ggAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGAT‑CTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGt                                                                                                                                                                                                     <  1:44481/251‑1 (MQ=255)

              cccAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGGTTGGGAGCGGCGGGTTTAGGTTTTGGCCCGCCCGGCTGGCCGGGATTTTCGGGTTTTTGCTGGGCGAATGGGTCTCCCTTTTGGTGGGGTGGGGGTTTTAAAAA‑CCGTTGCCTCGCTGGTCAttttctaagaa                                                                                                                                                                                        >  1:130313/1‑242 (MQ=255)

                             aTTCATTGGTGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGAT‑CTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCattact                                                                                                                                                                         <  2:262906/250‑1 (MQ=255)

                                     aTGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGAT‑CTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAAc                                                                                                                                                                <  1:75519/251‑1 (MQ=255)

                                       ggCGGTGGTGGTGGTGGCCTTTCAGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGGGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGGTTCTGGCACGGCTGGCTGGCAGGGGTTGTCGGGTGGCTGCTGGGCGACAGGAT‑CTCCTTCTGGGTGGGCTGGGGTTTTTAAAAACCGTTGCATTGCTGGGCATTTTTGAAGAAAAAAAAAAGAATACTCGATaaaaatg                                                                                                                                                              >  1:211933/1‑248 (MQ=255)

                                           gtggtgTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGAT‑CTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGGCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACa                                                                                                                                                          <  1:164887/251‑1 (MQ=255)

                                                  tggtggCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGAT‑CTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGtt                                                                                                                                                    <  2:89535/250‑1 (MQ=255)

                                                   ggtggCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGAATCTACCAGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGCTTGGCAGCGGCGAGTTAAGTTTCTGGAACGCCTGGCTGGCAGGGGATGTCGGGTGTCTGCTGGGCGACTGGGT‑CTCCTTCTGGCTGGGCTGGCGGTTTAAAAAGCCGCTGCATCGACGGTCATTTCTGGAGAAAAAACAAGCCTTACTCGGTAAAACAGAACATGAGTTGc                                                                                                                                                  >  1:44097/1‑251 (MQ=255)

                                                              tGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGAT‑CTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACGCGATAAAACTGAACATGCGTTGCATCAACACAg                                                                                                                                        <  2:106320/250‑1 (MQ=255)

                                                                          cgctggTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGAT‑CTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCAtt                                                                                                                           <  1:184893/251‑1 (MQ=255)

                                                                                  ggTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGGTGGGGGCGCTGATTGGCAGCGGCGAGTTTAGTTTCTGGGACGCCTGGCTGGCCGGGGTTGTCGGTTGGCTGCTGGGCGGCCGGAT‑CTCCTTCTGGCTGGGCTGGCGTTTTTAAAAGACGTTGCATCGCTGGGCATTTTTGTAGAAAAAACAAGCATTACTCGCTAAAACTGAACCTGCGTTGCATCAACACAGCCTGTTCACCATTCTtgtcgg                                                                                                                   >  2:111874/1‑251 (MQ=255)

                                                                                      tGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGAT‑CTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGgtcggtcg                                                                                                                <  1:132457/250‑1 (MQ=255)

                                                                                        gttCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGTATTGTCGGTTGTCTGCTGGTCGACTGGAT‑CTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAGACTGAACATGCGTTGCATCAACACAGCATGTTCTCCATTCTGgtcggtcgctt                                                                                                             <  2:80310/250‑4 (MQ=255)

                                                                                                                            gcgcTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGAT‑CTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGc                                                                         <  1:23455/251‑1 (MQ=255)

                                                                                                                                 gATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGAT‑CTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGTCGGGaa                                                                    <  1:188347/251‑1 (MQ=255)

                                                                                                                                 gATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGAT‑CTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGaa                                                                    <  2:90662/251‑1 (MQ=255)

                                                                                                                                           ggcgAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGAT‑CTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCggggaagttggggct                                                          >  2:2690/1‑239 (MQ=255)

                                                                                                                                                        tCTGGCACGCCTGGCTGGCAGGGATTGTCGGATGTCTGCTGGGCGACTGGAT‑CTCCTTCTGGCTGGGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGCCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTCGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGACAATGGAGGCGGGAATGCTGGATCTGCCGGTGGcaaa                                              >  1:158514/1‑247 (MQ=255)

                                                                                                                                                                        ggcagggattgtcggttgtctgctGGGCGACTGGAT‑CTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAat                              <  1:114148/227‑1 (MQ=255)

                                                                                                                                                                                     ggttgtctgctGGGCGACTGGAT‑CTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGctgcct                >  1:235612/11‑251 (MQ=255)

                                                                                                                                                                                                   cGACTGGAT‑CTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGt  >  1:66891/1‑251 (MQ=255)

GGAACACTTTATTACCCAATCCACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTATTGGCTGCTTGATGGGCGACTGGAT‑TTCTTTCTGGCTGGGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGT  >  NC_000913/71365‑71810

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