BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	865,305	T→G	10.5%	S207S (TCT→TCG)	est →	carboxylesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	865,305	0	T	G	10.5%	98.2 / 8.1	38	S207S (TCT→TCG)	est	carboxylesterase
Reads supporting (aligned to +/- strand): ref base T (11/23); new base G (2/2); total (13/25)
Fisher's exact test for biased strand distribution p-value = 5.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

TTAAAAGAACTAAGTGAAGCATTAGATACGATTAAAGAGCAAGTTGATGAAGTGTTGGATCCTATTTTAGTGATTCAAGCAGAAAACGACAATATGATTGATCCACAATCCGCAAATTATATATATGACCATGTAGATTCTGATGACAAAAATATCAAGTGGTACAGTGAATCTGGACATGTTATTACGATTGATAAAGAGAAAGAACAAGTATTTGAAGATATTTATCAATTTTTAGAGTCATTAGACTGGTCAGAATAAAAAGAGATTTTAACAT  >  CP000730/865165‑865441

ttAAAAGAACTAAGTGAAGCATTAGATACGATTAAAGAGCAAGTTGATGAAGTGTTGGATCCTATTTTAGTGATTCAAGCAGAAAACGACAATATGATTGATCCACAATCCGCAAATTATATATATGACCATGTAGATTCt                                                                                                                                          <  1:411630/141‑1 (MQ=255)

