BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,900,410	A→T	9.6%	P56P (CCT→CCA)	yieL ←	putative hydrolase YieL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,900,410	0	A	T	9.6%	95.2 / 2.5	40	P56P (CCT→CCA)	yieL	putative hydrolase YieL
Reads supporting (aligned to +/- strand): ref base A (15/20); new base T (2/2); total (17/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.04e-01

CAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGA‑TAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAG  >  NC_000913/3900271‑3900550

caccATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGA‑TAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAggg                                                                                                                                           <  1:8048/140‑1 (MQ=255)

  aCATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGA‑TAGGTGTGTGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCTGGGGc                                                                                                                                         <  2:342923/143‑1 (MQ=255)

     tgtgaaaaaTACTCGTACAGGTTGCCTTTCAGGA‑TAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGTGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCaaa                                                                                                                                      <  2:492372/139‑1 (MQ=255)

       gAAAAAATACTCGTACAGGTTGCCTTTCAGGA‑TAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGt                                                                                                                                    <  1:244563/143‑1 (MQ=255)

       gAAAAAATACTCGTACAGGTTGCCTTTCAGGA‑TAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTCGTCATCGGGTGAGTCTTGTCCGGAACCGGAACACCCCCCACTACAGAGACATTTTTTGCCCCAGGGGCaaact                                                                                                                                    >  2:291338/1‑141 (MQ=255)

        aaaaacTACTCGTCCATGTTGCCTTTCAGGA‑TAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATGTTTTGCCCCAGGGGCAAAGTa                                                                                                                                   <  1:134865/143‑1 (MQ=255)

        aaaaaaTACTCGTACAGGTTGCCTTTCAGGA‑TAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTa                                                                                                                                   >  1:561133/1‑143 (MQ=255)

                 tCGTACAGGTTGCCTTTCAGGA‑TAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGCAAGGt                                                                                                                          >  1:371662/1‑143 (MQ=255)

                   gTACAGGTTGCCTTTCAGGA‑TAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCa                                                                                                                        <  1:289171/143‑1 (MQ=255)

                      cAGGTTGCCGTTCAGGA‑TAGTGGTGCTCCACGACCAGCCTCCTGCTTCGTCTTTGGGCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGc                                                                                                                     <  2:145031/143‑1 (MQ=255)

                       aGGTTGCCTTTCAGGA‑TAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCt                                                                                                                    <  1:511396/143‑1 (MQ=255)

                             cttgtCAGGATTAGGTG‑GCGCCACGAACAGATTCCTGCTTCGTCGTTGGTCATCGGGAGAATATTGTCCGGAACCGGAATACCCACTACTACAGAGACATTTTTTGCCCCTGGGGTAAAGTAGCGAAAGGTGACGCTGTTATc                                                                                                              <  2:44469/139‑1 (MQ=255)

                                ttCAGGA‑TAGGTGTGGGGTGTGCCCAGACTCCTACTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGc                                                                                                           <  2:155076/143‑1 (MQ=255)

                                ttCAGGA‑TAGGTGTGCGCCACGACCAGACTCCTGCTTCGACTTGGGTCATCGGGTGGATATTGTCCGGACCCGGAACACCCACCACTACAGAGACATTTTTTGCGCCAGGGGCAACGTAGCGCCAGGTCGCGCTGTTATCGGc                                                                                                           >  1:553934/1‑143 (MQ=255)

                                    tggggaGGTGTGCGGCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACACTGTTGTCGGCAtt                                                                                                        <  2:430390/138‑1 (MQ=255)

                                      gatgtgtgtGCGCCACGACCAGACTCCTGCTTCGAATTGCCTTAGCGTGTGAATCTTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGa                                                                                                      <  2:283092/138‑1 (MQ=255)

                                                ctcgaCGACAAGACTCCTGCTTCGTCTGTAGTCATCCGCTGAATACTGTCCGGACCCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTcac                                                                                            <  1:291338/139‑1 (MQ=255)

                                                 gcCACGACCCGCCGCCTGCTTCGTCTTTGGTCATCCGGTGAATATTGTCCGCAACCGGAACACCCACCACTACAGAGACATTTTTTTCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTcaca                                                                                           <  2:443462/143‑1 (MQ=255)

                                                     cGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGCCATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACt                                                                                       >  1:405948/1‑143 (MQ=255)

                                                       aCCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGt                                                                                     <  1:465729/143‑1 (MQ=255)

                                                         cAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCGTTGACGTGAGTCAGGTGCTGttt                                                                                   >  1:80943/1‑143 (MQ=255)

                                                         cAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACAGACTGttt                                                                                   >  1:102462/1‑143 (MQ=255)

                                                            tcTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCAc                                                                                <  2:570171/142‑1 (MQ=255)

                                                            aCTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTAACGCCGTTATCGGCATTGACCTGAGTCACATACTGGTTCAc                                                                                >  1:236767/1‑143 (MQ=255)

                                                              tCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGGCACATACTGTTTCACCg                                                                              >  1:399924/1‑143 (MQ=255)

                                                               ccTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCgg                                                                             <  2:33545/143‑1 (MQ=255)

                                                                       gTCTTTGGTCCTCGGGTGAATATGGTCCGGAACCGCAACACCCACCACTACAGAGACATTTTTTGCCCCAGCGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTgg                                                                     <  1:369572/143‑1 (MQ=255)

                                                                              gTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACCACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGACAGTGGGGGCCgg                                                              >  2:564292/1‑143 (MQ=255)

                                                                               tCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCggg                                                             <  2:435465/143‑1 (MQ=255)

                                                                                    gggTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACCTACTGTTTCACCGGGATAGTGGGTGCTGggctggc                                                        >  1:375487/1‑143 (MQ=255)

                                                                                        gAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCtggctggc                                                    <  2:272553/143‑1 (MQ=255)

                                                                                                                ccaccaCTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGa                            >  1:35499/1‑143 (MQ=255)

                                                                                                                 caccaccacAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAg                           <  1:410898/143‑1 (MQ=255)

                                                                                                                   ccacTACAGACACATTTCTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCa                         <  1:80306/143‑1 (MQ=255)

                                                                                                                            ataCATTTTTTGCCCCCGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGTTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCg                <  1:254674/141‑1 (MQ=255)

                                                                                                                             gaCATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGc               >  1:319874/1‑143 (MQ=255)

                                                                                                                                ctactttGCCCCTGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCGCACTGATCAGAAATAGCGCTgg            >  2:443538/5‑143 (MQ=255)

                                                                                                                                  tttttGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGCAATACCGCTGGCa          >  1:397416/1‑143 (MQ=255)

                                                                                                                                   ttttGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAATGGATAATGTTCTTGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATAACGCCGGCaa         >  2:224969/1‑143 (MQ=255)

                                                                                                                                        ccccTGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCAAATCCGCTGCTATGGCCCACGGAGCAGAAATAACGCGGGCAAGAGcc    >  2:159312/1‑143 (MQ=255)

                                                                                                                                          ctaGGGGCAAAGTACCGAAAGGTCACGCTGTCATCGGCATTGACCTGAGTCACATACGGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAg  <  2:80943/141‑1 (MQ=255)

CAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGA‑TAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAG  >  NC_000913/3900271‑3900550

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