BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,985,309	T→C	23.4%	L266P (CTG→CCG)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,985,309	0	T	C	23.4%	1725.8 / inf	795	L266P (CTG→CCG)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand): ref base T (405/203); new base C (112/74); total (518/277)
Fisher's exact test for biased strand distribution p-value = 1.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 905 total aligned reads displayed.

TTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGT  >  CP009273/3985242‑3985376

 ttGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACAct                                                                     >  1:17433289/1‑67 (MQ=255)

  tGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACActgc                                                                   <  1:629613/68‑1 (MQ=255)

   ggTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACAccgct                                                                  <  1:9129107/68‑1 (MQ=255)

     tGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACActgctg                                                                 <  1:2517178/67‑1 (MQ=255)

     tGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACActgctg                                                                 <  1:23094877/67‑1 (MQ=255)

     tGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACActgct                                                                  >  1:5381525/1‑66 (MQ=255)

     tGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACActgct                                                                  >  1:19817870/1‑66 (MQ=255)

      gCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATAGAAAGCGGTACTGAAAACActgctgct                                                               <  1:23717956/68‑1 (MQ=255)

      gCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACActgct                                                                  >  1:27336084/1‑65 (MQ=255)

      gCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACActgct                                                                  >  1:28049700/1‑65 (MQ=255)

      gCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACActgct                                                                  >  1:22096749/1‑65 (MQ=255)

      gCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACAccgctgc                                                                <  1:9259185/67‑1 (MQ=255)

      gCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACAccgct                                                                  >  1:31144813/1‑65 (MQ=255)

       ccAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACActgctgct                                                               <  1:32603081/67‑1 (MQ=255)

            cTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAg                                                          <  1:4402200/67‑1 (MQ=255)

            cTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAg                                                          <  1:557871/67‑1 (MQ=255)

            cTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGc                                                         <  1:3769382/68‑1 (MQ=255)

                ggCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTAt                                                     <  1:5876125/68‑1 (MQ=255)

                ggCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTAt                                                     <  1:30381699/68‑1 (MQ=255)

                ggCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTAt                                                     <  1:16773617/68‑1 (MQ=255)

                ggCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTAt                                                     <  1:16418486/68‑1 (MQ=255)

                ggCAGCCCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTAt                                                     <  1:21201930/68‑1 (MQ=255)

                 gCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTAt                                                     >  1:1782165/1‑67 (MQ=255)

                 gCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTAt                                                     >  1:12383185/1‑67 (MQ=255)

                 gCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTAt                                                     >  1:29145905/1‑67 (MQ=255)

                  cAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTAtt                                                    >  1:18139631/1‑67 (MQ=255)

                  cAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTAtt                                                    <  1:8419845/67‑1 (MQ=255)

                  cAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTAtt                                                    <  1:24164962/67‑1 (MQ=255)

                  cAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTAtt                                                    <  1:12743931/67‑1 (MQ=255)

                      tctACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTcc                                                  <  1:8117360/65‑1 (MQ=255)

                      tctACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTcc                                                  <  1:20568923/65‑1 (MQ=255)

                      tctACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTgg                                               <  1:24929705/68‑1 (MQ=255)

                           aaGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAAt                                           >  1:10691293/1‑67 (MQ=255)

                           aaGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAAt                                           >  1:29651620/1‑67 (MQ=255)

                            agagTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATa                                          >  1:33523120/1‑67 (MQ=255)

                            agagTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATa                                          >  1:14101565/1‑67 (MQ=255)

                            agagTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATa                                          >  1:13093860/1‑67 (MQ=255)

                            agagTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATa                                          >  1:31714138/1‑67 (MQ=255)

                            agagTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATa                                          >  1:20986604/1‑67 (MQ=255)

                            agagTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATa                                          >  1:14270341/1‑67 (MQ=255)

                            agagTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTAGTCCTGGGAATa                                          >  1:22454943/1‑67 (MQ=255)

                            agagTATCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATAc                                         <  1:33822123/68‑1 (MQ=255)

                              agTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATAccc                                       >  1:21591304/1‑68 (MQ=255)

                              agTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATAccc                                       >  1:18928699/1‑68 (MQ=255)

                              agTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATAccc                                       >  1:6920644/1‑68 (MQ=255)

                              agTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATAccc                                       >  1:160470/1‑68 (MQ=255)

                              agTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATAccc                                       >  1:10005191/1‑68 (MQ=255)

                              agTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTCGGAATAccc                                       >  1:26558933/1‑68 (MQ=255)

                              agTATCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATAccc                                       >  1:3440082/1‑68 (MQ=255)

                              agTATCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATAccc                                       >  1:1164148/1‑68 (MQ=255)

                                 aTCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATACCCGa                                     <  1:20743562/67‑1 (MQ=255)

                                  tCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATACCCGAAc                                   <  1:5048038/68‑1 (MQ=255)

                                   cGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATACCCGAAc                                   >  1:777092/1‑67 (MQ=255)

                                    gATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATACCCGAAccc                                 >  1:9733555/1‑68 (MQ=255)

                                     aTTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAAccc                                 >  1:11936994/1‑67 (MQ=255)

                                      ttCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCa                                <  1:20208171/67‑1 (MQ=255)

                                      ttCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTACTGGGAATACCCGAACCCAc                               >  1:30585322/1‑68 (MQ=255)

                                       tCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACg                              >  1:7986268/1‑68 (MQ=255)

                                        ccccATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACg                              >  1:8754459/1‑67 (MQ=255)

                                          ccATACAAAGCGGTACTGAAAACACTGCTGCTGGAAg                                                          <  1:18705605/37‑1 (MQ=255)

                                              aCAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTctgctg                       >  1:25561168/1‑68 (MQ=255)

                                              aCAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTctgctg                       >  1:24652345/1‑68 (MQ=255)

                                              aCAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTctgctg                       >  1:23384827/1‑68 (MQ=255)

                                              aCAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTctgctg                       >  1:1486766/1‑68 (MQ=255)

                                              aCAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTctgctg                       >  1:1096437/1‑68 (MQ=255)

                                              aCAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTctgctg                       >  1:31510513/1‑68 (MQ=255)

                                              aCAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTctactg                       >  1:4223336/1‑68 (MQ=255)

                                              aCAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTctgctg                       <  1:7393375/68‑1 (MQ=255)

                                               cAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTctgctg                       >  1:30158522/1‑67 (MQ=255)

                                               cAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTgg                      <  1:25885525/68‑1 (MQ=255)

                                               cAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTgg                      <  1:28791132/68‑1 (MQ=255)

                                               cAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTctgctg                       >  1:26848544/1‑67 (MQ=255)

                                               cAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTctgctg                       >  1:32831901/1‑67 (MQ=255)

                                               cAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTgg                      <  1:21904973/68‑1 (MQ=255)

                                                aaaGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGc                     >  1:22681282/1‑68 (MQ=255)

                                                   gCGGTACTGAAAACACTGCTGCTGGAAGCTTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGaa                  >  1:27826665/1‑68 (MQ=255)

                                                   gCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGCATACCCGAACCCACGTCTGCTGGCGaa                  >  1:13658207/1‑68 (MQ=255)

                                                   gCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGaa                  >  1:34180599/1‑68 (MQ=255)

                                                   gCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGaa                  >  1:15758169/1‑68 (MQ=255)

                                                   gCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGaa                  >  1:13507855/1‑68 (MQ=255)

                                                   gCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACACGAACCCACGTCTGCTGGCGaa                  >  1:18321035/1‑68 (MQ=255)

                                                   gCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGaa                  >  1:31979036/1‑68 (MQ=255)

                                                        aCTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGata             >  1:4722652/1‑68 (MQ=255)

                                                        aCTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGat              <  1:27091692/67‑1 (MQ=255)

                                                        aCTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGat              <  1:19204251/67‑1 (MQ=255)

                                                        aCTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGat              <  1:3085750/67‑1 (MQ=255)

                                                        aCTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGata             >  1:10353775/1‑68 (MQ=255)

                                                         cTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGatat            <  1:154649/68‑1 (MQ=255)

                                                           gAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCa          <  1:14458038/68‑1 (MQ=255)

                                                             aaaCACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCaa         >  1:25381220/1‑67 (MQ=255)

                                                             aaaCACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCaa         >  1:15239712/1‑67 (MQ=255)

                                                             aaaCACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCaa         >  1:32153414/1‑67 (MQ=255)

                                                             aaaCACCGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCaa         >  1:26150324/1‑67 (MQ=255)

                                                               acacTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACa      >  1:28319854/1‑68 (MQ=255)

                                                               acacTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACa      >  1:29308239/1‑68 (MQ=255)

                                                               acacTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACa      >  1:11240007/1‑68 (MQ=255)

                                                               acacCGCTGCTGGAAGCCTATTCCTGGTAATACACGAACCCAATTCTGCTTGCGAAAGATATCAAACa      >  1:7211750/1‑68 (MQ=255)

                                                               acacCGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACa      <  1:17086153/68‑1 (MQ=255)

                                                                 acTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGc    <  1:29953326/68‑1 (MQ=255)

                                                                  ctgctgctgGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCg   <  1:33158632/68‑1 (MQ=255)

TTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGT  >  CP009273/3985242‑3985376

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