BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_pae_tpiA	264,098	G→T	L50L (CTC→CTA)	ykfH ←	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_pae_tpiA	264,098	0	G	T	97.2%	102.8 / NA	36	L50L (CTC→CTA)	ykfH	uncharacterized protein
Reads supporting (aligned to +/- strand): ref base G (0/1); new base T (14/21); total (14/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAATT  >  NC_000913_3_pae_tpiA/263954‑264242

ggCTCCCAGGCGCAGTTGAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTGTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGa                                                                                                                                              <  2:625761/149‑1 (MQ=255)

ggCTCACAGGCGCGGTGTAATCGTCCGCTGCAGGCGCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTGAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGa                                                                                                                                              <  2:462598/149‑1 (MQ=255)

                  aaTCTCCCGCTGCAGGCCCCTGGTAGGGTTGCGCATTACAGTTCTCCACTACTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGaa                                                                                                                            <  1:202944/149‑1 (MQ=255)

                    tctcCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTTTAGGTGATTTTCAGTTTCGGCTAATACTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGGTAGGACTCCGGCGATATCCGGAACg                                                                                                                          >  1:628895/1‑149 (MQ=255)

                    tctcCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACg                                                                                                                          >  1:433427/1‑149 (MQ=255)

                                           ggtttGCTCATTACAGTTCTCCGCTATTGTCAGTTGAGGATGCTGCTGATGAGGCAGGTATAGGGCTCATTGCTGTCCTGGCGCCGGTGGGCGTGTTCCGCTAGGAGTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGAc                                                                                                   <  2:380082/146‑1 (MQ=25)

                                                   cATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGt                                                                                           <  1:454528/149‑1 (MQ=255)

                                                   cATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGt                                                                                           <  1:609911/149‑1 (MQ=255)

                                                           ttCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTgg                                                                                   >  2:307136/1‑149 (MQ=255)

                                                             ctcCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACTTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTa                                                                                 >  2:220467/1‑149 (MQ=255)

                                                              tcCACTATTGTCAGTTCAGGGTGATGCTACTCCCAGATGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGGCGGTGTCATGACACACGCTACCCTGCCACTgcgat                                                                                >  2:304016/1‑145 (MQ=255)

                                                               ccACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTgctctt                                                                               >  2:537535/1‑144 (MQ=255)

                                                                  cTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGc                                                                            <  1:261014/149‑1 (MQ=255)

                                                                   tATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATc                                                                                                                                 >  1:190797/1‑95 (MQ=255)

                                                                   tATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATc                                                                                                                                 <  2:190797/95‑1 (MQ=255)

                                                                      tGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGt                                                                        >  1:168856/1‑149 (MQ=255)

                                                                       gTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTg                                                                       <  1:535163/149‑1 (MQ=255)

                                                                                 ggTGATGCTCATCAGGCAGGTATAGGGCCCATTGTGTTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGaa                                                             <  1:282704/149‑1 (MQ=255)

                                                                                 ggTGATGCTCATCAGGCAGGTAGAGTGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAAGGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGaa                                                             <  1:627459/149‑1 (MQ=255)

                                                                                  gTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGaaa                                                            <  1:307136/149‑1 (MQ=255)

                                                                                       gCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTc                                                       >  1:399263/1‑149 (MQ=255)

                                                                                         tcatcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTCCACAGATTGCACTTCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCgg                                                     <  1:470744/149‑1 (MQ=18)

                                                                                         tcatcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATCTCCGGAACGTCCCTGCCGGTGTAATGAAAGACGCTACCGTGCCACTGGTATTTGCCGGTCCCGTAGCGCAAGCTTCgg                                                     >  2:530795/1‑149 (MQ=255)

                                                                                           atcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGa                                                   >  2:293486/1‑149 (MQ=255)

                                                                                              aggcaggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTc                                                >  1:28783/1‑149 (MQ=255)

                                                                                                   ggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGAt                                           <  2:168856/149‑1 (MQ=255)

                                                                                                    ggagaGGGCCCATTGCGGGCCGGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATg                                          <  2:182264/145‑1 (MQ=255)

                                                                                                     tataGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACACTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGc                                         >  2:393302/1‑149 (MQ=255)

                                                                                                           gccccgggCGGTCGCGGGGGCGTGCGGCGTAAACCGTGAGGTCTCCGGCGAGTTCCGGAACGTCCCTGCCGGTGTGGTGACAGACGCTATCGCGCTCCTGGTATTGGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCgg                                   <  2:109254/142‑1 (MQ=17)

                                                                                                           gCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCgg                                   <  1:293486/149‑1 (MQ=255)

                                                                                                           gCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCgg                                   >  1:560714/1‑149 (MQ=255)

                                                                                                                 tGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATc                             <  2:28783/149‑1 (MQ=255)

                                                                                                                 gtggtgccTGGCGGCGTTCGGCGTATACCGATATGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATc                             <  2:322852/143‑1 (MQ=17)

                                                                                                                                 ttCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGtt             <  2:433427/149‑1 (MQ=255)

                                                                                                                                   cGGCGTAGACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTTCCACTGGTATTTGCCGGTGAAGTAGCGACAGATGCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGtttc           <  2:486861/149‑2 (MQ=17)

                                                                                                                                            ccGCTAGGACTCGGGCGATATACGGAACGTCCCTGCCGGTGTAATGACAGACGCTACTGTGCCACTGGTATTTGCAGGTCATGTCTCGAAAGGTTCGACGTTCACGTTCGTGGCGGTATATCGTCATGGCCCTGCGTTTACTGATAAtt  <  2:137367/149‑1 (MQ=9)

GGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAATT  >  NC_000913_3_pae_tpiA/263954‑264242

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