BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,483,608	A→G	V77A (GTC→GCC)	valS ←	valine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,483,608	0	A	G	86.7%	33.2 / ‑1.3	15	V77A (GTC→GCC)	valS	valine‑‑tRNA ligase
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (5/8); total (6/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.20e-01

TTGTCGATGAAAGCTTCGCGGCCGTAGTCGTGACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTC  >  NC_000913/4483500‑4483741

ttGTCGATGAAAGCTTCGCGGCCGTAGTCGTGTCGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCCGCGTGGTCAGTACCGGCCTGCCACAGGGTGTTTTTGCCCTGCATgcg                                                                                                        <  2:394689/140‑1 (MQ=255)

ttGTCGATGAAAGCTTCGCGGCCGTAGTCGTGACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGGCCTGGCACAGGGTGTTTTTGCCCTGCATgcg                                                                                                        <  2:470858/140‑1 (MQ=255)

ttGTCGATGAAAGCTCCGCGGCAGTGGTCTTGACGGGGTTTACCTTCTCCTGCGGCAATCTTGCGCTCAACGACCTTCTGGGTAGCGATCCCGGCGTGGTCAGTACCGGCCTGCCACAGGGTGTTTTTGCCCTGCATgcg                                                                                                        <  1:218710/140‑1 (MQ=255)

 tGTCGATGAAAGCTTCGCGGCCGTAGTCGTGACGGGTTTTACCTTCTCCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCAGCGTGGTCTGTACCGGCCTGCCACACGGTGTTTTTGCCCTGCATgcgc                                                                                                       >  2:218710/1‑140 (MQ=255)

     gATGAAAGCTTCGCGGCCGTAGTCGTGACGGGTTTTACCTTCTTCTGCGGCAATCTTGCG‑TCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGGCCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATa                                                                                                  <  2:86518/140‑1 (MQ=255)

              ttCGCGGCCGTAGTCGTGACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCAGCCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCa                                                                                          <  1:459722/140‑1 (MQ=255)

                  cgGCCGTAGTCGTGACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCTGGCGTGGTCAGTACCGGCCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGt                                                                                      <  1:513713/140‑1 (MQ=255)

                  cgGCCGTAGTCGTGACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGGCCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGt                                                                                      <  1:424424/140‑1 (MQ=255)

                                    gTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGGCCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTg                                                                    >  1:556488/1‑140 (MQ=255)

                                     ttttACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTgg                                                                   >  1:437181/1‑140 (MQ=255)

                                               ttctGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGGCCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGAc                                                         >  2:100537/1‑140 (MQ=255)

                                                tctGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGGCCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGAcc                                                        <  1:253956/140‑1 (MQ=255)

                                                    cGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGGCCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGCTAGCGGATAATGGTATACATGATGGGTTGCTGGAAGGCGTGACCCata                                                    >  2:233193/1‑140 (MQ=255)

                                                                              tGGGTAGCGATCCCGGCGTGGTCAGTACCGGCCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTcggcgg                          <  2:666942/140‑1 (MQ=255)

                                                                                 gTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTcggcggcgg                       <  2:437181/140‑1 (MQ=255)

                                                                                                  gTCAGTACCGGCCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAAc      >  1:531826/1‑140 (MQ=255)

                                                                                                      ggtccGGCCTGCCACAGTGTTTATTTGCCCTGCATCCGCAGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTc  <  2:190932/137‑1 (MQ=255)

TTGTCGATGAAAGCTTCGCGGCCGTAGTCGTGACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTC  >  NC_000913/4483500‑4483741

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