BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,634,019	T→C	I125V (ATT→GTT)	nohA ←	NohA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,634,019	0	T	C	79.7%	26.6 / 22.8	59	I125V (ATT→GTT)	nohA	NohA
Reads supporting (aligned to +/- strand): ref base T (6/6); new base C (22/25); total (28/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.03e-01

ATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCAC‑CGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGG  >  NC_000913/1633928‑1634111

atatCCCGTTTCAGGAAATCAACATGTCGGTTTTACAGTTCC‑GACAAAGCCGCCCAACTCTACAAGGGGAACCCGTCCACAATACTGG‑AAttt                                                                                            >  1:3885603/1‑93 (MQ=0)

     ccttttCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCAC‑CGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCAcc                                                                                        <  1:1295662/90‑1 (MQ=1)

                    aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAAC‑CGACAGAGGAACTCCATCACGAATGCTGGCTACTTCTCCGGCCCTCCGCGACAg                                                                         >  1:5013349/4‑93 (MQ=0)

                    aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAAC‑CGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAg                                                                         >  1:2214147/4‑93 (MQ=0)

                    aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAAC‑CGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAg                                                                         >  2:7136342/4‑93 (MQ=0)

                    aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAAC‑CGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAg                                                                         >  2:1079211/4‑93 (MQ=0)

                    aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAAC‑CGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAg                                                                         >  1:3908292/4‑93 (MQ=0)

                    aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAAC‑CGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAg                                                                         >  1:3139736/4‑93 (MQ=0)

                    aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAAC‑CGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAg                                                                         >  2:4081035/4‑93 (MQ=0)

                                    agctCCGGGAAGCGCCGCTGAAC‑CGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAAt                                                         <  1:3632982/90‑1 (MQ=0)

                                    agctCCGGGAAGCGCCGCTGAAC‑CGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAAt                                                         <  1:5577826/90‑1 (MQ=0)

                                    agctCCGGGAAGCGCCGCTGAAC‑CGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAAt                                                         >  1:8094508/4‑93 (MQ=0)

                                    agctCCGGGAAGCGCCGCTGAAC‑CGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAAt                                                         >  2:4812722/4‑93 (MQ=0)

                                    agctCCGGGAAGCGCCGCTGAAC‑CGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAAt                                                         >  2:8323188/4‑93 (MQ=0)

                                    agctCCGGGAAGCCCCGCTGAAC‑CGACAGAGGAACTCCATCAAGAATGCTGGCTATTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAAt                                                         <  2:2922592/90‑1 (MQ=0)

                                            gaagcgccgcTGAAC‑CGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTc                                                 >  2:3838494/5‑93 (MQ=0)

                                            gaagcgccgcTGAAC‑CGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTc                                                 >  1:1171938/5‑93 (MQ=0)

                                            gaagcgccgcTGAAC‑CGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTc                                                 >  2:1554131/5‑93 (MQ=0)

                                            gaagcgccgcTGAAC‑CGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTc                                                 <  2:5708663/89‑1 (MQ=0)

                                            gaagcgccgcTGAAC‑CGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTc                                                 <  2:7734907/89‑1 (MQ=0)

                                            gaagcgccgcTGAAC‑CGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTc                                                 >  2:8946125/5‑93 (MQ=0)

                                            gaagcgccgcTGAAC‑CGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTc                                                 <  1:5021187/89‑1 (MQ=0)

                                            gaagcgccgcTGAAC‑CGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTc                                                 <  1:1803465/89‑1 (MQ=0)

                                             aagcgccgcTGAAC‑CGACATAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTcc                                                <  1:7750695/90‑1 (MQ=0)

                                             aagcgccgcTGAAC‑CGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTcc                                                >  2:4317570/4‑93 (MQ=0)

                                             aagcgccgcTGAAC‑CGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTcc                                                <  2:463259/90‑1 (MQ=0)

                                             aagcgccgcTGAAC‑CGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTcc                                                <  2:2876776/90‑1 (MQ=0)

                                             aagcgccgcTGAAC‑CGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTcc                                                <  2:4714319/90‑1 (MQ=0)

                                                      tgaaC‑CGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCa                                       <  1:945715/90‑1 (MQ=0)

                                                                           atcaagaATTCTGGCTACTACTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCAttttt                   <  2:2636927/89‑1 (MQ=0)

                                                                           atcaagaATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCAttttt                   <  2:2058666/89‑1 (MQ=0)

                                                                           atcaagaATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCAttttt                   <  2:7320715/89‑1 (MQ=0)

                                                                            tctagaATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCTGAATGCGGTGTCCACCACCTCAGCGTATTCTTTTGCAtttttt                  >  2:8965858/4‑92 (MQ=0)

                                                                            tcaagaATGCTGGCTACTTCTCCGGCCATCCGCGCCAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGTATCTTTTGCAtttttt                  >  1:4014866/4‑92 (MQ=0)

                                                                            tcaagaATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCAtttttt                  <  1:2659324/90‑2 (MQ=0)

                                                                            tcaagaATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCAtttttt                  <  1:2659325/90‑2 (MQ=0)

                                                                            tcaagaATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCAtttttt                  <  2:5583767/90‑2 (MQ=0)

                                                                            tcaagaATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCAtttttt                  >  2:4751983/4‑92 (MQ=0)

                                                                            tcaagaATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCAtttttt                  >  2:7517864/4‑92 (MQ=0)

                                                                            tcaagaATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCAtttttt                  <  2:9878503/90‑2 (MQ=0)

                                                                            tcaagaATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCAtttttt                  >  1:8561269/4‑92 (MQ=0)

                                                                            tcaagaATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCCCAGCGGGATCTTTTGCAtttttt                  >  2:7225435/4‑92 (MQ=0)

                                                                                  atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTcc            <  1:647544/90‑1 (MQ=0)

                                                                                  atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTcc            <  2:97096/90‑1 (MQ=0)

                                                                                  atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTcc            <  1:1375327/90‑1 (MQ=0)

                                                                                  atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTcc            <  2:813397/90‑1 (MQ=0)

                                                                                  atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTcc            >  1:8722093/4‑93 (MQ=0)

                                                                                  atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTcc            <  1:9522514/90‑1 (MQ=0)

                                                                                  atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTcc            <  2:592915/90‑1 (MQ=0)

                                                                                  atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTcc            <  1:2804173/90‑1 (MQ=0)

                                                                                  atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTcc            <  2:2177975/90‑1 (MQ=0)

                                                                                  atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTcc            <  1:3960162/90‑1 (MQ=0)

                                                                                        gCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCg      <  1:3825496/93‑1 (MQ=31)

                                                                                         cAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGt     >  2:745874/1‑93 (MQ=32)

                                                                                          caTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTc    <  1:4799219/92‑1 (MQ=18)

                                                                                          aaTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTc    <  2:2158454/93‑1 (MQ=33)

                                                                                          aaTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTc    >  2:9711953/1‑93 (MQ=33)

                                                                                           aTTTCACGTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGTGATTTTTCAGCTCCTGTGCGTCg   >  2:6168764/1‑93 (MQ=255)

                                                                                           aTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCg   >  1:4485965/1‑93 (MQ=33)

                                                                                           aTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCg   >  2:1813895/1‑93 (MQ=33)

                                                                                           aTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCg   <  1:8091176/93‑1 (MQ=33)

                                                                                           aTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCg   >  1:3838734/1‑93 (MQ=33)

                                                                                            tttCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGAGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCgg  >  1:1019359/1‑93 (MQ=18)

ATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCAC‑CGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGG  >  NC_000913/1633928‑1634111

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