BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,421,797	C→T	7.4%	T220T (ACG→ACA)	ccmF ←	heme lyase, CcmF subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,421,797	0	C	T	7.4%	441.9 / 18.0	135	T220T (ACG→ACA)	ccmF	heme lyase, CcmF subunit
Reads supporting (aligned to +/- strand): ref base C (77/48); new base T (5/5); total (82/53)
Fisher's exact test for biased strand distribution p-value = 5.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

Only 100 of 135 total aligned reads displayed.

CCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCA  >  Exported/1421731‑1421858

ccAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGt                                                             >  1:3417972/1‑69 (MQ=255)

ccAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCa                                                           <  1:5918589/71‑1 (MQ=255)

ccAGAACCACCAGCCA‑CCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCa                                                           <  1:4195485/70‑1 (MQ=255)

   gAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAgg                                                         >  1:2460854/1‑70 (MQ=255)

   gAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGa                                                        <  1:6850390/71‑1 (MQ=255)

   gAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGa                                                        <  1:4205686/71‑1 (MQ=255)

   gAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGa                                                        <  1:6731566/71‑1 (MQ=255)

   gAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGa                                                        <  1:2226269/71‑1 (MQ=255)

   gAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGa                                                        <  1:2057082/71‑1 (MQ=255)

     accaccAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAg                                                      <  1:2733214/71‑1 (MQ=255)

     accaccAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAg                                                      <  1:2086605/71‑1 (MQ=255)

      ccaccaGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAgg                                                         >  1:670177/1‑67 (MQ=255)

      ccaccaGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGa                                                     <  1:4717489/71‑1 (MQ=255)

                 cccAGCCGAGTTCGTAATAGTCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATc                                                   <  1:5747492/62‑1 (MQ=255)

                  ccAGCCGAGTTCGTAATAGGCCCATGCGTAACCGAGCACGATGCCGAGCGTCAGGAAGATCCAc                                                >  1:3529048/1‑64 (MQ=255)

                  ccAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCAc                                                >  1:1600313/1‑64 (MQ=255)

                  ccAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCAc                                                >  1:882871/1‑64 (MQ=255)

                  ccAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCAc                                                >  1:2040681/1‑64 (MQ=255)

                  ccAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCAc                                                >  1:4101962/1‑64 (MQ=255)

                   cAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGt                                                             >  1:1581536/1‑50 (MQ=255)

                   cAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGt                                                             >  1:3278688/1‑50 (MQ=255)

                   cAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGt                                                             >  1:1831228/1‑50 (MQ=255)

                   cAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGt                                                             >  1:1252698/1‑50 (MQ=255)

                   cAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGt                                                             >  1:1273539/1‑50 (MQ=255)

                   cAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCa                                                 >  1:2349792/1‑62 (MQ=255)

                   cAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAg                                        >  1:4205865/1‑71 (MQ=255)

                   cAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAg                                        >  1:2825943/1‑71 (MQ=255)

                     gCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCg                                      >  1:5598507/1‑71 (MQ=255)

                         aGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGTGTCAGGa                                                        <  1:1283050/49‑1 (MQ=255)

                         aGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGTGTCAGGa                                                        <  1:6937982/49‑1 (MQ=255)

                         aGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGTGTCAGGa                                                        <  1:2814577/49‑1 (MQ=255)

                         aGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCa                                  >  1:3711853/1‑71 (MQ=255)

                           ttCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGTGTCAGGAAg                                                      >  1:3455199/1‑49 (MQ=255)

                           ttCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGTGTCAGGAAg                                                      >  1:3035145/1‑49 (MQ=255)

                           ttCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGTGTCAGGAAg                                                      >  1:1507828/1‑49 (MQ=255)

                           ttCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCa                                         <  1:5790960/62‑1 (MQ=255)

                           ttCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCg                                      <  1:4680722/65‑1 (MQ=255)

                           ttCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACg                                >  1:1694182/1‑71 (MQ=255)

                           ttCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACg                                >  1:5108428/1‑71 (MQ=255)

                           ttCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTCCCAGCg                                      <  1:5437860/65‑1 (MQ=255)

                           ttCGTAATAGGCCCATGAGGAACCGAGCACGATGCCGAGTGTCAGGAAg                                                      >  1:4573169/1‑49 (MQ=38)

                             cGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCAc                                                >  1:602411/1‑53 (MQ=255)

                             cGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCAc                                                >  1:2494623/1‑53 (MQ=255)

                             cGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCAc                                                >  1:5172222/1‑53 (MQ=255)

                             cGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGc                                       >  1:2796551/1‑62 (MQ=255)

                             cGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTc                                    >  1:2626860/1‑65 (MQ=255)

                             cGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGa                              >  1:4380880/1‑71 (MQ=255)

                             cGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGa                              >  1:3076102/1‑71 (MQ=255)

                               tAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCAc                                 >  1:5528530/1‑66 (MQ=255)

                               tAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCAc                                 >  1:4861134/1‑66 (MQ=255)

                               tAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCAc                                 >  1:5469508/1‑66 (MQ=255)

                                 aTAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAg                          >  1:3271665/1‑71 (MQ=255)

                                  tAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATGCACGCTGCCAGCGTCCACGGACGAGt                         <  1:3902424/71‑1 (MQ=255)

                                   aGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGt                         >  1:5633189/1‑70 (MQ=255)

                                    ggCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTaa                       <  1:3137387/71‑1 (MQ=255)

                                    ggCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTaa                       <  1:6716887/71‑1 (MQ=255)

                                     gCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTaaa                      <  1:5713825/71‑1 (MQ=255)

                                         aTGCGGAACCGAGCACGATGCCGAGCGTCAGGAAg                                                      <  1:3557765/35‑1 (MQ=255)

                                         aTGCGGAACCGAGCACGATGCCGAGCGTCAGGAAg                                                      <  1:1609694/35‑1 (MQ=255)

                                         aTGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTaa                       <  1:639753/66‑1 (MQ=255)

                                         aTGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTaa                       <  1:486316/66‑1 (MQ=255)

                                         aTGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTaa                       <  1:4313344/66‑1 (MQ=255)

                                         aTGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTaa                       <  1:3133229/66‑1 (MQ=255)

                                         aTGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAAcg                   >  1:6597475/1‑70 (MQ=255)

                                           gCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATc                                                   >  1:5750026/1‑36 (MQ=255)

                                           gCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATc                                                   >  1:4571980/1‑36 (MQ=255)

                                           gCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATc                                                   >  1:2590195/1‑36 (MQ=255)

                                           gCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTcc                                   >  1:3095764/1‑52 (MQ=255)

                                           gCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCa                                  >  1:6025410/1‑53 (MQ=255)

                                           gCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTa                        <  1:1513336/63‑1 (MQ=255)

                                           gCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTa                        <  1:1458654/63‑1 (MQ=255)

                                            cGGATCCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTaaaa                     >  1:1005477/1‑65 (MQ=255)

                                            cGGAACCGAGCACGATGCCGAGCGTCAGGAAGATc                                                   >  1:909314/1‑35 (MQ=255)

                                            cGGAACCGAGCACGATGCCGAGCGTCAGGAAGATc                                                   >  1:1205258/1‑35 (MQ=255)

                                            cGGAACCGAGCACGATGCCGAGCGTCAGGAAGATc                                                   >  1:3388038/1‑35 (MQ=255)

                                            cGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAg                          >  1:6120918/1‑60 (MQ=255)

                                            cGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTaaaa                     >  1:5624506/1‑65 (MQ=255)

                                            cGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTaaaa                     >  1:6902549/1‑65 (MQ=255)

                                             ggAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCAt              <  1:2459768/71‑1 (MQ=255)

                                              gAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATa             <  1:3267119/71‑1 (MQ=255)

                                              gAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATa             <  1:1179683/71‑1 (MQ=255)

                                                 ccGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAAcg                   >  1:5863977/1‑62 (MQ=255)

                                                  cGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTc                                    <  1:2868129/44‑1 (MQ=255)

                                                  cGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTc                                    <  1:4534374/44‑1 (MQ=255)

                                                  cGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTc                                    <  1:5713309/44‑1 (MQ=255)

                                                  cGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTc                                    <  1:3675580/44‑1 (MQ=255)

                                                  cGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTc                                    <  1:290917/44‑1 (MQ=255)

                                                  cGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTc                                    <  1:5815387/44‑1 (MQ=255)

                                                   gAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGc        >  1:2700982/1‑71 (MQ=255)

                                                    aGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCa                                  >  1:4075161/1‑44 (MQ=255)

                                                    aGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCa                                  >  1:1764635/1‑44 (MQ=255)

                                                    aGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCt       <  1:2431468/71‑1 (MQ=255)

                                                      cACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGt          >  1:4450653/1‑66 (MQ=255)

                                                        cGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCt       <  1:2691689/67‑1 (MQ=255)

                                                         gATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCa  <  1:3549807/71‑1 (MQ=255)

                                                           tGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCa               >  1:5513873/1‑56 (MQ=255)

                                                                aGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGa                              >  1:6071406/1‑36 (MQ=255)

                                                                aGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGa                              >  1:3732204/1‑36 (MQ=255)

                                                                aGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGa                              >  1:2784920/1‑36 (MQ=255)

                                                                aGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGa                              >  1:4457129/1‑36 (MQ=255)

CCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCA  >  Exported/1421731‑1421858

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