BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	225,499	T→A	intergenic (+42/‑1)	ileV → / → alaV	tRNA‑Ile/tRNA‑Ala

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	225,499	0	T	A	75.8%	71.8 / 20.0	33	intergenic (+42/‑1)	ileV/alaV	tRNA‑Ile/tRNA‑Ala
Reads supporting (aligned to +/- strand): ref base T (2/6); new base A (9/16); total (11/22)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

GTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTT  >  NC_000913/225396‑225637

cctttagctcagctggttagagcgcacccctgataagggtgaggtcggtggttcgagtccactaaggcccaccattatcaacacatgagtgttgatcccaTTAAGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGc                                                                                                        <  1:1900064‑M1/40‑1 (MQ=255)

 ctttagctcagctggttagagagcacccctgataagggtgaggtcggtggttcgagtccactaaggcccaccattatcaacacatgagtgttgatcccaTTAAGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCa                                                                                                       >  2:373424‑M1/100‑140 (MQ=255)

  tttagctcagctggttagagcgcacccctgataagggtgaggtcggtggttcgagtccactaaggcccaccattatcaacacatgagtgttgatcccaTTAAGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGcaa                                                                                                      <  2:1346434‑M1/42‑2 (MQ=255)

   ttagctcagctggttagagcgcacccctgataagggtgaggtcggtggttcgagtccactaaggcccgccattatcaacacatgagtgttgatcccaTTAAGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTgcaag                                                                                                     <  2:2169636‑M1/43‑3 (MQ=255)

   ttagctcagctggttagagcgcacccctgataagggtgaggtcggtggttcgagtccactaaggcccaccattatcaacacatgagtgttgatcccaTTAAGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTgcaag                                                                                                     >  2:1121051‑M1/98‑138 (MQ=255)

          agctggttagagcgcacccctgataagggtgaggtcggtggttcgagtccactaaggcccaccattatcaacacatgagtgttgatcccaTTAAGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCAAGCaggagg                                                                                              <  2:1384288‑M1/50‑1 (MQ=255)

              ggttagagcgcacccctgataagggtgaggtcggtggttcgagtccactaaggcccaccattatcaacacatgagtgttgatcccaTTAAGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCAAGCAGGAGgtcag                                                                                          >  1:1372127‑M1/87‑138 (MQ=255)

                    agcgcacccctgataagggtgaggtcggtggttcgagtccactaaggcccaccattatcaacacatgagtgttgatcccaTTAAGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCAAGCAGGAGGTCAGCGGTTc                                                                                    <  1:1376237‑M1/60‑1 (MQ=255)

                      cgcacccatgataagggtgaggtcggtggttcgagtccactaaggcccaccattatcaacacatgagtgttgatcccaTTAAGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCAAGCAGGAGGTCAGCGGTTCGa                                                                                  <  1:881194‑M1/62‑1 (MQ=255)

                       gcacccctgataagggtgaggtcggtggttcgagtccactaaggcccaccattatcaacacatgagtgttgatcccaTTAAGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCAAGCAGGAGGTCAGCGGTTCGAt                                                                                 >  2:435597‑M1/78‑140 (MQ=255)

                            cctgataagggtgaggtcggtggttcgagtccactaaggcccaccattatcaacacatgagtgttgatcccaTTAAGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCATTGCAAGCAGGAGGTCAGCGGTTCGATCCCGc                                                                            >  1:673176‑M1/73‑140 (MQ=255)

                                ataagggtgaggtcggtggttcgagtccactaaggcccaccattatcaacacatgagtgttgatcccaTTAAGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCAAGCAGGAGGTCAGCGGTTCGATcccgcttgg                                                                        <  1:540404‑M1/72‑5 (MQ=255)

                                 taatggtgaggtcggtggttcgagtccactaaggcccaccattatcaacacatgagtgttgatcccaTTAAGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCAAGCAGGATGTCAGCGGTTCGAtcccgcttggc                                                                       >  2:396223‑M1/68‑135 (MQ=255)

                                 taagggtgaggtcggtggttcgagtccactaaggcccaccattatcaacacatgagtgttgatcccaTTAAGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCAAGCAGGAGGTCAGCGGTTCGAtcccgcttggc                                                                       >  2:110482‑M1/68‑135 (MQ=255)

                                     ggtgaggtcggtggttcgagtccactaaggcccaccattatcaacacatgagtgttgatcccaTTAAGGGGCCATAGCTCAGCTGGGATAGCGCCTGCCTTGCAAGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTcca                                                                   <  2:922173‑M1/77‑1 (MQ=255)

                                     ggtgaggtcggtggttcgagtccactaaggcccaccattatcaacacatgagtgttgatcccaTTAAGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCAAGCAGGAGGTCAGCGGTTCGATcc                                                                               <  1:533366‑M1/65‑1 (MQ=255)

                                     ggtgaggtcggtggttcgagtccactaaggcccaccattatcaacacatgagtgttgatcccaTTAAGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCAAGCAGGAGGTCAGCGGTTCGATcc                                                                               >  2:533366‑M1/64‑128 (MQ=255)

                                     ggtgaggtcggtggttcgagtccactaaggcccaccattatcaacacatgagtgttgatcccaTTAAGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCAAGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTcca                                                                   <  1:524313‑M1/77‑1 (MQ=255)

                                       tgaggtcggtggttcgagtccactaaggcccaccattatcaacacatgagtgttgatcccaTTAAGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCAAGCAGGAGGTCAGCGGTTCGATCCCGCTTGGctccatt                                                                 <  2:508768‑M1/79‑3 (MQ=255)

                                       tgaggtcggtggttcgagtccactaaggcccaccattatcaacacatgagtgttgatcccaTTAAGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCAAGCAGGAGGTCAGCGGTTCGATCCCGCTTGGctccatt                                                                 <  1:1444867‑M1/79‑3 (MQ=255)

                                       tgaggtcggtggttcgagtccactaaggcccaccattatcaacacatgagtgttgatcccaTTAAGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCAAGCAGGAGGTCAGCGGTTCGATCCCGCTTGGctccatt                                                                 <  2:283561‑M1/79‑3 (MQ=255)

                                       tgaggtcggtggttcgagtccactaaggcccaccattatcaacacatgagtgttgatcccaTTAAGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCAAGCAGGAGGTCAGCGGTTAGATCCCGCTTGGctccatt                                                                 <  2:36239‑M1/79‑3 (MQ=255)

                                         aggtcggtggttcgagtccactaaggcccaccattatcaacacatgagtgttgatcccaTTAAGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCAAGCAGGAGGTCAGCGGTTCGATCCCGCTTggctccatttg                                                               <  2:1549939‑M1/81‑5 (MQ=255)

                                          ggtcggtggttcgagtccactaaggcccaccattatcaacacatgagtgttgatcccaTTAAGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCAAGCAGGAGGTCAGCGGTTCGATCCCGCTtggctccatttgt                                                              <  2:1114860‑M1/82‑6 (MQ=255)

                                            tcggtggttcgagtccactaaggcccaccattatcaacacatgagtgttgatcccaTTAAGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCAAGCAGGAGGTCAGCGGTTCGATCCCGcttggctccatttgttc                                                            >  1:1159044‑M1/57‑133 (MQ=255)

                                                      gCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGTACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGtt                                                  >  1:1830812/1‑140 (MQ=17)

                                                               aTTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATaaaaa                                         <  1:1810169/140‑1 (MQ=33)

                                                                     aCGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACtt                                   <  1:2316965/140‑1 (MQ=33)

                                                                                      gTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTCAGAGTGTACCTGCaaa                  >  1:1689869/1‑140 (MQ=32)

                                                                                          tAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATTCTTCAGAGTGTACCTGCAAAGGtt              <  1:1611121/140‑1 (MQ=17)

                                                                                                 aTGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCCCCATCTCTGTAGTGGTTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCg       <  2:1897241/140‑1 (MQ=17)

                                                                                                 aTGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCg       <  2:2161450/140‑1 (MQ=33)

                                                                                                      aTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGtt  <  2:2093185/140‑1 (MQ=33)

GTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTT  >  NC_000913/225396‑225637

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