ttAAAAGAACTAAGTGAAGCATTAGATACGATTAAAGAGCAAGTTGATGAAGTGTTGGATCCTATTTTAGTGATTCAAGCAGAAAACGACAATATGATTGATCCACAATCCGCAAATTATATATATGACCATGTAGATTCt                                                                                                                                          >  2:73922/1‑141 (MQ=255)

    aaGAACTAAGTGAAGCATTAGATACGATTAAAGAGCAAGTTGATGAAGTGTTGGATCCTATTTTAGTGATTCAAGCAGAAAACGACAATATGATTGATCCACAATCCGCAAATTATATATATGACCATGTAGATTCtgatg                                                                                                                                      >  2:319935/1‑141 (MQ=255)

       aaCTAAGTGAAGCATTAGATACGATTAAAGAGCAAGTTGATGAAGTGTTGGATCCTATTTTAGTGATTCAAGCAGAAAACGACAATATGATTGATCCACAATCCGCAAATTATATATATGACCATGTAGATTCTGATGACa                                                                                                                                   >  2:83001/1‑141 (MQ=255)

           aaGTGAAGCATTAGATACGATTAAAGAGCAAGTTGATGAAGTGTTGGATCCTATTTTAGTGATTCAAGCAGAAAACGACAATATGATTGATCCAAACTCCGCAAATTATATATATGACCATGTAGATTCTGATGACaaaaa                                                                                                                               <  1:366825/141‑1 (MQ=255)

             gTGAAGCATTAGATACGATTAAAGAGCAAGTTGATGAAGTGTTGGATCCTATTTTAGTGATTCAAGCAGAAAACGACAATATGATTGATCCACAATCCGCAAATTATATATATGACCATGTAGATTCTGATGACAAAAata                                                                                                                             >  2:78505/1‑141 (MQ=255)

             gTGAAGCATTAGATACGATTAAAGAGCAAGTTGATGAAGTGTTGGATCCTATTTTAGTGATTCAAGCAGAAAACGACAATATGATTGATCCACAATCCGCAAATTATATATATGACCATGTAGATTCTGATGACAAAAata                                                                                                                             >  2:361525/1‑141 (MQ=255)

              tGAAGCATTAGATACGATTAAAGAGCAAGTTGATGAAGTGTTGGATCCTATTTTAGTGATTCAAGCAGAAAACGACAATATGATTGATCCACAATCCGCAAATTATATATATGACCATGTAGATTCTGATGACAAAAatat                                                                                                                            >  2:111131/1‑141 (MQ=255)

                        gATACGATTAAAGAGCAAGTTGATGAAGTGTTGGATCCTATTTTAGTGATTCAAGCAGAAAACGACAATATGATTGATCCACAATCCGCAAATTATATATATGACCATGTAGATTCTGATGACAAAAATATCAAGTGGTAc                                                                                                                  <  1:245274/141‑1 (MQ=255)

                        gATACGATTAAAGAGCAAGTTGATGAAGTGTTGGATCCTATTTTAGTGATTCAAGCAGAAAACGACAATATGATTGATCCACAATCCGCAAATTATATATATGACCATGTAGATTCTGATGACAAAAATATCAAGTGGTAc                                                                                                                  <  2:20293/141‑1 (MQ=255)

                         aTACGATTAAAGAGCAAGTTGATGAAGTGTTGGATCCTATTTTAGTGATTCAAGCAGAAAACGACAATATGATTGATCCACAATCCGTATATTATATATATGACCATGTAGATTCTGATGACAAAAATATCAAGTGGTACa                                                                                                                 <  1:239103/141‑1 (MQ=255)

                         aTACGATTAAAGAGCAAGTTGATGAAGTGTTGGATCCTATTTTAGTGATTCAAGCAGAAAACGACAATATGATTGATCCACAATCCGCAAATTATATATATGACCATGTAGATTCTGATGACAAAAATATCAAGTGGTACa                                                                                                                 <  1:85129/141‑1 (MQ=255)

                          tACGATTAAAGAGCAAGTTGATGAAGTGTTGGATCCTATTTTAGTGATTCAGGCAGAAAACGACAATATGATTGATCCACAATCCGCAAATTATATATATGACCATGTAGATTCTGATGACAAAAATATCAAGTGGTACAg                                                                                                                <  2:417202/141‑1 (MQ=255)

                          tACGATTAAAGAGCAAGTTGATGAAGTGTTGGATCCTATTTTAGTGATTCAAGCAGAAAACGACACTATGATTGTTCCACAATCCGCAAATTATATATATGACCCTGTAGATTCTGATGACAAAAATTTCAAGTGGTACAg                                                                                                                <  2:157122/141‑1 (MQ=255)

                          tACGATTAAAGAGCAAGTTGATGAAGTGTTGGATCCTATTTTAGTGATTCAAGCAGAAAACGACAATATGATTGATCCACAATCCGCAAATTATATATATGACCATGTAGATTCTGATGACAAAAATATCAAGTGGTACAg                                                                                                                <  1:82535/141‑1 (MQ=255)

                            cGATTAAAGAGCAAGTTGATGAAGTGTTGGATCCTATTTTAGTGATTCAAGCAGAAAACGACAATATGATTGATCCACAATCCGCAAATTATATATATGACCATGTAGATTCTGATGACAAAAATATCAAGTGGTACAGTg                                                                                                              >  1:238533/1‑141 (MQ=255)

                              aTTAAAGAGCAAGTTGATGAAGTGTTGGATCCTATTTTAGTGATTCAAGCAGAAAACGACAATATGATTGATCCACAATCCGCAAATTATATATATGACCATGTAGATTCTGATGACAAAAATATCAAGTGGTACAGTGaa                                                                                                            <  2:101214/141‑1 (MQ=255)

                                        aaGTTGATGAAGTGTTGGATCCTATTTTAGTGATTCAAGCAGAAAACGACAATATGATTGATCCACAATCCGCAAATTATATATATGACCATGTAGATTCTGATGACAAAAATATCAAGTGGTACAGTGAATCTGGACATg                                                                                                  <  1:105363/141‑1 (MQ=255)

                                           ttGATGAAGTGTTGGATCCTATTTTAGTGATTCAAGCAGAAAACGACAATATGATTGATCCACAATCCGCAAATTATATATATGACCATGTAGATTCTGATGACAAAAATATCAAGTGGTACAGTGAATCTGGACATGtta                                                                                               >  1:154132/1‑141 (MQ=255)

                                            tgatgaAGTGTAGGATCCCATTGTAGTGATTCAAGCAGAAAACGACAATATGATTGATCCACAATACGCAAATTATATATATGACCATGTAGATTCTGATGACAAAAATATCAAGTGGTATAGTGAAACTTGACATGttat                                                                                              >  2:490059/1‑141 (MQ=255)

                                                gaAGTGTTGGATCCTATTTTAGTGATTCAAGCAGAAAACGACAATATGATTGATCCACAATCCGCAAATTATATATATGACCATGTAGATTCTGATGACAAAAATATCAAGTGGTACAGTGAATCTGGACATGTTATTACg                                                                                          <  2:275514/141‑1 (MQ=255)

                                                  aGTGTTGGATCCTATTTTAGTGATTCAAGCAGAAAACGACAATATGATTGATCCACAATCCGCAAATTATATATATGACCATGTAGATTCTGATGACAAAAATATCAAGTGGTACAGTGAATCTGGACATGTTATTACGAt                                                                                        <  1:78505/141‑1 (MQ=255)

                                                   gtgaaggATCCAATTTAAGTGAATCAAGCAGAAAACGACAATATGATTGATCCACAATCCGCAAATTACATATATGACCATGTAGATTCTGATGACAAAAAAATCAAGTGGTACAGTGAATCTCCACTTGTTAATACGAtt                                                                                       >  2:26750/6‑141 (MQ=255)

                                                                    aGTGATTCAAGCAGAAAACGACAATATGATTGATCCACAATCCGCAAATTATATATATGACCATGTAGATTCTGATGACAAAAATATCAAGTGGTACAGTGAATCTGGACATGTTATTACGATTGATAAAGAGAAAGAACa                                                                      >  1:84235/1‑141 (MQ=255)

                                                                     gTGATTCAAGCAGAAAACGACAATATGATTGATCCACAATCCGCAAATTATATATATGACCATGTAGATTCTGATGACAAAAATATCAAGTGGTACAGTGAATCTGGACATGTTATTACGATTGATAAAGAGAAAGGACaa                                                                     <  1:82623/141‑1 (MQ=255)

                                                                       gATTCAAGCAGAAAACGACAATATGATTGATCCACAATCCGCAAATTATATATATGACCATGTAGATTCTGATGACAAAAATATCAAGTGGTACAGTGAATCTGGACATGTTATTACGATTGATAAAGAGAAAGAACAAGt                                                                   <  1:25616/141‑1 (MQ=255)

                                                                        aTTCAAGCAGAAAACGACAATATGATTGATCCACAATCCGCAAATTATATATATGACCATGTAGATTCTGATGACAAAAATATCAAGTGGTACAGTGAATCTGGACATGTTATTACGATTGATAAAGAGAAAGAACAAGTa                                                                  <  1:269121/141‑1 (MQ=255)

                                                                           cAAGCAGAAAACGACAATATGATTGATCCACAATCCGCAAATTATATATATGACCATGTAGATTCTGATGACAAAAATATCAAGTGGTACAGTGAATCTGGACATGTTATTACGATTGATAAAGAGAAAGAACAAGTAttt                                                               >  2:104504/1‑141 (MQ=255)

                                                                             aGCAGAAAACGACAATATGATTGATCCACAATCCGCAAATTATATATATGACCATGTAGATTCTGATGACAAAAATATCAAGTGGTACAGTGAATCTGGACATGTTATTACGATTGATAAAGAGAAAGAACAAGTATTTGa                                                             >  1:172363/1‑141 (MQ=255)

                                                                                    aaCGACAATATGATTGATCCACAATCCGCAAATTATATATATGACCATGTAGATTCTGATGACAAAAATATCAAGTGGTACAGTGAATCTGGACATGTTATTACGATTGATAAAGAGAAAGAACAAGTATTTGAAGATAtt                                                      <  1:175580/141‑1 (MQ=255)

                                                                                            tatGATTGATCCACAATCCGCAAATTATATATATGACCATGTAGATTCTGATGACAAAAATATCAAGTGGTACAGTGAATCTGGACATGTTATTACGATTGATAAAGAGAAAGAACAAGTATTTGAAGATATTTATCAAtt                                              <  1:344218/141‑1 (MQ=255)

                                                                                                     tCCACAATCCGCAAATTATATATATGACCATGTAGATTCTGATGACAAAAATATCAAGTGGTACAGTGAATCTGGACATGTTATTACGATTGATAAAGATAAAGAACAAGTATTTGAAGATATTTATCAATTTTTAGAGTc                                     <  1:439602/141‑1 (MQ=255)

                                                                                                          aaTCCGCAAATTATATATATGACCATGTAGATTCGGATGACAAAAATATCAAGTGGTACAGTGAATCTGGACATGTTATTACGATTGATAAAGAGAAAGAACAAGTATTTGAAGATATTTATCAATTTTTAGAGTCATTAg                                >  2:112779/1‑141 (MQ=255)

                                                                                                          aaTCCGCAAATTATATATATGACCATGTAGATTCGGATGACAAAAATATCAAGTGGTACAGTGAATCTGGACATGTTATTACGATTGATAAAGAGAAAGAACAAGTATTTGAAGATATTTATCAATTTTTAGAGTCATTAg                                >  2:112780/1‑141 (MQ=255)

                                                                                                            tCCGCAAATTATATATATGACCATGTAGATTCTGATGACAAAAATATCAAGTGGTACAGTGAATCTGGACATGTTATTACGATTGATAAAGAGAAAGAACAAGTATTTGAAGATATTTATCAATTTTTAGAGTCATTAGAc                              <  1:73922/141‑1 (MQ=255)

                                                                                                                  aaTTATATATATGACCATGTAGATTCTGATGACAAAAATATCAAGTGGTACAGTGAATCTGGACATGTTATTACGATTGATAAAGAGAAAGAACAAGTATTTGAAGATATTTATCAATTTTTAGAGTCATTAGACTGGTCa                        <  1:218837/141‑1 (MQ=255)

                                                                                                                             tGACCATGTAGATTCTGATGACAAAAATATCAAGTGGTACAGTGAATCTGGACATGTTATTACGATTGATAAAGAGAAAGAACAAGTATTTGAAGATATTTATCAATTTTTAGAGTCATTAGACTGGTCAGAATAAAAaga             <  1:151218/141‑1 (MQ=255)

                                                                                                                             tGACCATGTAGATTCGGATGACAAAAATATCAAGTGGTACAGTGAATCTGGACATGTTATTACGATTGATAAAGAGAAAGAACAAGTATTTGAAGATATTTATCAATTTTTAGAGTCATTAGACTGGTCAGAATAAAAaga             <  1:112779/141‑1 (MQ=255)

                                                                                                                             tGACCATGTAGATTCGGATGACAAAAATATCAAGTGGTACAGTGAATCTGGACATGTTATTACGATTGATAAAGAGAAAGAACAAGTATTTGAAGATATTTATCAATTTTTAGAGTCATTAGACTGGTCAGAATAAAAaga             <  1:112780/141‑1 (MQ=255)

                                                                                                                                       gATTCTGATGACAAAAATATCAAGTGGTACAGTGAATCTGGACATGTTATTACGATTGATAAAGAGAAAGAACAAGTATTTGAAGATATTTATCAATTTTTAGAGTCATTAGACTGGTCAGAATAAAAAGAGATTTTAACa   <  1:286742/141‑1 (MQ=255)

                                                                                                                                        aTTCTGATGACAAAAATATCAAGTGGTACAGTGAATCTGGACATGTTATTACGATTGATAAAGAGAAAGAACAAGTATTTGAAGATATTTATCAATTTTTAGAGTCATTAGACTGGTCAGAATAAAAAGAGATTTTAACAt  <  1:148393/141‑1 (MQ=255)

TTAAAAGAACTAAGTGAAGCATTAGATACGATTAAAGAGCAAGTTGATGAAGTGTTGGATCCTATTTTAGTGATTCAAGCAGAAAACGACAATATGATTGATCCACAATCCGCAAATTATATATATGACCATGTAGATTCTGATGACAAAAATATCAAGTGGTACAGTGAATCTGGACATGTTATTACGATTGATAAAGAGAAAGAACAAGTATTTGAAGATATTTATCAATTTTTAGAGTCATTAGACTGGTCAGAATAAAAAGAGATTTTAACAT  >  CP000730/865165‑865441

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