BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,176,209	C→A	D209Y (GAT→TAT)	cpdA ←	cAMP phosphodiesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,176,209	0	C	A	97.6%	101.9 / NA	42	D209Y (GAT→TAT)	cpdA	cAMP phosphodiesterase
Reads supporting (aligned to +/- strand): ref base C (1/0); new base A (18/22); total (20/22)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATG  >  NC_000913/3176075‑3176341

cGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCAGCGACCATTCCAGTCGAGATCGACCTcca                                                                                                                              >  2:81413/1‑142 (MQ=255)

 gggcgaTCAGCATTTAACCCGAGAGTGCGCCAGACGGGCCCGATGGTATCCAGCGGACAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCAACGGCGTTGCCAGCAGGCGGTGACCATTCCAGTAGAGATACAGCTCCtg                                                                                                                             <  2:185742/138‑1 (MQ=255)

 gggCCGTCAGCTTGTAGCTCGAGCGTACGCCAGCCGGGCGCGATGGTATCCAGCATCCTGTGTGGTCTGTGCGGCTTAAACGGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCtg                                                                                                                             <  2:276321/141‑1 (MQ=255)

  tGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGACGACGGCGTTGCCAGCAGGCGGCGACCAATCCAGCCGAGAAACAGCCCCaga                                                                                                                            >  2:479146/1‑143 (MQ=255)

      aTCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGATGGAACAGTGCGGCTGAAACTGCACACAGGCCGACGGCGTTGCCAGCAGGCGGCGACCATACCAGTCGAGATACAGCTCCTCAGGaa                                                                                                                        >  2:534920/1‑143 (MQ=255)

         aGCATTTAATGCGAGAGTACGCCTCTCGGGTGCGATGGTATCGAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAatat                                                                                                                     <  1:271375/143‑1 (MQ=255)

         aGCATGTAACTCGAGAGTACGCCAGCCTGGCGCGATGTTATCCAGCGTAATGTTGGAACAGTGCGGTTTAATTTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAatat                                                                                                                     <  1:196644/143‑1 (MQ=255)

          gCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCACCAGCCGTCGACCATTCCAGTCGAGATACAGCTCCTGGCGAATCTg                                                                                                                    >  2:578384/1‑143 (MQ=255)

              gTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCCAGATACAGCTCCTGATGAATATGACCg                                                                                                                >  2:19228/1‑143 (MQ=255)

               tAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGc                                                                                                               >  2:234135/1‑143 (MQ=255)

               tAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCGGATGAATATGACCGc                                                                                                               >  2:339673/1‑143 (MQ=255)

               tAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCAATCCAGTCGAGATACAGCTCCGGATGAAGACGACCGc                                                                                                               >  2:38392/1‑143 (MQ=255)

                 aCTCGAGAGTACGCCAGCCGGGCGCGATACAATTTAGCAAAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGATGCCAGCAGGCGGCGACCATACCAGCCGAGATACAGCCCCCGAGGAAGATGACCGccc                                                                                                             >  2:677541/1‑141 (MQ=255)

                      agagTACGCCAGCCTGGCTTTTTGTTTCCTAACGTAAATCTGTAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCaa                                                                                                        <  1:330037/143‑1 (MQ=255)

                      agagTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGTTTAAACTGCACACAGGTCGACGGGGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCaa                                                                                                        <  1:768382/143‑1 (MQ=255)

                                   cGGGCGCGATGGTATCCAGCGTAATGTTGGAACATTGTGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGc                                                                                           <  1:276986/143‑1 (MQ=255)

                                   cGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGc                                                                                           <  1:473623/143‑1 (MQ=255)

                                    gggCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGCCGAGATACAGCTCCTGATGAAAATGGCCGGACAGCAAGTATTTGACGTGCg                                                                                          >  2:804010/1‑143 (MQ=255)

                                       cgcgATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCATGATGAATATGACCGCACAACAAGTATATGACGTGCGGaa                                                                                       >  1:810546/1‑143 (MQ=255)

                                        gcgATGGTATCCAGCGTAAAGTTGGACCAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGaaa                                                                                      <  1:1606/143‑1 (MQ=255)

                                        gcgATGCTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCTACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGaaa                                                                                      <  1:1601/143‑1 (MQ=255)

                                            tGGTTTTCATCGTAAAGTTGGAACAGTTCGGCTTAAACTGTACACATGTCGACGGTGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTc                                                                                  <  1:144229/143‑1 (MQ=255)

                                            tGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCAGAAGAAAATGACCGCCCAGCAAGTATTAGACGTGCGGAAACtta                                                                                  >  2:93022/1‑142 (MQ=255)

                                                 tCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAg                                                                             >  1:326794/1‑143 (MQ=255)

                                                   cAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATTCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCa                                                                           >  2:394499/1‑143 (MQ=255)

                                                        tAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCCGTCGAGATACAGCTCCTGATGAATATGACCACACAGCAAGTATTTGACGTGCGGAAACTACGCCAGCACGGTa                                                                      >  1:116819/1‑143 (MQ=255)

                                                         acaGTTCGAACAGTGCGGCTTCAACTGCACACAGGTTGACGGTGTTGCCTGCAGGCGGCGTCCATTCCAGTCGAGATACAGCTCCGGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTAt                                                                     <  1:269146/141‑1 (MQ=255)

                                                          aaGTTGGAACAGTGCGGCTTAAACTTCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATc                                                                    <  1:547310/143‑1 (MQ=255)

                                                                 aaCTGTGCGGCTTAATCTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCGGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCg                                                             <  1:578384/143‑1 (MQ=255)

                                                                          gCTTAAACTGCACACACGCCGACGGCCTTGCGAGCAGGCCGCTCCCATTCCAGCCTCGATACATCTCCTTATCAATATAACCGCACAGCAATTATTTGACTCGCTCAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTa                                                    <  2:442944/143‑1 (MQ=255)

                                                                                             cgagtgtgTTGCCAGCGGCCTGCGACCATTCCAGTCGAGTTACAGCTCCTGATTAATATGACCTCACAGCACGTAGTTGCCGTGCGAAAATTTCTTCAGCACTGTATCTAGTTCTCTCGCGTTACGCAGACTGTGTTGAGCGa                                 <  2:247176/136‑1 (MQ=255)

                                                                                             cGACGGCGGTGGCAGCGGGCGGCGACGATTCCAGTCGAGAGACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTGCGCCCGCGTTACGCAGACTGTGTTGATCGa                                 <  2:452994/143‑1 (MQ=255)

                                                                                                   cGTTGCCAGCAGGCGGCGACCATTCCAGTCTATATACTGCTCCTGTTTAATATGTCCGTACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAAc                           <  1:259348/143‑1 (MQ=255)

                                                                                                         cagcatGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCTCACATCAAGTATTTGATTTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAAc                     <  1:719376/143‑1 (MQ=255)

                                                                                                         cagcagGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCTCAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTAGCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAAc                     <  1:212619/143‑1 (MQ=255)

                                                                                                                     aCCATCCCAGTCGATATACAGCTCCTTATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCg         <  1:801562/143‑1 (MQ=255)

                                                                                                                       cATTCTAGTTTAGATACAGCTCCTGATGATTATGACCGCACAGCAAGTATTTGACGTGCTGAAACTTCGCCTGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGga       <  2:50652/143‑1 (MQ=255)

                                                                                                                       cATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTGGATCGAGCCAACTACAACCCGCAGGTAGCGga       >  2:235371/1‑143 (MQ=255)

                                                                                                                          tCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGgatga    >  1:433284/1‑143 (MQ=255)

                                                                                                                          tCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAATTAGACCGCGATACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGgatga    >  1:667803/1‑143 (MQ=255)

                                                                                                                            cAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCCAGTATTTGACGTGCGGATACTTCGCCAGCACCGTCTCCAGTTCGCTCGCGTTACGCAAAATGTATAGATCGAGCCAAATAAAAACCGCAGGTAGCGgatgatg  >  1:89851/1‑143 (MQ=255)

                                                                                                                            cAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTATGATAATCAAGAAAACTACAACCCGCAGGTAGCGgatgatg  >  1:515388/1‑143 (MQ=255)

CGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATG  >  NC_000913/3176075‑3176341

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