BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,085,483	A→T	I83F (ATC→TTC)	gshB →	glutathione synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,085,483	0	A	T	94.7%	135.7 / ‑2.4	38	I83F (ATC→TTC)	gshB	glutathione synthase
Reads supporting (aligned to +/- strand): ref base A (1/1); new base T (23/13); total (24/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GGAGATGGGCGATCTGTATCTGATCAATGGTGAAGCCCGCGCCCATACCCGCACGCTGAACGTGAAGCAGAACTACGAAGAGTGGTTTTCGTTCGTCGGTGAACAGGATCTGCCGCTGGCCGATCTCGATGTGATCCTGATGCGTAAAGACCCGCCGTTTGATACCGAGTTTATCTACGCGACCTATATTCTGGAACGTGCCGAAGAGAAAGGGACGCTGATCGTTAACAAGCCGCAGAGCCTGCGCGACTGTAACGAGAAACTGTT  >  NZ_CP009273/3085350‑3085616

ggAGATGGGCGATCTGTATCTGATCAATGGTGAAGCCCGCGCCCATACCCGCACGCTGAACGTGAAGCAGAACTACGAAGAGTGGTTTTCGTTCGTCGGTGAACAGGATCTGCCGCTGGCCGATCTCGATGTGTTCCt                                                                                                                                   >  1:583808/1‑138 (MQ=255)

  agaTGGGCGATCTGTATCTGATCAATGGTGAAGCCCGCGCCCATACCCGCACGCTGAACGTGAAGCAGAACTACGAAGAGTGGTTTTCGTTCGTCGGTGAACAGGATCTGCCGCTGGCCGATCTCGATGTGTTCCTGAt                                                                                                                                >  2:490072/1‑139 (MQ=255)

       ggCGATCTGTATCTGATCAATGGTGAAGCCCGCGCCCATACCCGCACGCTGAACGTGAAGCAGAACTACGAAGAGTGGTTTTCGTTCGTCGGTGAACAGGATCTGCCGCTGGCCGATCTCGATGTGTTCCTGATGCGt                                                                                                                            >  1:539247/1‑138 (MQ=255)

       ggCGATCTGTATCTGATCAATGGTGAAGCCCGCGCCCATACCCGCACGCTGAACGTGAAGCAGAACTACGAAGAGTGGTTTTCGTTCGTCGGTGAACAGGATCTGCCGCTGGCCGATCTCGATGTGTTCCTGATGCGt                                                                                                                            >  1:55952/1‑138 (MQ=255)

       ggCGATCTGTATCTGATCAATGGTGAAGCCCGCGCCCATACCCGCACGCTGAACGTGAAGCAGAACTACGAAGAGTGGTTTTCGTTCGTCGGTGAACAGGATCTGCCGCTGGCCGATCTCGATGTGTTCCTGATGCGt                                                                                                                            >  1:513984/1‑138 (MQ=255)

       ggCGATCTGTATCTGATCAATGGTGAAGCCCGCGCCCATACCCGCACGCTGAACGTGAAGCAGAACTACGAAGAGTGGTTTTCGTTCGTCGGTGAACAGGATCTGCCGCTGGCCGATCTCGATGTGTTCCTGATGCGt                                                                                                                            >  1:36316/1‑138 (MQ=255)

                tATCTGATCAATGGTGAAGCCCGCGCCCATACCCGCACGCTGAACGTGAAGCAGAACTACGAAGAGTGGTTTTCGTTCGTCGGTGAACAGGATCTGCCGCTGGCCGATCTCGATGTGTTCCTGATGCGTAAAGACccg                                                                                                                   <  2:285231/138‑1 (MQ=255)

                     gATCAATGGTGAAGCCCGCGCCCATACCCGCACGCTGAACGTGAAGCAGAACTACGAAGAGTGGTTTTCGTTCGTCGGTGAACAGGATCTGCCGCTGGCCGATCTCGATGTGTTCCTGATGCGTAAAGACCCGCCGttt                                                                                                             >  2:796977/1‑139 (MQ=255)

                      aTCAATGGTGAAGCCCGCGCCCATACCCGCACGCTGAACGTGAAGCAGAACTACGAAGAGTGGTTTTCGTTCGTCGGTGAACAGGATCTGCCGCTGGCCGATCTCGATGTGTTCCTGATGCGTAAAGACCCGCCGttt                                                                                                             >  1:847675/1‑138 (MQ=255)

                           tGGTGAAGCCCGCGCCCATACCCGCACGCTGAACGTGAAGCAGAACTACGAAGAGTGGTTTTCGTTCGTCGGTGAACAGGATCTGCCGCTGGCCGATCTCGATGTGTTCCTGATGCGTAAAGACCCGCCGTTTGATAcc                                                                                                       >  2:361819/1‑139 (MQ=255)

                                     cgcgCCCATACCCGCACGCTGAACGTGAAGCAGAACTACGAAGAGTGGTTTTCGTTCGTCGGTGAACAGGATCTGCCGCTGGCCGATCTCGATGTGTTCCTGATGCGTAAAGACCCGCCGTTTGATACCGAGTTTATc                                                                                              >  1:658320/1‑138 (MQ=255)

                                     cgcgCCCATACCCGCACGCTGAACGTGAAGCAGAACTACGAAGAGTGGTTTTCGTTCGTCGGTGAACAGGATCTGCCGCTGGCCGATCTCGATGTGTTCCTGATGCGTAAAGACCCGCCGTTTGATACCGAGTTTATc                                                                                              >  1:671271/1‑138 (MQ=255)

                                            aTACCCGCACGCTGAACGTGAAGCAGAACTACGAAGAGTGTTTTTCGTTCGTCGGTGAACAGGATCTGCCGCTGGCCGATCTCGATGTGTTCCTGATGCGTa                                                                                                                           <  1:154372/102‑1 (MQ=255)

                                            aTACCCGCACGCTGAACGTGAAGCAGAACTACGAAGAGTGGTTTTCGTTCGTCGGTGAACAGGATCTGCCGCTGGCCGATCTCGATGTGTTCCTGATGCGTa                                                                                                                           >  2:154372/1‑102 (MQ=255)

                                            aTACCCGCACGCTGAACGTGAAGCAGAACTACGAAGAGTGGTTTTCGTTCGTCGGTGAACAGGATCTGCCGCTGGCCGATCTCGATGTGTTCCTGATGCGTAAAGACCCGCCGTTTGATACCGAGTTTATCTACGCGAc                                                                                      <  1:244049/139‑1 (MQ=255)

                                            aTACCCGCACGCTGAACGTGAAGCAGAACTACGAAGAGTGGTTTTCGTTCGTCGGTGAACAGGATCTGCCGCTGGCCGATCTCGATGTGTTCCTGATGCGTAAAGACCCGCCGTTTGATACCGAGTTTATCTACGCGAc                                                                                      >  1:680708/1‑139 (MQ=255)

                                               cccGCACGCTGAACGTGAAGCAGAACTACGAAGAGTGGTTTTCGTTCGTCGGTGAACAGGATCTGCCGCTGGCCGATCTCGATGTTTTCCTGATGCTTAAAGACCCGCCGTTTGATACCGAGTTTATCTACGCGACCt                                                                                    <  2:36316/138‑1 (MQ=255)

                                               cccGCACGCTGAACGTGAAGCAGAACTACGAAGAGTGGTTTTCGTTCGTCGGTGAACAGGATCTGCCGCTGGCCGATCTCGATGTGTTCCTGATGCGTAAAGACCCGCCGTTTGATACCGAGTTTATCTACGCGACCta                                                                                   <  2:658320/139‑1 (MQ=255)

                                               cccGCACGCTGAACGTGAAGCAGAACTACGAAGAGTGGTTTTCGTTCGTCGGTGAACAGGATCTGCCGCTGGCCGATCTCGATGTGTTCCTGATGCGTAAAGACCCGCCGTTTGATACCGAGTTTATCTACGCGACCt                                                                                    <  1:653344/138‑1 (MQ=255)

                                               cccGCACGCTGAACGGGAAGCAGAACTACGAAGAGTGGTTTTCGTTCGTCGGTGAACAGGATCTGCCGCTGGCCGATCTCGATGTGTTCCTGATGCGTAAAGACCCGCCGTTTGATACCGAGTTTATCTACGCGACCt                                                                                    <  2:55952/138‑1 (MQ=255)

                                                   cACGCTGAACGTGAAGCAGAACTACGAAGAGTGGTTTTCGTTCGTCGGTGAACAGGATCTGCCGCTGGCCGATCTCGATGTGTTCCTGATGCGTAAAGACCCGCCGTTTGATACCGAGTTTATCTACGCGACCTATAtt                                                                               >  2:35661/1‑139 (MQ=255)

                                                        tGAACGTGAAGCAGAACTACGAAGAGTGGTTTTCGTTCGTCGGTGAACAGGATCTGCCGCTGGCCGATCTCGATGTGTTCCTGATGCGTAAAGACCCGCCGTTTGATACCGAGTTTATCTACGCGACCTATATTCTGGa                                                                          <  2:847675/139‑1 (MQ=255)

                                                                aaGCAGAACTACGAAGAGTGGTTTTCGTTCGTCGGTGAACAGGATCTGCCGCTGGCCGATCTCGATGTGTTCCTGATGCGTAAAGACCCGCCGTTTGATACCGAGTTTATCTACGCGACCTATATTCTGGAACGTGCCg                                                                  <  1:35661/139‑1 (MQ=255)

                                                                                agTGGTTTTCGTTCGTCGGTGAACAGGATCTGCCGCTGGCCGATCTCGATGTGTTCCTGATGCGTAAAGACCCGCCGTTTGATACCGAGTTTATCTACGCGACCTATATTCTGGAACGTGCCGAAGAGAACGGGACGCt                                                  >  2:779586/1‑139 (MQ=255)

                                                                                agTGGTTTTCGTTCGTCGGTGAACAGGATCTGCCGCTGGCCGATCTCGATGTGTTCCTGATGCGTAAAGACCCGCCGTTTGATACCGAGTTTATCTACGCGACCTATATTCTGGAACGTGCCGAAGAGAAAg                                                         >  1:383935/1‑132 (MQ=255)

                                                                                agTGGTTTTCGTTCGTCGGTGAACAGGATCTGCCGCTGGCCGATCTCGATGTGTTCCTGATGCGTAAAGACCCGCCGTTTGATACCGAGTTTATCTACGCGACCTATATTCTGGAACGTGCCGAAGAGAAAg                                                         <  2:383935/132‑1 (MQ=255)

                                                                                agTGGTTTTCGTTCGTCGGTGAACAGGATCTGCCGCTGGCCGATCTCGATGTGTTCCTGATGCGTAAAGACCCGCCGTTTGATACCGAGTTTATCTACGCGACCTATATTCTGGAACGTGCCGAAGAGAAAGGGACGCt                                                  >  1:385521/1‑139 (MQ=255)

                                                                                agTGGTTTTCGTTCGTCGGTGAACAGGATCTGCCGCTGGCCGATCTCGATGTGTTCCTGATGCGTAAAGACCCGCCGTTTGATACCGAGTTTATCTACGCGACCTATATTCTGGAACGTGCCGAAGAGAAAGGGACGCt                                                  >  2:58101/1‑139 (MQ=255)

                                                                                  tGGTTTTCGTTCGTCGGTGAACAGGATCTGCCGCTGGCCGATCTCGATGTGTTCCTGATGCGTAAAGACCCGCCGTTTGATACCGAGTTTATCTACGCGACCTATATTCTGGAACGTGCCGAAGAGAAAGGGCCGCTGa                                                >  2:430974/1‑139 (MQ=255)

                                                                                  tGGTTTTCGTTCGTCGGTGAACAGGATCTGCCGCTGGCCGATCTCGATGTGTTCCTGATGCGTAAAGACCCGCCGTTTGATACCGAGTTTATCTACGCGACCTATATTCTGGAACGTGCCGAAGAGAAAGGGACGCTGa                                                >  1:61300/1‑139 (MQ=255)

                                                                                               tcgGTGAACAGGATCTGCCGCTGGCCGATCTCGATGTGTTCCTGATGCGTAAAGACCCGCCGTTTGATACCGAGTTTATCTACGCGACCTATATTCTGGAACGTGCCGAAGAGAAAGGGACGCTGATCGTTAACAAGcc                                   >  1:130472/1‑139 (MQ=255)

                                                                                                         ggATCTGCCGCTGGCCGATCTCGATGTGTTCCTGATGCGAAAAGCCCCGCCGTTTGATACCGAGTTTATCTACGCGACCTATATTCTGGAACGTGCCGAAGAGAAAGGGACGCTGATCGTTAACAAGCCGCAGAGCCTg                         <  1:490072/139‑1 (MQ=255)

                                                                                                               gCCGCTGGCCGATCTCGATGTGTTCCTGATGCGTAAAGACCCGCCGTTTGATACCGAGTTTATCTACGCGACCTATATTCTGGAACGTGCCGAAGAGAAAGGGACGCTGATCGTTAACAAGCCGCAGAGCCTGCGCGAc                   >  2:68142/1‑139 (MQ=255)

                                                                                                                   cTGGCCAACTTCGATGTGTTCGTGATGCGAAAAGGCCCGCCGTTTGATACCGAGTTTATCTACGCGACCTATATTCTGGAACGTGCCGAAGAGAAAGGGACGCTGATCGTTAACAAGCCGCAGAGCCTGCGCGACTGTa               <  2:804091/139‑1 (MQ=255)

                                                                                                                               gATGTTTTCCTGAGGCGTAAAGACCCGCCGTTTGATACCGAGTTTATCTACGCGACCTATATTCTGGAACGTGCCGAAGAGAAAGGGACGCTGATCGTTAACAAGCCGCAGAGCCTGCGCGACTGTAACgaga         <  2:437255/133‑1 (MQ=255)

                                                                                                                                aTGTGTTCCTGATGCGTAAAGACCCGCCGTTTGATACCGAGTTTATCTACGCGACCTATATTCTGGAACGTGCCGAAGAGAAAGGGACGCTGATCGTTAACAAGCCGCAGAGCCTGCGCGACTGTAACGAGAAACTGtt  >  2:39196/1‑139 (MQ=255)

                                                                                                                                aTGTGATCCTGATGCGTAAAGACCCGCCGTTTGATACCGAGTTTATCTACGCGACCTATATTCTGGAACGTGCCGAAGAGAAAgg                                                        >  2:123798/1‑85 (MQ=255)

                                                                                                                                aTGTGATCCTGATGCGTAAAGACCCGCCGTTTGATACCGAGTTTATCTACGCGACCTATATTCTGGAACGTGCCGAAGAGAAAgg                                                        <  1:123798/85‑1 (MQ=255)

GGAGATGGGCGATCTGTATCTGATCAATGGTGAAGCCCGCGCCCATACCCGCACGCTGAACGTGAAGCAGAACTACGAAGAGTGGTTTTCGTTCGTCGGTGAACAGGATCTGCCGCTGGCCGATCTCGATGTGATCCTGATGCGTAAAGACCCGCCGTTTGATACCGAGTTTATCTACGCGACCTATATTCTGGAACGTGCCGAAGAGAAAGGGACGCTGATCGTTAACAAGCCGCAGAGCCTGCGCGACTGTAACGAGAAACTGTT  >  NZ_CP009273/3085350‑3085616

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