BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,985,890	A→G	94.3%	N460D (AAC→GAC)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,985,890	0	A	G	94.3%	2158.5 / 80.7	629	N460D (AAC→GAC)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand): ref base A (5/31); new base G (136/457); total (141/488)
Fisher's exact test for biased strand distribution p-value = 3.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

Only 100 of 631 total aligned reads displayed.

CCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCT  >  CP009273/3985824‑3985957

  gCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGAcc                                                                  <  1:32845825/68‑1 (MQ=255)

  gCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGAcc                                                                  <  1:17271713/68‑1 (MQ=255)

  gCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGAcc                                                                  <  1:7358278/68‑1 (MQ=255)

    aGATTTCACCCTATCTCTCGGAACCGAATCTGACCTTTCTTTATGTGCCGCCGGGCCGGGCTGACCg                                                                 <  1:13937875/67‑1 (MQ=255)

       tttCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGt                                                                <  1:33224925/65‑1 (MQ=255)

       tttCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGt                                                                <  1:26719565/65‑1 (MQ=255)

       tttCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGt                                                                <  1:26088625/65‑1 (MQ=255)

       tttCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGt                                                                <  1:28225530/65‑1 (MQ=255)

       tttCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGt                                                                <  1:19194195/65‑1 (MQ=255)

       tttCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGt                                                                <  1:15030058/65‑1 (MQ=255)

       tttCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGt                                                                <  1:12851426/65‑1 (MQ=255)

         tCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAgg                                                           <  1:22104399/68‑1 (MQ=255)

         tCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAgg                                                           <  1:8834463/68‑1 (MQ=255)

         tCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAgg                                                           <  1:19664853/68‑1 (MQ=255)

          cACCCGATCTCTCGGAACCGTACCTGACCTTTATTTATGGGCCGCCGGGGCGGGCTGACCGTTCAGGt                                                          <  1:11879062/68‑1 (MQ=255)

          cACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGt                                                          <  1:23537460/68‑1 (MQ=255)

          cACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGt                                                          <  1:11397988/68‑1 (MQ=255)

          cACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGt                                                          <  1:34682237/68‑1 (MQ=255)

          cACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGt                                                          <  1:5790393/68‑1 (MQ=255)

            cccGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGCCCGTTCAGGTTg                                                        >  1:31663055/1‑68 (MQ=255)

            cccGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTg                                                        >  1:30650552/1‑68 (MQ=255)

            cccGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTg                                                        >  1:8454085/1‑68 (MQ=255)

            cccGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTg                                                        >  1:3652520/1‑68 (MQ=255)

             ccGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTgg                                                       <  1:1361256/68‑1 (MQ=255)

              cGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGt                                                      <  1:18778527/68‑1 (MQ=255)

              cGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGt                                                      <  1:16902724/68‑1 (MQ=255)

              cGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGt                                                      <  1:29903711/68‑1 (MQ=255)

              cGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGt                                                      <  1:4699514/68‑1 (MQ=255)

              cGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGt                                                      <  1:279227/68‑1 (MQ=255)

              cGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGt                                                      <  1:23768256/68‑1 (MQ=255)

               gATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTa                                                     <  1:20109759/68‑1 (MQ=255)

               gATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTa                                                     <  1:6172234/68‑1 (MQ=255)

                 tctctcGTAACCCAATCTGTCCTTTATTTATGTGGCGGCGGGACGGGCTGACCGTTCAGGTTGGTAt                                                    <  1:10892391/67‑1 (MQ=255)

                 tctctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTa                                                     >  1:24474088/1‑66 (MQ=255)

                  ctctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCTGGCTGACCGTTCAGGTTGGTATCt                                                  <  1:7018868/68‑1 (MQ=255)

                  ctctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGGCCGTTCAGGTTGGTATCt                                                  <  1:18381973/68‑1 (MQ=255)

                  ctctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCt                                                  <  1:8118797/68‑1 (MQ=255)

                  ctctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCt                                                  <  1:4995413/68‑1 (MQ=255)

                  ctctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCt                                                  <  1:25161417/68‑1 (MQ=255)

                  ctctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCt                                                  <  1:27290437/68‑1 (MQ=255)

                  ctctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCt                                                  <  1:28674332/68‑1 (MQ=255)

                  ctctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCt                                                  <  1:29092145/68‑1 (MQ=255)

                  ctctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCt                                                  <  1:30384702/68‑1 (MQ=255)

                  ctctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCt                                                  <  1:30942069/68‑1 (MQ=255)

                  ctctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCt                                                  <  1:31280608/68‑1 (MQ=255)

                  ctctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCt                                                  <  1:33972587/68‑1 (MQ=255)

                  ctctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCt                                                  <  1:14123615/68‑1 (MQ=255)

                  ctctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCt                                                  <  1:893966/68‑1 (MQ=255)

                  ctctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCt                                                  <  1:14467277/68‑1 (MQ=255)

                  ctctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCt                                                  <  1:229640/68‑1 (MQ=255)

                  ctctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCt                                                  <  1:10526127/68‑1 (MQ=255)

                  ctctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCt                                                  <  1:21663130/68‑1 (MQ=255)

                  ctctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCt                                                  <  1:10755370/68‑1 (MQ=255)

                  ctctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCt                                                  <  1:6442759/68‑1 (MQ=255)

                  ctctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCt                                                  <  1:19383766/68‑1 (MQ=255)

                  ctctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCt                                                  <  1:6719782/68‑1 (MQ=255)

                  ctctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCt                                                  <  1:17975261/68‑1 (MQ=255)

                  ctctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCt                                                  <  1:16320212/68‑1 (MQ=255)

                    ctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGt                                                <  1:24765599/68‑1 (MQ=255)

                    ctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGt                                                <  1:12540482/68‑1 (MQ=255)

                    ctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGt                                                <  1:15278956/68‑1 (MQ=255)

                    ctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGt                                                <  1:34582331/68‑1 (MQ=255)

                    ctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGt                                                <  1:318169/68‑1 (MQ=255)

                    ctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGt                                                <  1:2632970/68‑1 (MQ=255)

                    ctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGt                                                <  1:23991016/68‑1 (MQ=255)

                    ctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGt                                                <  1:20851380/68‑1 (MQ=255)

                      cGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGtat                                              >  1:5501584/1‑68 (MQ=255)

                               aTCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGTATAACcgcgcg                                     <  1:21277110/68‑1 (MQ=255)

                               aTCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGTATAACcgcgcg                                     <  1:5935129/68‑1 (MQ=255)

                                tCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGTATAACCgcgcgc                                    <  1:13392377/68‑1 (MQ=255)

                                tCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGTATAACCgcgcgc                                    <  1:9265879/68‑1 (MQ=255)

                                tCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGTATAACCgcgcgc                                    <  1:21464491/68‑1 (MQ=255)

                                tCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGTATAACCgcgcgc                                    <  1:33816725/68‑1 (MQ=255)

                                tCTGACCTTTATTTATGTGCCACCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGTATAACCgcgcgc                                    <  1:9442687/68‑1 (MQ=255)

                                 cTGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGTATAACCGCGCGcc                                   <  1:25830284/68‑1 (MQ=255)

                                  tGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCa                                  <  1:15842065/68‑1 (MQ=255)

                                  tGACCTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCa                                  <  1:22565581/68‑1 (MQ=255)

                                     ccTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCaaa                                <  1:7685796/67‑1 (MQ=255)

                                     ccTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAat                               <  1:30284907/68‑1 (MQ=255)

                                     ccTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAat                               <  1:29662436/68‑1 (MQ=255)

                                     ccTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAat                               <  1:28430087/68‑1 (MQ=255)

                                     ccTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAat                               <  1:27492139/68‑1 (MQ=255)

                                     ccTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAat                               <  1:27107978/68‑1 (MQ=255)

                                     ccTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAat                               >  1:25490320/1‑68 (MQ=255)

                                     ccTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAat                               <  1:20569435/68‑1 (MQ=255)

                                     ccTTTATTTATGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAat                               >  1:18156095/1‑68 (MQ=255)

                                            ttatGTGCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGt                        <  1:1775226/68‑1 (MQ=255)

                                            ttatGTGCCGCCGGGCCGGGCTGACCGTTCAGGCTGGTATCTGTATAACCGCGCGCCAAATATTGAg                         >  1:3262384/1‑67 (MQ=255)

                                               tgtgCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGa                     >  1:32069015/1‑68 (MQ=255)

                                               tgtgCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGa                     >  1:19794219/1‑68 (MQ=255)

                                               tgtgCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGa                     >  1:2423797/1‑68 (MQ=255)

                                               tgtgCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGa                     >  1:4287976/1‑68 (MQ=255)

                                               tgtgCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGa                     <  1:34280511/68‑1 (MQ=255)

                                                gtgCCGCCGGGCCGGGCTGTCCGTTCAGGGTGTTATCTGTATAACCGCGCTCCAAATATTGAGTCGat                    <  1:10310920/68‑1 (MQ=255)

                                                 tgCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGatc                   >  1:10087257/1‑68 (MQ=255)

                                                 tgCCGCCGGGCCGGGCTGACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGatc                   >  1:9578506/1‑68 (MQ=255)

                                                    cgccgGGCCGGGCTGACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGatcatc                <  1:3935171/68‑1 (MQ=255)

                                                       cgggccgggcTGACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATcatcagc             <  1:30780476/68‑1 (MQ=255)

                                                       cgggccgggcTGACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATcatcagc             <  1:23301191/68‑1 (MQ=255)

                                                       cgggccgggcTGACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATcatcagc             <  1:12246278/68‑1 (MQ=255)

CCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCT  >  CP009273/3985824‑3985957

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: `‑`