BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	3,104,817	A→G	12.1%	M76T (ATG→ACG)	leuS ←	leucyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	3,104,817	0	A	G	12.1%	78.2 / 5.3	33	M76T (ATG→ACG)	leuS	leucyl‑tRNA synthetase
Reads supporting (aligned to +/- strand): ref base A (12/17); new base G (2/2); total (14/19)
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.28e-01

CGATCCCAGTCATATGAGAAGCCAAGCGCTTGAATTTGGCGGCGGAAGTTATCAATATTCTGCTTCGTAAACACAGCGGGGTCGTTCCCTGTGTCAAGCGCGTACTGTTCAGCTGGCAGGCCGAATGCGTCCCAGCCCATTGGATGAAGGACATCATAGCCCTGCATGCGCTTCATGCGGGACAGAATATCCGTAGCTGTGTATCCTTCAGGATGGCCGACGTGCAACCCAGCTCCAGACGGATAAGGAAACATGTCCAGCGCGTAA  >  BS168_reference_genome/3104679‑3104945

cgATCCCAGTCATATGAGAAGCCAAGCGCTTGAATTTGGCGGCGGAAGTTATCAATATTCTGCTTCGTAAACACAGCGGGGTCGTTCCCTGTGTCAAGCGCGTACTGTTCAGCTGGCAGGCCGAATGCGTCCCAGCCCa                                                                                                                                  <  4:128118/139‑1 (MQ=255)

  aTCCCAGTCATATGAGAAGCCAAGCGCTTGAATTTGGCGGCGGAAGTTATCAATATTCTGCTTCGTAAACACAGCGGGGTCGTTCCCTGTGTCAAGCGCGTACTGTTCAGCTGGCAGGCCGAATGCGTCCCAGCCCAtt                                                                                                                                <  1:138074/139‑1 (MQ=255)

       aGTCATATGAGAAGCCAAGCGCTTGAATTTGGCGGCGGAAGTTATCAATATTCTGCTTCGTAAACACAGCGGGGTCGTTCCCTGTGTCAAGCGCGTACTGTTCAGCTGGCAGGCCGAATGCGTCCCAGCCCATTGGATg                                                                                                                           <  6:173021/139‑1 (MQ=255)

           atatGAGAAGCCAAGCGCTTGAATTTGGCGGCGGAAGTTATCAATATTCTGCTTCGTAAACACAGCGGGGTCGTTCCCTGTGTCAAGCGCGTACTGTTCAGCTGGCAGGCCGAATGCGTCCCAGCCCATTGGATGAAg                                                                                                                        <  3:245838/138‑1 (MQ=255)

           atatGAGAAGCCAAGCGCTTGAATTTGGCGGCGGAAGTTATCAATATTCTGCTTCGTAAACACAGCGGGGTCGTTCCCTGTGTCAAGCGCGTACTGTTCAGCTGGCAGGCCGAATGCGTCCCAGCCCATTGGATGAAg                                                                                                                        <  7:205602/138‑1 (MQ=255)

            tatGAGAAGCCAAGCGGTTGAATTTGGCGGCGGAAGTTATGAATATTCTGCTTCGTAAACACAGCGGGGTCGTTCCCTGTGTGAAGCGCGTACTGTTCAGCTGGCAGGACGATTGCGTCCCAGCCCATTGGATGAAGGa                                                                                                                      <  4:136026/139‑1 (MQ=255)

            tatGAGAAGCCAAGCGCTTGAATTTGGCGGCGGAAGTTATCAATATTCTGCTTCGTAAACACAGCGGGGTCGTTCCCTGTGTCAAGCGCGTACTGTTCAGCTGGCAGGCCGAATGCGTCCCAGCCCATTGGATGAAGGa                                                                                                                      <  2:190469/139‑1 (MQ=255)

            tatGAGAAGCCAAGCGCTTGAATTTGGCGGCGGAAGTTATCAATATTCTGCTTCGTAAACACAGCGGGGTCGTTCCCTGTGTCAAGCGCGTACTGTTCAGCTGGCAGGCCGAATGCGTCCCAGCCCATTGGATGAAGGa                                                                                                                      <  7:2820/139‑1 (MQ=255)

                 gaAGCCAAGCGCTTGAATTTGGCGGCGGAAGTTATCAATATTCTGCTTCGTAAACACAGCGGGGTCGTTCCCTGTGTCAAGCGCGTACTGTTCAGCTGGCAGGCCGAATGCGTCCCAGCCCATTGGATGAAGGAcatca                                                                                                                 >  4:200504/1‑139 (MQ=255)

                   aGCCAAGCGCTTGAATTTGGCGGCGGAAGTTATCAATATTCTGCTTCGTAAACACAGCGGGGTCGTTCCCTGTGTCAAGCGCGTACTGTTCAGCTGGCAGGCCGAATGCGTCCCAGCCCATTGGATGAAGGAcatca                                                                                                                 >  2:228460/1‑137 (MQ=255)

                    gCCAAGCGCATGAAGTTGGCGGCGGAAGTTATCAATATTCTGCTTCGTAACCACAGCGGGTTCGTTCCCTGTGTCAAGCGCGTACTGTTCAGCTGGCAGGCCGAATCCGGCCCAGCCCATAGGATGAAGGACATCATAg                                                                                                              <  8:22813/139‑1 (MQ=255)

                     ccAAGCGCTTGAATTTGGCGGCGGAAGTTATCAATATTCTGCTTCGTAAACACAGCGGGGTCGTTCCCTGTGTCAAGCGCGTACTGTTCAGCTGGCAGGCCGAATGCGTCCCAGCCCATTGGATGAAGGACATCATAGc                                                                                                             <  3:200504/139‑1 (MQ=255)

                     ccAAGCGCTTGAATTTGGCGGCGGAAGTTATCAATATTCTGCTTCGTAAACACAGCGGGGTCGTTCCCTGTGTCAAGCGCGTACTGTTCAGCTGGCAGGCCGAATGCGTCCCAGCCCATTGGATGAAGGACATCATAGc                                                                                                             >  7:193294/1‑139 (MQ=255)

                               gAATTTGGCGGCGGAAGTTATCAATATTCTGCTTCGTAAACACAGCGGGGTCGTTCCCTGTGTCAAGCGCGTACTGTTCAGCTGGCAGGCCGAATGCGTCCCAGCCCATTGGATGAAg                                                                                                                        <  7:30333/118‑1 (MQ=255)

                               gAATTTGGCGGCGGAAGTTATCAATATTCTGCTTCGTAAACACAGCGGGGTCGTTCCCGGTGTCAAGCGCGTACTGTTCAGCTGGCAGGCCGAATGCGTCCCAGCCCATTGGATGAAg                                                                                                                        >  8:30333/1‑118 (MQ=255)

                                               acagtCAATATTCTGCTTCGTAAACAGAGCGGGGTCGTTCCCTGTGTAAAGCGCGTACTGTTCAGCTGGCAGGCCGAATGCGTCCCAGCCCATTGGATGAAGGTCATCATAGCCCTGCATGCGCTTCATGAGGGACAGa                                                                                   <  7:73981/135‑1 (MQ=255)

                                                     aaTATTCTGCTTCGTAAACACAGCGGGGTCGTTCCCTGTGTCAAGCGCGTACTGTTCAGCTGGCAGGCCGAATGCGTCCCAGCCCATTGGATGAAGGACATCATAGCCCTGCATGCGCTTCATGCGGGACAGAATATcc                                                                             >  3:274496/1‑139 (MQ=255)

                                                      atatTCTGCTTCGTAAACACAGCGGGGTCGTTCCCTGTGTCAAGCGCGTACTGTTCAGCTGGCAGGCCGAATGCGTCCCAGCCCATTGGATGAAGGACATCATAGCCCTGCATGCGCTTCATGCGGGACAGAATATCCg                                                                            <  1:228460/139‑1 (MQ=255)

                                                      atatTCTGCTTCGTAAACACAGCGGGGTCGTTCCCTGTGTCAAGCGCGTACTGTTCAGCTGGCAGGCCGAATGCGTCCCAGCCCATTGGATGAAGGACATCATAGCCCTGCATGCGCTTCATGCGGGACAGAATATCCg                                                                            <  3:81682/139‑1 (MQ=255)

                                                          tCTGCTTCGTAAACACAGCGGGGTCGTTCCCTGTGTCAAGCGCGTACTGTTCAGCTGGCAGGCCGAATGCGTCCCAGCCCGTTGGATGAAGGACATCATAGCCCTGCATGCGCTTCATGCGGGAc                                                                                      <  7:284116/125‑1 (MQ=255)

                                                          tCTGCTTCGTAAACACAGCGGGGTCGTTCCCTGTGTCAAGCGCGTACTGTTCAGCTGGCAGGCCGAATGCGTCCCAGCCCGTTGGATGAAGGACATCATAGCCCTGCATGCGCTTCATGCGGGAc                                                                                      <  7:287952/125‑1 (MQ=255)

                                                          tCTGCTTCGTAAACACAGCGGGGTCGTTCCCTGTGTCAAGCGCGTACTGTTCAGCTGGCAGGCCGAATGCGTCCCAGCCCGTTGGATGAAGGACATCATAGCCCTGCATGCGCTTCATGCGGGAc                                                                                      >  8:284116/1‑125 (MQ=255)

                                                          tCTGCTTCGTAAACACAGCGGGGTCGTTCCCTGTGTCAAGCGCGTACTGTTCAGCTGGCAGGCCGAATGCGTCCCAGCCCGTTGGATGAAGGACATCATAGCCCTGCATGCGCTTCATGCGGGAc                                                                                      >  8:287952/1‑125 (MQ=255)

                                                              cTTCGTAAACACAGCGGGGTCGTTCCCTGTGTCAAGCGCGTACTGTTCAGCTGGCAGGCCGAATGCGTCCCAGCCCATTGGATGAAGGACATCATAGCCCTGCATGCGCTTCATGCGGGACAGAATATCCGTAGCt                                                                       <  1:75306/136‑1 (MQ=255)

                                                                            cGGGGTCGTTACCTGTGTCAAGCGCGTACTGTTCAGCTGGCACGCCGAATGCGTCCCAGCCCATTGGAAGAAGGACATCATGGCCCTCCATGCGCGTCCTGCGGGCCAGAAAATCCGTAGCTGTGTATCCTTCAGGATg                                                      >  6:25871/1‑139 (MQ=255)

                                                                                gTCGTTCCCTGTGTCAAGCGCGTACTGTTCAGCTGGCAGGCCGAATGCGTCCCAGCCCATTGGATGAAGGACATCATAGCCCTGCATGCGCTTCATGCGGGTCAGAATATCCGTAGCTGTGTATCCTTCAGGATGGCCg                                                  <  5:34428/139‑1 (MQ=255)

                                                                                      cccTGTGTCAAGCGCGTACTGTTCAGCTGGCAGGCCGAATGCGTCCCAGCCCATTGGATGAAGGACATCATAGCCCTGCATGCGCTTCATGCGGGACAGAATATCCGTAGCTGTGTATCCTTCAGGATGGCCGACGTg                                             <  8:193294/138‑1 (MQ=255)

                                                                                         tgtgTCAAGCGCGTACTGTTCAGCTGGCAGGCCGAATGCGTCCCAGCCCATTGGATGAAGGACATCATAGCCCTGCATGCGCTTCATGCGGGACAGAATATCCGTAGCTGTGTATCCTTCAGGATGGCCGACGTGCacc                                         >  8:7601/1‑137 (MQ=255)

                                                                                                  cgcgTACTGTTCAGCTGGCAGGCCGAATGCGTCCCAGCCCATTGGATGAAGGACATCATAGCCCTGCATGCGCTTCATGCGGGACAGAATATCCGTAGCTGTGTATCCTTCAGGATGGCCGACGTGCAACCCAGCTCCa                                <  6:91593/139‑1 (MQ=255)

                                                                                                                cTGGCAGGCCGAATGCGTCCCAGCCCATTGGATGAAGGACATCATAGCCCTGCATGCGCTTCATGCGGGACAGAATATCCGTAGCTGTGTATCCTTCAGGATGGCCGACGTGCAACCCAGCTCCAGACGGATAAGGaaa                  >  3:173905/1‑139 (MQ=255)

                                                                                                                     aGGCCGAATGCGTCCCAGCCCATTGGATGAAGGACATCATAGCCCTGCCTGCGCGTCATGCGGGACAGAATATCCGTAGCGGTGTATCCTTCAGGATGGCCGACGTGCAACCCAGCTCCAGACGGATAAGGACACATGt             >  5:42214/1‑139 (MQ=255)

                                                                                                                     aGGCCGAATGCGTCCCAGCCCATTGGATGAAGGACATCATAGCCCTGCATGCGCTTCATGCGGGACAGAATATCCGTAGCTGTGTATCCTTCAGGATGGCCGACGTGCACCCCAGCTCCAGACGGATAAGGAAACATGt             >  4:163845/1‑139 (MQ=255)

                                                                                                                     aGGCCGAATGCGTCCCAGCCCATTGGATGAAGGACATCATAGCCCTGCATGCGCTTCATGCGGGACAGAATATCCGTAGCTGTGTATCCTTCAGGATGGCCGACGTGCAACCCCGCTCCAGACGGATAAGGAAACATGt             >  6:93193/1‑139 (MQ=255)

                                                                                                                                 tCCCAGCCCATTGGATGAAGGACATCATAGCCCTGCATGCGCTTCATGCGGGACAGAATATCCGTAGCTGTGTATCCTTCAGGATGGCCGACGTGCAACCCAGCTCCAGACGGATAAGGAAACATGTCCAGCGCGTaa  >  4:150312/1‑138 (MQ=255)

CGATCCCAGTCATATGAGAAGCCAAGCGCTTGAATTTGGCGGCGGAAGTTATCAATATTCTGCTTCGTAAACACAGCGGGGTCGTTCCCTGTGTCAAGCGCGTACTGTTCAGCTGGCAGGCCGAATGCGTCCCAGCCCATTGGATGAAGGACATCATAGCCCTGCATGCGCTTCATGCGGGACAGAATATCCGTAGCTGTGTATCCTTCAGGATGGCCGACGTGCAACCCAGCTCCAGACGGATAAGGAAACATGTCCAGCGCGTAA  >  BS168_reference_genome/3104679‑3104945

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: `‑`