BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,947,465	Δ1 bp	12.5%	intergenic (+2/‑32)	metZ → / → metW	tRNA‑Met/tRNA‑Met

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,947,465	0	T	.	12.5%	81.7 / 7.6	32	intergenic (+2/‑32)	metZ/metW	tRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand): ref base T (8/19); new base . (2/2); total (10/22)
Fisher's exact test for biased strand distribution p-value = 5.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCC  >  NC_000913/2947319‑2947566

aaaTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAAtt                                                                                                      >  5:233979/1‑148 (MQ=255)

  aTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTaaa                                                                                                   <  7:390857/149‑1 (MQ=255)

   tCCGTGATCTGGTTAAAAAAAAGGTTGCATGAAATCGATTTCGGAGTATAGTGCGCATCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCaaat                                                                                                  <  8:415630/149‑2 (MQ=18)

     cGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAtt                                                                                                <  6:233979/149‑1 (MQ=255)

         aTCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgat                                                                                            <  2:90380/149‑1 (MQ=255)

         aTCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgat                                                                                            <  3:224965/149‑1 (MQ=255)

          tCGGTGTAAAAAAAAGGTTGCATGAAAACGCGAGCGTAGTATAGTTCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgatg                                                                                           <  2:163930/149‑1 (MQ=255)

          tCGGGGTAAAAACAAGGTGGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGGGGAGCAGCCTGGTAGCTCGGCGGGCTCATCACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgatg                                                                                           <  1:286030/149‑1 (MQ=25)

          tCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgatg                                                                                           <  1:150436/149‑1 (MQ=255)

             gggTAAAAAAAAGGTTGCATGAAATCTCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTAGGGCTCTTAACCCGAAGGTCGTCGGTTCAATTCCGGCCCCCGCAACCAATTAAAATTTGATGAAg                                                                                        <  7:376852/149‑1 (MQ=21)

             gggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAg                                                                                        <  1:116799/149‑1 (MQ=255)

                   aaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAACCTCGTAACACCTCCGGCTAATAAACCAAAGGTCGTCCGTTAAAATACAGCCCCCGAAACCAATTAAAAATAAAAAAAGTAAAg                                                                                   >  2:279143/1‑148 (MQ=255)

                    aaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAttt                                                                                               <  8:144321/135‑1 (MQ=255)

                    aaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAttt                                                                                               >  7:144321/1‑135 (MQ=255)

                         ggTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACg                                                                            <  4:284530/149‑1 (MQ=255)

                                    aaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAATGCAGTACGGTGACGCgggt                                                                 >  5:202676/1‑148 (MQ=255)

                                     aCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGt                                                                <  2:421724/149‑1 (MQ=255)

                                     aCGCGAGCGGAGTATAGTGCGCATCCACGGAAGCGGGGTGGAGCAGCCTTGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCTGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGt                                                                <  4:229347/149‑1 (MQ=255)

                                         gAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTtgatg                                                                                           >  1:137172/1‑117 (MQ=255)

                                         gAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTtgatg                                                                                           <  2:137172/117‑1 (MQ=255)

                                          cgCGGAGTATAGTTCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGTCTCTTAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGagc                                                           <  1:279143/148‑1 (MQ=255)

                                          aGCGGAGTATAGGGCTAATTCACGGACGCTGGGTGGAGCAGCCTGGTAGCTCGTCGGTCTCATATCCCGAAGGTGTTCGGTTCAAATCATGCCCCCTCAACCAAATAAAATTTGATGAAGATAAGCAGTACCGTGACGCGGGGTGGagc                                                           <  6:63857/149‑1 (MQ=14)

                                                        cgcATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTaaaa                                                                                   <  1:392974/110‑2 (MQ=2)

                                                        cgcATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTaaaa                                                                                   >  2:392974/1‑109 (MQ=2)

                                                                                    cTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtga                                                                                             <  2:226049/73‑1 (MQ=255)

                                                                                    cTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtga                                                                                             >  1:226049/1‑73 (MQ=255)

                                                                                     tGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGa                                                             >  6:29763/1‑104 (MQ=255)

                                                                                     tGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGa                                                             <  5:29763/104‑1 (MQ=255)

                                                                                                                             aaaTCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGc                                                >  6:350833/1‑77 (MQ=255)

                                                                                                                                        ccGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTg                                                               >  5:353677/1‑51 (MQ=255)

                                                                                                                                        ccGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTg                                                               <  6:353677/51‑1 (MQ=255)

                                                                                                                                                aaTTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGccccc  >  7:281440/1‑104 (MQ=40)

                                                                                                                                                aaTTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGccccc  <  8:281440/104‑1 (MQ=40)

AAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCC  >  NC_000913/2947319‑2947566

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