BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,793,898	Δ1 bp	coding (212/552 nt)	btuE ←	glutathione peroxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,793,897	0	T	.	94.6%	146.3 / 2.6	38	coding (213/552 nt)	btuE	glutathione peroxidase
Reads supporting (aligned to +/- strand): ref base T (1/1); new base . (20/15); total (22/16)
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

GCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGG  >  NC_000913/1793751‑1794039

gcgACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCGATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCtgt                                                                                                                                              >  1:900843/1‑147 (MQ=255)

     tGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAg                                                                                                                                        >  2:924719/1‑149 (MQ=255)

      gCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTGATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGa                                                                                                                                       >  1:642928/1‑149 (MQ=255)

         gTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGGCCAGTAAGTTTTAATCTCTTCATCCGTGCCCGGTTCTT‑TTCCAGAAAc                                                                                                                                    >  2:257868/1‑149 (MQ=255)

          tCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCCGAAACt                                                                                                                                   >  2:796856/1‑149 (MQ=255)

                    ggCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACgg                                                                                                                         <  1:212956/149‑1 (MQ=255)

                       aatcaatTTTTGATACAGCGGATGGCGTCCTTCGCCATTCAGTTGAATGTTGCTGAACATCGGGAACGTCACCCCCCATGTGGTTGTACAGTTAGTTTTAATCTCTTCATAGCCCCCCGGGT‑ATGTTCCCGCACCGGGTTGCACGGGaa                                                                                                                      >  1:373389/1‑149 (MQ=255)

                       aatcaatTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGGGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGTTTGCACGGGaa                                                                                                                      >  1:130197/1‑149 (MQ=255)

                       aatcaatTTTTGATACAGCGGATAGCGTCCTTCGCCAGTAACTTCAATCTTGCTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGCGaa                                                                                                                      >  1:510949/1‑149 (MQ=255)

                         tcaatTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATc                                                                                                                    <  2:556528/149‑1 (MQ=255)

                           aatTGTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATccc                                                                                                                  <  1:283231/149‑1 (MQ=255)

                                   aTACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCAt                                                                                                          <  1:497129/149‑1 (MQ=255)

                                   aTACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCAt                                                                                                          >  1:501817/1‑149 (MQ=255)

                                    tACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATa                                                                                                         >  1:348989/1‑149 (MQ=255)

                                    tACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATa                                                                                                         <  1:373297/149‑1 (MQ=255)

                                       aGCGGAGGGCCTCCTTCGCCATTAACTTCAATCGTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATaaaa                                                                                                      <  2:91775/149‑1 (MQ=255)

                                          ggACGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCt                                                                                                   <  1:234613/149‑1 (MQ=255)

                                                     ttCGCCATTAACTTCAATCTTACTGAAAATCGGGAACGTCACCCGCCATGTGGTGGTACAGTAAGTGTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCCGAAACTGGTTGCACGGGACTCCCAGCACCCTAAAGCCGCGATCAAccac                                                                                        >  2:237587/1‑147 (MQ=255)

                                                           aTTAACTTCAAGTTTGCTGAACTTCGGGAACGTCACCCCCCCTGTGGTGGGCCAGTAATTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCttt                                                                                  <  2:900875/149‑1 (MQ=255)

                                                               aCTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTGAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCCTTCTGa                                                                              >  1:351066/1‑149 (MQ=255)

                                                                cTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGGCCCAGGCTTTCTGaa                                                                             >  1:70653/1‑149 (MQ=255)

                                                                     aTCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATtc                                                                        <  2:485377/149‑1 (MQ=255)

                                                                          aCTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCaa                                                                   <  1:797770/149‑1 (MQ=255)

                                                                            tGAACACCGGGAACCGCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACt                                                                 <  1:408848/149‑1 (MQ=255)

                                                                                   cGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATACAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCata                                                          >  2:4787/1‑149 (MQ=255)

                                                                                       aaCGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTgc                                                      <  1:924719/149‑1 (MQ=255)

                                                                                               cccccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATAGTTCGGCGTTaa                                              >  1:159991/1‑149 (MQ=255)

                                                                                                         gtggtACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACttt                                    >  1:860851/1‑149 (MQ=255)

                                                                                                            gtACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGCTCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAg                                 <  1:459947/149‑1 (MQ=255)

                                                                                                              aCAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCCCACTTTGAGGc                               >  1:807635/1‑149 (MQ=255)

                                                                                                                 gTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCGCAGCACCATAAAACCGCGAACGACCCAGGCTTACTGCATATTCTCGAACTGCCCATATTGCGGCGTTAGGCCACAATTTGAACCGAc                            >  2:255602/1‑149 (MQ=255)

                                                                                                                  tAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAACCCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACa                           >  2:392674/1‑149 (MQ=255)

                                                                                                                      ttttAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGa                       >  1:827198/1‑149 (MQ=255)

                                                                                                                       tttAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCCCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGAc                      >  1:457122/1‑149 (MQ=255)

                                                                                                                       tttAATCTCTGCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTGGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAACCCACACTTTAAGGAGACATTGaa                      >  2:181423/1‑148 (MQ=255)

                                                                                                                           aTCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATc                  <  1:353949/149‑1 (MQ=255)

                                                                                                                                  cctCGCTGCCCGGTTCTT‑TTCCAGAGACTGGTTGCCCGGGAATCCCAGCACCATAAAACCTCGATCGTCCCGGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGca           <  2:313842/147‑1 (MQ=255)

                                                                                                                                         gCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTAcc    <  2:351066/149‑1 (MQ=255)

                                                                                                                                           ccGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTCCCgg  >  1:510655/1‑149 (MQ=255)

GCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGG  >  NC_000913/1793751‑1794039

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